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TECHNICAL MEMORANDUM 


TRACE CONTAMINANT CONTROL SIMULATION COMPUTER 
PROGRAM— VERSION 8.1 


INTRODUCTION 


The Trace Contaminant Control Simulation (TCCS) computer program development began with 
the efforts of Olcott which were documented in 1972. Since then, improvements in the user interface and 
more up-to-date information on activated charcoal loading characteristics and high temperature oxida- 
tion catalysts have become available. The progression of the program from version 1.0 through 8.0 and 
acknowledgments for its development are documented by reference 1. The descriptive material on the 
computer program subroutines has been extracted in its entirety from reference 1; however, some modi- 
fications have been made where necessary. 


VERSION 8.0 MODIFICATIONS 


The TCCS computer program version 8.0 which is documented in reference 1 has been modified 
as part of the International Space Station (ISS) design process. These modifications were made initially 
by personnel at Lockheed Missiles and Space Co. in Sunnyvale, CA, as part of the design activities for 
the ISS trace contaminant control subassembly. During these activities, new information on activated 
charcoal loading capacity, influence of relative humidity on charcoal loading capacity, and preliminary 
information on poisoning characteristics of high temperature oxidation catalysts was obtained. This 
information was assessed by both Lockheed and NASA Marshall Space Flight Center personnel and 
integrated into the TCCS computer program as version 8.1. These modifications were officially accepted 
by NASA in March 1994. Specific details on the modifications can be found in reference 2. A listing of 
the source files for version 8.1 is contained in appendix A. Version 8.1 supersedes all other versions of 
the program. 

A brief summary of changes to version 8.0 which have resulted in the new version 8.1 are the 
following: 

1. The main program, MAIN. FOR, was modified in the following ways: 

(a) Matrix TT was increased from 300 by 7 to 750 by 7 to accommodate larger time- 
dependent data files 

(b) Code which sets the time increment size was moved to precede the code which checks 
for changes in the basic time increment 

(c) The code was modified to accept changes in the basic time increment size as long as the 
change occurs between the beginning and end of the current time step rather than at the 
beginning of the time step only 



(d) Changes to the cabin volume now result in a recalculation of all contaminant concen- 
trations 

(e) Mission duration output formats have been changed to accommodate number sizes up 
to five digits. 

2. Testing of Bamebey-Sutcliffe types AC and 3032 activated carbon with and without 10 
weight percent phosphoric acid impregnation has resulted in new charcoal capacity equa- 
tions. The programs ACHBD.FOR and RCHBD.FOR were modified in the following ways: 

(a) Capacity for water soluble contaminants, those with a Henry’s Law constant between 0 
and 5, was determined to be a function of the adsorption potential factor, A, only and 
not a function of relative humidity. Activated carbon capacity for this case is expressed 
by the following equations: 

q- 2.1<r°- 3M for A > 8 , (1) 

q = 0.5-0.0405 A for A < 8 . (2) 

(b) Insoluble contaminants were found to be a function of adsorption potential factor. A, 
and relative humidity, Hr. Activated carbon capacity for this case is expressed by the 
following equations: 

q = (-1.28x10“ 6 )//^-(2.64x10- 3 * )//^+0.5+[(1.12x10- 6 )//|+(2.08x 10^^-0.0405]A for A < 8 , (3) 
q = [(-9.6xlO _5 )//^-(1.88xlO _2 )//j ? -2.11]e -03M for A >8 \H R < 50 percent , (4) 

^ = [(9.6xl0- 5 )//^-(1.88xl0- 2 )// J? -2.11]H°-25+o.ooi2W*vt for A >% ; Hr > 50 percent . (5) 

(c) The programs were modified to read cabin percent relative humidity from the device 
definition matrix, DD, row l column 14 

(d) Carbon chemisorption capacity for ammonia at its spacecraft maximum allowable con- 
centration (SMAC) was changed to 0.0061 grams of ammonia per gram of carbon to 
reflect the latest phosphoric acid impregnated charcoal performance data. The previous 
number was based on theoretical estimates rather than experimental results. 

3. Subprogram CATBNR.FOR was modified to reflect poisoning of 0.5-percent palladium on 
alumina catalyst by halocarbons and sulfide compounds. The result of this poisoning is a 
decrease in the removal efficiency, 7), for methane. The efficiency calculation was modified 
to account for the total mass of halocarbon and sulfide compounds, P, in milligrams by using 
the following equations: 

q = 97.506x 10-° 00010507 p for P< 5,500 mg , (6) 

77 = 31.453-0. 151xl0- 3 )P+(1.9045xl0- 8 )/> 2 -(1.0389xl0- 13 )/> 3 for P> 5,500 mg . (7) 
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4. The following changes were made to subprogram CNRSUB.FOR: 

(a) Modifications were made to allow transferring the relative humidity value to 
ACHBD.FOR and RCHBD.FOR 

(b) Partial catalytic oxidizer efficiency restoration for methane removal is set to coincide 
with axial and radial charcoal bed regeneration 

(c) Cumulative masses of halocarbons and sulfides removed by the catalytic oxidizer are 
transferred to CATBNR.FOR. The cumulative mass is reinitialized at charcoal bed 
regeneration. 

(d) Code was added to allow for reinitialization of methane oxidation efficiency in 
CATBNR.FOR if an upstream adsorption device is regenerated. 


DETAILED COMPUTER PROGRAM DESCRIPTION 


A detailed description of the TCCS computer program source files is provided to acquaint the 
user with the main program and each significant subroutine. Flow charts of these routines are provided 
and discussion of the theoretical basis for some routines is provided where appropriate. A listing of the 
program source files is provided in appendix A. This description is paraphrased from a description pro- 
duced by Lockheed Missiles and Space Co, Inc., under contract NAS8-36406. This work served as the 
primary reference for this section, and all block flow diagrams were adapted from this document. 3 


Program Editing, Compiling, and Linking 

This program was edited, compiled, and linked using the Ryan-McFarland RM/FORTRAN™ 
version 2.42 which include the RM/FORTE™ project manager. This FORTRAN compiler is recom- 
mended for making changes to the source files. 


Main Program 

The main program, MAIN, is a simple program with no branching and two loops. A flow dia- 
gram is shown in figure 1. Each subroutine required for the particular program run is called during each 
pass of the main calculation loop until the end of the simulation. 

Subroutines CAFILL and RAFILL, which write zeros into all the calculation matrices are called 
initially to initialize each calculation matrix. Next, CRIN and PRIN are called to read the contaminant, 
device definition, and time-dependent input data into matrices NN, CDI, DD, and TT. The input data are 
printed line by line, if desired, by calling subroutines CROUT2 and RROUT2. All initial variables such 
as time increment beginning time, time increment ending time, and the increment counter are zeroed. 

The precalculation setup routine, PCSET, is called next. This routine calculates the initial 
removal efficiency for each removal device, the equilibrium cabin concentration, and the final cabin 
concentration for a cabin concentration of lxlO -20 mg/m 3 for all contaminants. Intermediate and final 
calculation results are stored in matrices CC and DD. 
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Figure 1. Main program block flow diagram. 


The calculation loop is now entered, and the iterative process of determining the cabin concen- 
tration for each contaminant at the end of a time increment is begun. Since cabin concentration is a 
function of the contaminant mass removed and the contaminant mass removed is a function of the cabin 
concentration, it is important that the same value be used in the mass balance and removal efficiency 
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calculation routines. A solution is achieved by assuming an increment concentration, calculating an 
average contaminant concentration, and comparing the two concentrations. This procedure continues 
until the assumed and calculated concentrations are equal. 

Calculation for each time increment is initiated by increasing the increment counter, setting the 
increment size, and listing the increment number, starting time, and ending time to the computer screen. 
Subroutine RINCDD is called to read time-dependent data from matrix TT at this time, and subroutine 
REGEN is called to check for regeneration of any devices during the time increment and to set the 
adsorbed contaminant masses and device flow rates to zero if necessary. The main calculation loop sub- 
routine, MCALC, is called next to calculate the removal efficiencies, average calculated concentration, 
and final concentration for each contaminant based on the sum of the mass removed during the previous 
time increment. Subroutine SLIOH then calculates the amount of lithium hydroxide (LiOH) used during 
the time increment if a LiOH bed is specified in the device definition file. Next, the original device flow 
rate is restored for any device that was being regenerated during the increment. Subroutine DATOUT is 
called to print the calculated data at the end of a time increment, if necessary, to both the standard for- 
matted and plot data output devices. The simulation and mission duration times are then compared to 
determine whether the mission simulation has ended. If the mission simulation has not ended, another 
pass through the calculation loop begins by setting the new increment beginning time and initial cabin 
concentration equal to the previous increment ending time and concentration. If the mission simulation 
has ended, subroutine DATOUT is called to write the final answers to the appropriate output devices 
as specified by the user. The output files are closed, and the program loops to the beginning to begin 
another run if the user wishes. If the user has no other runs to make, the program execution is termi- 
nated, otherwise, the calculation matrices are zeroed and new input data is supplied to the program for 
the next run. 

Brief descriptions for each major TCCS computer program subroutine are provided in the order 
that they are called by MAIN. Table 1 lists the subroutines as they are called and provides a brief 
description each subroutine’s purpose. Block flow diagrams are provided for the most significant sub- 
routines. 


Calculation Loop Subroutines 

The following subroutines comprise the principal calculation framework for the TCCS Computer 
Program. 

Subroutine CAFILL 

The subroutine CAFILL is called by MAIN and fills the matrix NN with blanks. Matrix NN 
contains the contaminant names during the simulation run. 

Subroutine R A FILL 

The subroutine RAFILL is called by MAIN and fills the matrices CC, TT, CDI, and DD with 
zeros. Matrices CDI and CC contains contaminant input and calculation data, matrix DD contains device 
calculation data, and matrix TT contains time-dependent data. This routine is used at the beginning of a 
computer simulation to initialize these matrices in the event a previous run has been made. 
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Table 1. TCCS computer program subroutine listing and description. 


1 

MAIN 


Subrouting Level 


2 3 


SLIOH 

DATOUT 


PRAFIL 

CNRSUB 


MASBAL 


PREDCT 


CONVRG 


PRFANS 


GROUP 


ACHBD 

RCHBD 

ALIOH 

COOXID 

CATBNR 

CONDHX 

CALCM 

LDIGEN 

PCAVCF 


PRAFIL 

LODEFF 

MASBAL 


PRAFIL 

CNRSUB 


MASBAL 


HEADGS 


CALCM 

LDIGEN 

PCAVCF 


ACHBD 

RCHBD 

ALIOH 

COOXID 

CATBNR 

CONDHX 

CALCM 

LDIGEN 

PCAVCF 


Description 


MAIN PROGRAM 
ZERO MATRIX NN 
ZERO MATRICES CDI, CC, DD, TT 
INPUT DATA INTO NN AND CDI 
INPUT FROM FILE TO DD AND TT 
PRINT DATA FROM NN AND CDI 
PRINT DATA FROM DD AND TT 
PRECALCULATION SETUP FOR ALL CONT 
ZERO MATRIX DD COLUMNS 17-21 
CALCULATE REMOVAL EFFICIENCY 
AXIAL CHARCOAL BED EFFICIENCY 
RADIAL CHARCOAL BED EFFICIENCY 
AXIAL LiOH BED EFFICIENCY 
CO OXIDIZER EFFICIENCY 
CATALYTIC OXIDIZER EFFICIENCY 
CONDENSATE EFFICIENCY 
CONTAMINANT MATERIAL BALANCE 
SUM OF MASS REMOVED BY DEVICES 
LOAD GENERATION INTO DD COL 19 
CALCULATE FINAL AND AVERAGE CONC 
READ INCREMENT DEPENDENT DATA 
CALCULATE REGENERATION/CHANGEOUT 
MAIN CALCULATION ROUTINE 
CALCULATE PREDICTED AVERAGE CONC 
ZERO MATRIX DD COLUMNS 17-21 
PUT LAST INCREMENT EFFICIENCY IN DD 
MASS BALANCE ROUTINE 
SUM OF MASS REMOVED BY DEVICES 
LOAD GENERATION INTO DD COL 19 
CALCULATE FINAL AND AVERAGE CONC 
SOLVE FOR NEW REMOVAL EFFICIENCY 
ZERO MATRIX DD COLUMNS 17-21 
CALCULATE REMOVAL EFFICIENCIES 
AXIAL CHARCOAL BED EFFICIENCY 
RADIAL CHARCOAL BED EFFICIENCY 
AXIAL LiOH BED EFFICIENCY 
CO OXIDIZER EFFICIENCY 
CATALYTIC OXIDIZER EFFICIENCY 
CONDENSATE EFFICIENCY 
CONTAMINANT MATERIAL BALANCE 
SUM OF MASS REMOVED BY DEVICES 
LOAD GENERATION INTO DD COL 19 
CALCULATE FINAL AND AVERAGE CONC 
CALCULATE LiOH USED IN INCREMENT 
PRINT DATA TO THE SPECIFIED DEVICE(S) 
PRINT CONCENTRATION DATA ANSWERS 
PRINT DATA HEADINGS 
PRINT TOXIC HAZARD INDEX ANSWERS 
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Subroutine PCSET 

The subroutine PCSET is the precalculation setup routine. PCSET gets calculations started by 
assuming an initial cabin concentration before the program enters the time calculation loop. Figure 2 
shows a flow diagram of PCSET. PCSET sets the initial time increment ending time to V 240 of the basic 
time increment specified in the device definition input file. Subroutine PRAFIL is then called and 
columns 17 to 21 are zeroed. These columns are used to store the results of subsequent calculations. 
CNRSUB is called to calculate each device removal efficiency for an assumed initial contaminant con- 
centration of 1x10 ^ 20 mg/m 3 . Contaminant removal rates and predicted, equilibrium, and final cabin 
concentrations are calculated by subroutine MASBAL. These calculated values are copied from matrix 
DD to the calculation matrix, CC, and printed out by subroutines CROUT and RROUT if required. 



Figure 2. Subroutine PCSET block flow diagram. 
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The subroutine CNRSUB calculates the removal efficiency of each device for each contaminant 
in the simulation during every time increment. This calculation is based on the average calculated cabin 
concentration. Figure 3 shows a block flow diagram for CNRSUB. These calculations are conducted by 
device type rather than the relative positions of each device with respect to each other. 
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This routine sets the cabin removal efficiency to zero and the leakage efficiency to the maximum 
of 1.0. Removal efficiencies for all devices with no flow are also set to zero. The remaining device effi- 
ciencies are calculated by calling the subroutines ACHBD, RCHBC, ALIOH, COOXID, CATBNR, and 
CONDHX. These calculated efficiencies are stored in matrix DD. 

Subroutine MASBAL 

Using the device efficiencies calculated by CNRSUB, MASBAL determines the mass removed, 
the calculated cabin concentration, the equilibrium cabin concentration, and final cabin concentration for 
each contaminant during a time increment. This calculation is conducted for all removal devices in 
parallel and in series. Figure 4 shows a block flow diagram for MASBAL. 

MASBAL uses the mass of contaminant removed and the net mass to the cabin to determine the 
final cabin concentration for each contaminant. The mass of the contaminant removed is defined as the 
product of the removal device flow rate, contaminant concentration, and device removal efficiency. The 
net mass of contaminant to the cabin is defined as the difference between the mass generated and mass 
removed. At steady state or equilibrium, the mass removed equals the mass generated. The mass gener- 
ated is the sum of all generation sources which includes the cabin generation rate and the generation rate 
in each device. The steady-state concentration is defined according to the following equation: 


Css ~ (Wnet-to-cabinV(t]r^'l3) * 

where m ne t-to-cabin is the mass of contaminant, x] r is the overall removal efficiency for all devices, and Q 
is the atmospheric flow rate through the removal devices. 

MASBAL is composed of two parts to determine the steady-state concentration. The first part of 
MASBAL determines the product of the overall efficiency and flow rate by setting the device generation 
rates to zero, assuming an arbitrary value for average cabin concentration (100 mg/m 3 ) and cabin gen- 
eration rate (50 mg/h), and calling CALCM to determine the sum of mass removed for all the removal 
devices. The second part of MASBAL evaluates the net mass to the cabin by setting the average cabin 
concentration equal to zero, restoring the contaminant and device generation rates to the values specified 
in the contaminant data matrix, and calculating the mass removed using CALCM. The m ne t-to-cabin equals 
the difference between the masses generated in the cabin and removal devices and the mass removed. 
From these values, C„ is calculated according to equation (8). After calculating the steady-state concen- 
tration, the final and average cabin concentrations are calculated by calling PCAVCF, and CALCM is 
called to calculate the mass removed by the cabin and each device using the average calculated cabin 
concentration. 

Subroutine CALCM 

The removal device inlet and outlet concentrations and the total mass removed by the cabin and 
the specified removal devices is calculated CALCM by using the removal efficiencies, generation rates, 
and average cabin concentration. This calculation is sequential from one device to another and uses the 
outlet concentration of an upstream device as the inlet concentration for a downstream device. This cal- 
culation requires the device definition input data to be arranged to allow calculations for all upstream 
devices to be completed before calculations for the downstream devices. Figure 5 shows a block flow 
diagram for CALCM. 

The subroutine sets the cabin and leakage device inlet and outlet concentrations equal to the 
average cabin concentration. All other devices are tested for zero flow. Devices with zero flow have 
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Figure 4. Subroutine MASBAL block flow diagram. 
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(^BEGIN^) 

I 


UPSTREAM 


INPUT C av cabin 


LOAD DEV 1 CABIN 
Cjn = C av ca bjn DD(1,17) 
^out = < ^av cabin DD(1 ,18) 
M rem = 0DD(1,21) 

I 


i 


DEVICE REF 1=0 

^ 


,, N 

o o 
II II 

d d 

Ql “ Qrefl 
DD(DD(J.4),2) 
Cl ” CreM 
DD(DD(J,4),18) 




l: 




LOOP FOR DEVICE 2-15 
CALCULATE SUM M rem 
DD(J,21) 

* 

LOOP FOR DEVICE 1-15 
CALCULATE SUM M r 
DD(J,19) 


r| gen 


^ av in = 

((Ci Q-) + C2Q2 + C3Q3)/ 

(Ql + Q2+Q3)) 

+ Mg en dev/Qdev 


C in =0DD(J,17) 
C()ut = 0 DD(J,18) 
M rem =0DD(J,21) 


PRINT WARNING 


SUM M rem cabin - SUM Mg en 
cabin”SUM M rem dev 2-15 

DD(1 ,21) 


RETURN 




Figure 5. Subroutine CALCM block flow diagram. 















their inlet concentration, outlet concentration, and mass removed set equal to zero. Upstream devices for 
each removal device are identified. If the upstream device type is 1 or 2, the inlet concentration is set 
equal to the average cabin concentration plus any internal device generation rate divided by the device 
flow rate. Upstream device types other than 1 or 2 cause the device inlet concentration to be based on the 
flow rates and outlet concentrations of all the upstream devices. 

Outlet concentration and the mass removed by the devices are calculated according to the follow- 
ing equations: 


Cout = Cind-TJr) , (9) 

m rem = C in (e)(7)r) . ( 1Q ) 

The inlet concentration for a device with an upstream device is set equal to the outlet concentra- 
tion for the upstream device. A device with multiple upstream devices requires the mixing of streams 
with varying concentrations to be considered. For example, the inlet concentration for a device with 
three upstream devices must be calculated according to the following equation: 

C 4 = {C\Q\+C2Q.2+Ct,Qi)IQ.4 • 

The sum of the mass removed and mass generated is calculated by adding the masses removed 
and masses generated by all the devices. The difference between the sum of the mass generated and the 
sum of the mass removed gives the mass remdved by the cabin. 

Subroutine LPIGEN 

Subroutine LDIGEN is called by MASBAL to load the generation rates from matrix CDI column 
1 and columns 10 through 22 into matrix DD column 19. 

Subroutine PCAVCF 

Subroutine PCAVCF is called by MASBAL to calculate the increment final and average cabin 
concentrations for each contaminant. 

Subroutine RINCDD 

Subroutine RINCDD is used at the beginning of each time increment to input and operate on the 
time-dependent data. A flow diagram of RINCDD is shown by figure 6. This subroutine checks the 
time-dependent data to determine whether any changes occur during the current time increment. Vari- 
ables in matrix TT are identified. If a contaminant generation rate is indicated, the new rate is placed in 
the calculation matrix, CDI. Likewise, if a change in removal device flow rate or any other device 
change is indicated, the new information is placed in the appropriate device definition matrix, DD, loca- 
tion. 

Subroutine REGEN 

Figure 7 shows a block flow diagram for subroutine REGEN. This subroutine determines 
whether any charcoal or LiOH beds will be regenerated during the current time increment If regenera- 
tion occurs, the mass of contaminants stored in the beds is set equal to zero. Similarly, if the regenera- 
tion duration lasts for the entire time increment, the device flow rate is set equal to zero. 
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Figure 6. Subroutine RINCDD block now diagram. 









Figure 7. Subroutine REGEN block flow diagram. 











The first check conducted by the subroutine is for device type. Only charcoal and LiOH beds 
may be regenerated. Only regeneration cases which begin at the time increment beginning or 
regeneration cases which last for one or more complete time increments are treated. Data concerning the 
regeneration interval, duration, and first regeneration time are obtained trom matrix DD. For the LiOH 
bed regeneration, the duration is set equal to zero since bed changeout is assumed to occur quickly. The 
initial time, regeneration time, and regeneration duration are then checked to determine whether they are 
exact multiples of the basic time increment. If they are not, they are rounded to the next lowest multiple 
of the time increment and a warning is written to the screen. 

The next checks conducted by the routine determine whether regeneration occurs at the 
beginning of a time increment and whether the regeneration lasts for the entire increment. Regeneration 
for the entire increment causes the program to deactivate this device for that increment by setting the 
device flow rate equal to zero. Regeneration at the beginning ot an increment causes the sum of the mass 
removed by that device to be set equal to zero. For an LiOH bed, the total mass of LiOH used is also set 
equal to zero. 

Subroutine MCALC 

Calculation of the removal efficiency, mass removed, and calculated, equilibrium, and final cabin 
concentrations for each contaminant and each removal device is controlled by MCALC. These calcula- 
tions are based on the cumulative mass removed for each contaminant during the previous time incre- 
ment. Figure 8 shows a block flow diagram for MCALC. 

MCALC calls the subroutine PREDCT to calculate the average predicted cabin concentration 
based on the removal efficiency and the cumulative mass of contaminant removed during the previous 
time increment and the generation rate during the present time increment. The average predicted concen- 
tration is used by the subroutine CONVRG to calculate a new removal efficiency, mass removed, and 
average calculated, equilibrium, and final cabin concentrations. The predicted and calculated concentra- 
tions are compared in CONVRG and recalculated until the difference between them is less than the con- 
vergence error specified in the device definition data file, matrix DD. This recalculation and comparison 
continues for 20 iterations with the full time increment or until the difference is less than the conver- 
gence error. 

If the convergence error is still exceeded after 20 iterations, another loop with a maximum of 20 
iterations is entered which uses one-twentieth the basic time increment for the calculation. This loop 
ends as soon as the difference between the predicted and calculated values is less than the convergence 
error or 20 iterations have been completed. If convergence is not attained after this loop, the program 
writes a warning to the screen indicating that the calculation for the contaminant did not converge. The 
loop using the one-twentieth time increment is used only for a contaminant that does not converge dur- 
ing the first 20 iterations. This is more efficient than reducing the time increment for all the contaminant 
calculations. 

Subroutine PREDCT 

MCALC calls the subroutine PREDCT to calculate the average predicted cabin concentration for 
each contaminant during each time increment. This calculation is based on the removal efficiency and 
sum of contaminant mass removed in the previous increment and the generation rate during the present 
increment. Figure 9 shows a block flow diagram for PREDCT. 
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Figure 9. Subroutine PREDCT block flow diagram. 

PREDCT calls the subroutine PRAFIL to zero the part of matrix DD required for storing the 
calculation results. Data from the previous time increment are obtained by LODEFF, and MASBAL is 
called to calculate the efficiency and concentration. 

Subroutine LODEFF 

Subroutine LODEFF loads the efficiency calculated in the preceding increment from matrix CC 
to matrix DD. 

Subroutine CONVRG 

CONVRG is the main convergence loop subroutine. This subroutine calculates the average cabin 
concentration and compares it with the predicted cabin concentration for each contaminant during every 
time increment. Figure 10 shows a block flow diagram for CONVRG. 

CNRSUB is called by CONVRG to calculate the removal efficiency for each device using the 
predicted cabin concentration. Based on this removal efficiency, MASBAL calculates the average, final, 
and equilibrium cabin concentrations and the mass of contaminant removed by each device. The pre- 
dicted and calculated cabin concentrations are compared by calculating the absolute value of the differ- 
ence of the predicted and calculated concentration divided by the predicted concentration and comparing 
that value to the convergence error. If the absolute value of the comparison is less than the convergence 
error, convergence has been achieved and the iteration stops for that contaminant. If convergence has not 
been achieved, a new cabin concentration is calculated using a bisection technique after the first iteration 
and a Newton-Raphson technique for each additional increment. The loop counter value passes back to 
MCALC which determines whether convergence has been reached within 20 iterations. 
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Subroutine SLIOH 


SLIOH is the subroutine which calculates the cumulative mass of LiOH used during the simula- 
tion run. This calculation is cumulative since the mass of LiOH consumed during the present increment 
is added to the mass consumed in all the previous increments. 


Contaminant Removal Device Calculation Subroutines 

The contaminant removal device subroutines are supported by a substantial amount of theoretical 
and experimental data. A brief description of each subroutine is provided in addition to a discussion of 
the supporting theory and experimental data. 

Subroutine ACHBD 

Subroutine ACHBD calculates the removal device efficiency for an axial flow charcoal bed. This 
routine simulates the physical adsorption of contaminants onto the surface of the charcoal. Specially 
treated charcoals are also considered which include chemical reaction between the surface treatment and 
the contaminant in addition to adsorption. 

A charcoal bed is composed of two zones during the adsorption process. These zones are desig- 
nated as the saturated zone and the adsorption zone. All contaminant removal takes place in the adsorp- 
tion zone. The saturated zone provides no net removal since it is in equilibrium with the vapor phase 
contaminant composition. Figure 1 1 illustrates the zones simulated by the program graphically. Physical 
adsorption is an equilibrium process which depends on variables such as the contaminant vapor pressure, 
inlet concentration, molar volume, and cabin temperature. Studies conducted by Robell investigated the 
thermodynamics of adsorption dynamics and developed a correlation between the physical properties of 
a contaminant and the charcoal saturation capacity. This correlation is based on the Polanyi Potential 
Theory and the Gibbs equation. From this study, a correlation factor, called the adsorption potential 
factor was developed. This factor is defined according to the following equation: 


A — (T/V m )log w (p v /p c ) » 


( 12 ) 


where T is the cabin temperature in Kelvin, V m is the contaminant liquid molar volume in cm 3 /gram- 
mole, p v is the contaminant vapor pressure at the cabin temperature expressed in concentration units of 
mg/m 3 , and p c is the cabin contaminant partial pressure expressed in concentration units of mg/m 3 . This 
factor was plotted as a function of experimentally determined charcoal saturation capacities to obtain the 
plot shown by figure 12. 4 The plot in this Figure was constructed for Barnebey-Sutcliffe type BD granu- 
lar activated charcoal. This correlation is not only sensitive to charcoal impregnation and contaminant 
solubility, but also to relative humidity as shown by figure 13. 5 Additional information may be obtained 
on charcoal capacity and performance from references 6 and 7. Based on potential plots, empirical 
equations are obtained which relate the potential factor to the charcoal saturation capacity. The equations 
used in this program are functions of the adsorption potential factor, contaminant solubility, and cabin 
relative humidity. Specific equations used are found in subroutine FQI in the ACHBD. FOR listing found 
in appendix A. The general form of the equations is the following: 

q = ae-P A , (13) 

where q is the charcoal saturation capacity in cm 3 of liquid contaminant per gram of charcoal and A is 
the adsorption potential factor in Kelvin-gmol/cm 3 . As new information concerning adsorption capacity 
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LEGEND 

tj = Bed profile at the time when a bed segment reaches steady state 

tb = Bed profile at the service time when the bed outlet concentratio equals Cb 

I = Active adsorption zone length 

Cj = Bed inlet concentration 

Cb = Bed penetration concentration 

Figure 11. Charcoal saturation and adsorption zone distribution. 

is obtained, these equations can be modified accordingly. In addition, this technique can be applied to 
other adsorbent materials to simulate other packing materials besides charcoal. 

The adsorption zone length for ninety percent removal is determined from experimental data 
obtained by Olcott at a 0.0066 m/s (1.3 ft/min) flow rate. 8 This data is plotted in figure 14 and the com- 
puter program uses the following equation to calculate the adsorption zone length: 

fads = aadsa l 1.3ft/m 1 n)(V/1.3)0- 8 . (14) 

The adsorption zone length study conducted by Olcott shows that the adsorption zone length 
increases with velocity to the 0.8 power as indicated in the equation. The saturated zone length is based 
on the charcoal capacity at the prevailing cabin conditions and the amount of contaminant already 
adsorbed. For a given contaminant mass retained in the bed, the saturation zone length equals the mass of 
contaminant adsorbed divided by the saturation capacity, q. The total bed length minus the saturated zone 
length equals the adsorption zone length. The adsorption zone length is the length of the bed actually 
available for contaminant removal. The program calculates the saturation zone length as calculated based 
on the bed geometry, the amount of contaminant removed by the bed, and the saturation capacity, q. 
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Figure 14. Adsorption zone length as a function of the potential factor. 

Adsorption of multiple contaminants by charcoal involves some interaction between the contam- 
inants. This interaction, called blockage or coexistence, reduces the capacity of the charcoal to hold 
other contaminants. Although the theory for coexistence is complex, experimental data indicates that an 
additional twenty percent can be added to the saturation zone when the calculation for the adsorption 
zone is conducted. 

Reaction With Specially Treated Charcoals . Some contaminants are not readily removed by 
granular activated charcoal but can be removed by charcoal which has been specially treated with 
chemicals that react with the contaminant after adsorption onto the charcoal surface. Two commonly 
used treated charcoals target ammonia and formaldehyde. 

Ammonia Removal . Ammonia is removed by treating granular activated charcoal with phos- 
phoric acid. Usually, phosphoric acid loading is 1.22 mmol/g of charcoal. This results in a requirement 
of 0.0061 g of charcoal per gram of ammonia to be removed if the reaction goes to completion. The bed 
removal efficiency is typically 90 to 99 percent per pass for a fresh bed. As the phosphoric acid is 
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depleted, the efficiency drops, eventually reaching zero. This routine assumes that the removal 
efficiency is 100 percent if the bed is less than 80 percent utilized. The efficiency for last 20 percent of 
the bed is calculated using the following sine relationship: 

Or = Sin (^charcoal - ^treated charcoal usedV(0.2)/W c h a rcoal > (15) 

where mcharcoai is the mass of charcoal in the bed and w trea ted charcoaj used is the mass of treated charcoal 
used. 


Formaldehyde Removal . Formaldehyde is removed most efficiently by chromate impregnated 
charcoal. Manufacturer’s data indicates that this charcoal can chemisorb a total amount of formaldehyde 
equivalent to 5 percent of its weight. 9 Testing at Lockheed Missiles and Space Company, Inc., showed 
this material’s efficiency to drop linearly from 100 to 90 percent for an amount of formaldehyde 
chemisorbed from 0 to 0.12 percent of the bed weight. Also, if the bed residence time is less than 0.25 s, 
the removal efficiency drops linearly. 

Subroutine RCHBP 

The subroutine RCHBD uses the same logic as ACHBD for simulating charcoal adsorption. 
However, this routine accommodates the geometry of a radial flow charcoal bed. 

Subroutine ALIOH 

ALIOH simulates removal of acidic contaminants by granular lithium hydroxide and lithium 
carbonate. The amount of lithium hydroxide consumed per weight of contaminant is calculated from the 
reaction stoichiometry. This number is input with the contaminant data. Reaction of lithium hydroxide 
with carbon dioxide to produce lithium carbonate has no noticeable effect on the contaminant removal. 
The removal efficiency for this device is 100 percent unless the bed is less than 1.905-cm thick or more 
than 80 percent utilized. The drop in efficiency as the bed is utilized is approximated by the following 
sine relationship: 


J] r = sin (WLiOtr m LiOH usedV(0.2)mLiOH . (16) 

where mu oh is the mass of the lithium hydroxide bed and muon used is the mass of lithium hydroxide 
utilized. Efficiency for a bed less than 1.905-cm thick decreases linearly with thickness. 

Subroutine COOXID 

Ambient temperature catalytic oxidation of carbon monoxide and hydrogen is simulated by the 
COOXID subroutine. This routine simulates ambient temperature catalytic oxidation using a granular 
activated charcoal with 2 weight percent platinum loading. This simulation is effective only for carbon 
monoxide and hydrogen. Efficiency remains constant at 100 percent per pass unless the residence time 
falls below 0.2 s. For residence times below 0.2 s, the efficiency decreases linearly according to figure 
15. 10 


Subroutine CATBNR 

The CATBNR subroutine simulates the destruction of hydrogen, carbon monoxide, methane, and 
other low molecular weight organic contaminants into carbon dioxide and water vapor using high tem- 
perature catalytic oxidation. The degree of oxidation in the oxidizer must be input by the user in the 
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contaminant data input file. Typically, oxidation efficiency is based on experimental oxidation 
performance testing. On average, operating the oxidizer at 400 °C (750 °F) provides removal efficiency 
of 100 percent for most contaminants. 

Subroutine CONDHX 

The subroutine CONDHX simulates the removal of contaminants by absorption into humidity 
condensate in a condensing heat exchanger. Some contaminants are removed by this route not only by 
absorption but also by chemical reaction in the condensate. Ammonia is treated in this manner since it 
dissociates in water and reacts with dissolved carbon dioxide. 11 All other contaminant removal is simu- 
lated using Henry’s Law. Using Henry’s Law is justified for trace contaminants since their concentra- 
tions in the atmosphere approach infinite dilution. Henry’s Law correlates the concentration of a 
contaminant in the atmosphere to its concentration in the liquid phase. The correlation coefficient is the 
Henry’s Law Constant, H, which has units of atmospheres per mole fraction. Equation (17) shows the 
Henry’s Law relationship in which p c is the contaminant partial pressure in atmospheres, H is the 
Henry’s Law constant in atmospheres per mole fraction, and x is the liquid phase mole fraction. 

p c = Hx. (17) 

The simulation assumes that the absorption process is concurrent and that equilibrium is very 
closely approached. A material balance on this process provides a relationship for the condensate mole 
fraction shown by: 


x = yl[(C/A)+(H/P )] . (18) 

In this equation, x is the liquid phase mole fraction, y is the vapor phase mole fraction, C is the conden- 
sate mass molar flow rate in mol/h, A is the atmospheric molar flow rate in mol/h, H is the Henry’s Law 
constant in atmospheres, and P is the total pressure in atmospheres. Figure 16 illustrates the absorption 
process. Based on the cabin concentration, the program calculates the inlet mole fraction based on a 1 
atmosphere total pressure. The condensate flow rate and atmosphere flow rate are entered in the device 
definition data and converted to molar flow rates based on 1 atmosphere pressure and 294 K absolute 
temperature. The mole fraction of contaminant leaving in the condensate is used to determine the mass 
of contaminant removed. The removal efficiency is calculated from the ratio of the difference in mass of 
contaminant entering and mass of contaminant removed to the mass of contaminant entering. 

Ammonia removal is treated separately since it reacts chemically with dissolved carbon dioxide 
in the humidity condensate. According to reference 22, data correlating ammonia partial pressure to 
liquid phase ammonia concentration for several carbon dioxide atmospheric partial pressures was used 
to obtain an equation relating liquid and gas phase ammonia composition. This data was obtained by 
sparging a gas mixture through a volume of water. Figure 17 shows a plot of the result. 12 The carbon 
dioxide curve corresponding to 666.6 Pa (5 mm Hg) was used to obtain equation (19) which relates 
ammonia mass per kilogram of condensate to the entering ammonia concentration. 

m a = 189.6C c 0 - 535 . (19) 

In this equation, m a is the mass of ammonia in milligrams per kilogram of condensate and C c is heat 
exchanger inlet ammonia concentration in mg/m 3 . This equation is used to determine the ammonia 
removal efficiency from the mass of ammonia entering and leaving the condensing heat exchanger 
assembly. 
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Data Input and Output Subroutines 


The following subroutines regulate the data input and output for each computer simulation run. 
These subroutines have been designed to allow maximum flexibility for calculated data output to facili- 
tate data analysis reporting. 

Subroutine CRIN 

Subroutine CRIN is called by MAIN and reads contaminant input data into the contaminant 
name matrix, NN, and the main calculation matrix, CDI. 

Subroutine RRIN 

Subroutine RRIN is called by MAIN and reads device definition data and time-dependent data 
into matrices DD and TT, respectively. 

Subroutine CROUT2 

CROUT2 is called by MAIN and controls output of the contaminant input data to the printer or 
computer terminal screen. One row at a time without headings is written to these output devices for the 
user to review before entering the calculation loop. 

Subroutine RROUT2 

Subroutine RROUT2 is called by MAIN and controls output of the device definition data and 
time-dependent data to the printer or computer terminal screen. One row at a time without headings is 
written to these output devices for the user to review before entering the calculation loop. 

Su b rou ting CRO U T 

CROUT is called by MAIN and regulates output of matrix CC data during each time increment 
for diagnostic purposes. This subroutine is called only when print switch No. 5 is set equal to 1 . 

Subroutine RROUT 

Subroutine RROUT is called by MAIN and regulates output of matrix DD for diagnostic purposes. 
This subroutine is called only when print switch No. 5 is set equal to 1. 

Subroutine PATOUT 

DATOUT is called by MAIN and serves as the master output regulation routine. Routines contained 
within PRFANS are called from DATOUT according to the print switch designations made by the user. 

Subroutine GROUP 

Subroutine GROUP is called by DATOUT and calculates the toxic hazard index according to 
appendix B. This subroutine also regulates the output for the toxic hazard index for both the standard 
formatted output and the plot data output. 
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Subroutine PRFANS 

PRFANS contains several subroutines that are called by DATOUT which regulate the output for 
contaminant concentration data, sum of contaminant masses removed data, and removal device effi- 
ciency data. This subroutine regulates output for both the standard formatted output and the plot data 
output. 

Subroutine HEAPGS 

Subroutines within HEADGS are called by PRFANS subroutines to regulate standard formatted 
data output headings. Headings are provided for contaminant concentration data, contaminant removal 
rate data, sum of contaminant masses removed data, and removal device efficiency data. 
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This appendix contains listings for each major subroutine and the main TCCS computer 
program. The main program is listed first followed by listings of each subroutine. The subroutines 
listings are arranged in alphabetical order by name to provide easy reference. 


RM/ FORTRAN Compiler (V2.42) 

Source File: C:\FORTRAN\TCC\MAIN.F Options: /C 80 /L /BIJY 
1 C FILE: MAIN. FOR 


03/15/94 13:23:35 


15 C 

16 C 


19 C 

20 C 

21 C 

22 C 

23 C 

24 

25 C 

26 C 


29 C 

30 C 

31 C 

32 C 


35 C 

36 C 

37 C 

38 C 

39 C 

40 C 

41 

42 

43 

44 C 

45 


PROGRAM TCCS 81 

***************************************************************** 

* ★★★★•A:*********************************************************** 

* PROGRAM TO MODEL REMOVAL OF SPACECRAFT * 

* GASEOUS CONTAMINANTS * 

* VERSION 8.1 Alpha * 

* March 15, 1994 * 

* ★★★★★I*******************************************************.*.*** 

a************:**************************************************** 

SUBROUTINES REQUIRED: 

CAFILL-FILL MATRIX WITH ZEROS 
RAFILL-FILL MATRIX WITH ZEROS 
CRIN-READ IN INPUT DATA 
RR IN -READ IN INPUT DATA 
CROUT- PRINT OUT INPUT DATA 
RROUT- PRINT OUT INPUT DATA 
PCSET- PRECALCULATION SET UP ROUTINE 
MCALC-MAIN CALCULATION ROUTINE 
DATOUT-DATA PRINTOUT ROUTINE 
XXXXX-TIME DEPENDENT DATA ROUTINE 
REGEN- REGENERATION OF DEVICES ROUTINE 
SLIOH-SUM LI OH USED IN TIME INCREMENT 

NOTE : SUBROUTINES USE ADJUSTABLE SIZE ARRAYS 
WATCH COMPILER OPTIONS /DIMENSIONING IF 
ANY ARRAY IS LARGER THAN 64K BYTES 

****************** DIMENSI0N majn PROGRAM MATRICES**************** 

NN=CONTAMINANT NAME MATRIX 

CDI=CONTAMINANT INPUT DATA MATRIX 

CC=CALCULATON MATRIX 

DD=DEVICE DEFINITION MATRIX 

TT=TIME DEPENDENT DATA MATRIX 

PRESENTLY SET TO HANDLE MAXIMUM OF 150 CONTAMINANTS IN MATRICES 
THIS VALUE =NROW AND IS USED IN ADDRESSING ADJUSTABLE SIZE 
ARRAYS IN SUBROUTINES 

NTTROW IS USED FOR MATRIX TT MAXIMUM LENGTH 

CHARACTER NN(300)*30 
REAL CDI (300,23) 

REAL CC ( 300 , 48 ) 

REAL DD (15,21) 

REAL DD (15,23) 

REAL TT (750, 7 ) 


51 C 

52 C 


NOTE: MUST COMPILE SUBROUTINES PROPERLY FOR ADJUSTABLE SIZE 
ARRAYS IF A MAIN MATRIX EXCEEDS 65536 BYTES (REALS =4 BYTES) 


************ 


DECLARE OTHER TERMS USED IN MAIN PROGRAM ********* 


P-AGfc 




‘ u * 1 ^ i y' \ 
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53 

54 

55 

56 

57 

58 C 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 C 

69 C 

70 C 

71 C 

72 C 

73 C 

74 C 

75 C 

76 C 

77 C 

78 C 

79 C 

80 

81 C 

82 

83 C 

84 

85 

86 C 

87 

88 

89 C 

90 

91 

92 C 

93 

94 C 

95 

96 C 

97 

98 C 

99 

100 C 

101 

102 C 

103 C 

104 C 

105 C 

106 

107 C 

108 

109 

110 
111 
112 C 


CHARACTER FNAME*24 , DES*1 , FCPLOT*24 , FTPLOT*24 , FEPLOT*24 
LOGICAL EX 

INTEGER PRTSWl , PRTSW2 , PRTSW3 , PRTSW4 , PRTSW5 , PRTSW6 , PRTSW7 , PRTSW8 , 
+ TVAL , IDEVNO, IDEVN1 , IDEVN2 , IDEVN3 , PRTSW9 


******* print WELCOME AND PROGRAM VERSION NUMBER 
WRITE ( * , 9 ) 


*************** 


FORMAT (IX 1 ************************************************** 1 / 


IX, 1 * 
IX, ' * 
IX, ' * 
IX, • * 
IX, ' * 
IX, ' * 


WELCOME TO THE WORLD 
OF THE 

SPACECRAFT ATMOSPHERIC TRACE CONTAMINATION 
CONTROL SIMULATION PROGRAM 
-VERSION 8 . 1 Alpha- 
March 15, 1994 


IX i ************************************************** • ) 


******* DEFINE PROGRAM VARIABLES ***********’ 
LIN=NO. OF LINES OF DATA IN MAT NN & MAT CDI 
LINl=NO . LINES OF DATA IN MAT TT 
LIN2 =NO . LINES OF DATA IN MAT DD 
TN= INCREMENT END TIME (HRS) 

TN1= INCREMENT BEGINNING TIME (HRS) 

TMIS=TOTAL MISSION TIME (HRS) 

NINC =NUMBER OF TIME INCREMENTS ELAPSED 


MAT NN , CC , AND CDI MUST HAVE SAME NO. OF ROWS 

DIMENSIONS OF MAT DD 

NROW=15 

NCOL=2 1 

NCOL=23 

DIMENSIONS OF MAT CC & ROWS IN MAT NN 

NROWl=300 

NC0L1=48 

DIMENSIONS OF MAT CDI 
NROW2 =NROWl 
NCOL2=23 

DIMENSIONS OF MAT TT 
NTTROW=7 50 
NTTCOL=7 

DEVICE NUMBER FOR OUTPUT DATA (SET TO 6 FOR FORM FEED ON OUTPUT) 
IDEVNO=6 

DEVICE NUMBER FOR MESSAGE OUTPUT 
IMSGDN=2 

DEVICE NUMBER FOR CONTAMINANT PLOT DATA 
IDEVNl=10 

DEVICE NUMBER FOR T-VALUE PLOT DATA 
IDEVN3=11 

DEVICE NUMBER FOR EFFICIENCY PLOT DATA 
IDEVN2=12 

*********** END OF DEFINITION SECTION ************************* 


*************** 2ER0 MATRICES* *********************************** 
PUT BLANKS IN NAME MATRIX 
Oil CALL CAFILL (NN, 1 ,NROWl ) 

PUT ZEROS IN OTHER MATRICES 
CALL RAFILL (CDI , NROW2 , NCOL2 ) 

CALL RAFILL (CC , NROW1 , NCOLl ) 

CALL RAFILL ( DD , NROW, NCOL ) 

CALL RAFILL (TT, NTTROW, NTTCOL) 
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113 C 

114 C 

115 

116 

117 

118 

119 

120 
121 

122 C 

123 

124 C 

125 

126 

127 

128 

129 

130 

131 C 

132 

133 

134 

135 

136 

137 

138 

139 C 

140 C 

141 C 

142 

143 

144 C 

145 C 

146 C 

147 C 

148 C 

149 C 

150 C 

151 C 

152 C 

153 C 
PLO 

154 C 

155 C 

156 

157 

158 

159 

160 
161 
162 

163 

164 C 

165 C 

166 

167 

168 

169 

170 

171 


******* READ in data FROM FILES AND PRINT IT IF DESIRED ******** 

010 WRITE (*,*) 'INPUT CONTAMINANT DATA FILE NAME: 

CALL CRIN (NN, CDI , NROW2 ,NCOL2 , LIN) 

012 WRITE (*,*)' PRINT CONTAMINANT INPUT DATA? (Y/N) 1 

READ ( * , ' (A) 1 ) DES 

IF( (DES.NE. ' Y' ) .AND. (DES.NE. 'N* ) ) GOTO 12 
IF {DES.EQ.'N') GOTO 20 

CALL CROUT2 (NN , CDI , NROW2 , NCOL2 , 1 , NCOL2 ,LIN, 1 , LIN, IMSGDN) 

020 WRITE (*,*) 'INPUT DEVICE DEFINITION TABLE FILE NAME: 

NOTE: ONLY 16 COLUMNS ARE IN THE INPUT FILE 
CALL RRIN (DD,NROW,NCOL, 16, LIN2 ) 

022 WRITE (*,*) 'PRINT DEVICE DEFINITION TABLE? (Y/N) ' 

READ ( * , ' (A) ' ) DES 

IF( (DES.NE. ' Y ' ) .AND. (DES.NE. ' N ' ) ) GOTO 22 
IF (DES . EQ. 'N' ) GOTO 30 

CALL RROUT2 ( DD , NROW, NCOL , 1 , 1 6 , LIN2 , IMSGDN) 

030 WRITE ( * , * ) 'INPUT TIME DEPENDENT DATA FILE NAME: 

CALL RRIN ( TT , NTTROW , NTTCOL , NTTCOL , LINl ) 

032 WRITE (*,*) 'PRINT TIME DEPENDENT DATA? (Y/N) ' 

READ ( * , ' (A) ' ) DES 

IF ( (DES.NE. ' Y' ) .AND. (DES . NE . ' N ' ) ) GOTO 32 
IF (DES.EQ.'N') GOTO 40 

CALL RROUT2 ( TT , NTTROW , NTTCOL , 1 , NTTCOL , LINl , IMSGDN ) 


************* READ IN MISSION DATA VARIABLES 
READ IN MISSION TOTAL TIME (HRS) 




040 WRITE (*,*) 'INPUT TOTAL MISSION TIME IN HOURS: 


READ ( * , * ) TMIS 


********** print SWITCH DEFINITION ***************************** 
1=RESULTS FOR ONE CONTAMINANT IN PCSET 

2=RESULTS FOR 1 CONT & INCR IN 1/10 INCR CONV ROUTINE (IN MCALC ) 

3 CONVERGENCE VALUES IN CONVRG 
4 =RESULTS FOR 1 CONT IN MCALC AFTER CAV CALC 
5=MAT CC AND MAT DD AT END OF TIME INCREMENT 

6= PRINT CONC+M . REM+SUM MASS REM+REM EFF (OTHERWISE ONLY CONC DATA) 

7= PRINT OUTPUT WITH NO FORM FEEDS 

8=PRINT ANSWERS DURING EACH ITERATION (IN MAIN PROGRAM) AND CONTROL 


9 CONTROL PLOT FILE OUTPUT 

TVAL = C ONTROL OUTPUT OF GROUP CONTRIBUTION T-VALUE DATA 
PRTSWl =NINT ( DD ( 2 , 9 ) ) 

PRTSW2 =NINT ( DD ( 2 , 1 0 ) ) 

PRTSW3 =NINT ( DD ( 2 , 1 1 ) ) 

PRTSW4 =NINT ( DD ( 2 , 12 ) ) 

PRTSW5=NINT ( DD (2,13) ) 

PRTSW6=NINT ( DD (2,14)) 

PRTSW7 =NINT ( DD ( 2 , 1 5 ) ) 

PRTSW8 =NINT ( DD ( 2 , 1 6 ) ) 


**** MAKE DECISION ON CONCENTRATION AND EFFICIENCY PLOT DATA **** 
WRITE (*,*) 'DO YOU WISH TO WRITE INCREMENT DATA TO A PLOT FILE?' 
WRITE (*,*) ' 1. Concentration Data (C)' 


WRITE ( 
WRITE ( 
WRITE ( 


* ) 
*) 
*) 


2. Efficiency Data (E) ' 

3. Both Concentration and Efficiency Data (B) ' 

4. Neither (N) ' 


WRITE 


ENTER YOUR SELECTION: 
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172 

173 

174 

175 

176 

177 

178 

179 

180 
181 

182 C 

183 C 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 C 

198 C 

199 

200 
201 
202 

203 C 

204 C 

205 C 

206 

207 C 

208 

209 

210 C 

211 

212 C 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 


READ (*, ' (A) ' ) DES 
IF ( DES . EQ . ' C ' ) THEN 
PRTSW9=1 

ELSEIF (DES.EQ.'E’) THEN 
PRTSW9=2 

ELSEIF (DES . EQ . ' B ' ) THEN 
PRTSW9 = 3 
ELSE 

PRTSW9=0 

ENDIF 

*************** make DECISION ON T-VALUE OUTPUT ***************** 
WRITE ( * , * ) ' PRINT GROUP CONTRIBUTION T-VALUE DATA? 1 

WRITE (*,*) 1 1. Print to Normal Output (Y) ' 

WRITE (*,*) ' 2. Print to Normal Output and Plot File (P) 1 

WRITE (*,*) ' 3. Do Not Print (N) ' 

WRITE (*,*) 'ENTER YOUR SELECTION: 

READ ( * , ' (A) ' ) DES 
IF (DES . EQ . 1 Y ' ) THEN 
TVAL=1 

ELSEIF (DES . EQ . 1 P ' ) THEN 
TVAL=2 
ELSE 
TVAL=3 
ENDIF 

CHANGE TO NO FORM FEED IF PRTSW7=1 
IF ( PRTSW7 . EQ . 1 ) THEN 
IDEVNO=7 
ENDIF 

*************** make DECISIONS ON DATA OUTPUT ****************** 
THIS IS WHERE ALL PROGRAM OUTPUT DATA FILES ARE OPENED 
THEY MUST BE CLOSED AT THE END OF THE PROGRAM 

SECTION WHICH CHECKS FOR EXISTANCE OF OUTPUT FILE & OPENS IT 
050 WRITE ( * , * ) 1 WRITE OUTPUT TO FILE, PRINTER, SCREEN, OR END?' 

WRITE (*,*)’ (FILE NAME/LPT1/CON/END) ' 

*********** NOTE: LPT1 OUTPUT REQUIRES 132 COLUMNS *********** 

READ ( * , ' (A) ’ ) FNAME 

QUIT IF FNAME=END 

IF ( FNAME . EQ . ' END 1 ) GOTO 999 

IF ( (FNAME. NE. 1 LPT1 ' ) .AND. (FNAME. NE. 'CON' ) ) THEN 
INQUIRE ( FILE=FNAME , EXIST=EX) 

IF (EX) THEN 

WRITE (*,*)' FILE EXISTS - OVERWRITE? (Y/N) ' 

READ ( * , ' (A) ' ) DES 
IF (DES . NE . 1 Y ' ) THEN 
GOTO 50 
ELSE 

OPEN ( IDEVNO, FILE=FNAME , STATUS= ' OLD ' , IOSTAT=IOVAL ) 

ENDIF 

ELSE 

OPEN (IDEVNO, FILE = FNAME, STATUS = 'NEW' , IOSTAT=IOVAL ) 

ENDIF 

ENDIF 

IF ( (FNAME. EQ. 'LPT1 ') .OR. (FNAME. EQ. 'CON' ) ) THEN 
OPEN ( IDEVNO, FILE=FNAME, IOSTAT=IOVAL) 

ENDIF 

IF ( I OVAL .NE . 0 ) THEN 
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CLOSE ( IDEVNO ) 

WRITE ( * , * ) ' IOERROR= ' , IOVAL 
GOTO 50 
END IF 

*** OPEN FILE FOR CONCENTRATION PLOT DATA IF PRTSW9=1 OR 3 * 

IF ( (PRTSW9.EQ.1) .OR. (PRTSW9.EQ.3) ) THEN 

WRITE (*,*) 'FILE NAME FOR CONCENTRATION PLOT DATA OUTPUT? 

READ ( * , ' (A) 1 ) FCPLOT 
INQUIRE ( FILE=FCPLOT , EXIST=EX) 

IF (EX) THEN 

WRITE (*,*) 'PLOT FILE EXISTS - OVERWRITE? (Y/N) ' 

READ ( * , ' (A) 1 ) DES 
IF (DES.NE. ' Y' ) THEN 
GOTO 52 
ELSE 

OPEN (UNIT=10 , FILE=FCPLOT, STATUS= ' OLD' , IOSTAT=IOVAL) 
END IF 
END IF 

OPEN (UNIT=10 , FILE=FCPLOT, STATUS= ' NEW' , IOSTAT= IOVAL) 

END IF 

IF ( IOVAL. NE.0) THEN 
CLOSE (UNIT=10 ) 

WRITE ( * , * ) ' IOERROR= 1 , IOVAL 

GOTO 52 
ENDIF 

******* OPEN FILE FOR EFFICIENCY PLOT DATA IF PRTSW9=2 OR 3 

IF ( (PRTSW9.EQ.2) .OR. (PRTSW9.EQ.3) ) THEN 

WRITE (*,*) 'FILE NAME FOR EFFICIENCY PLOT DATA OUTPUT? ' 
READ < * , 1 (A) * ) FEPLOT 
INQUIRE ( FILE=FEPLOT , EXIST=EX) 

IF (EX) THEN 

WRITE (*,*) 'PLOT FILE EXISTS - OVERWRITE? (Y/N) ' 

READ ( * , 1 (A) ' ) DES 
IF (DES.NE. 'Y') THEN 
GOTO 54 
ELSE 

OPEN (UNIT=12 , FILE=FEPLOT, STATUS= 1 OLD 1 , IOSTAT= IOVAL) 
ENDIF 
ENDIF 

OPEN (UNIT=12 , FILE=FEPLOT, STATUS= ' NEW' , IOSTAT=IOVAL ) 

ENDIF 

IF ( IOVAL. NE.0) THEN 
CLOSE (UNIT=12 ) 

WRITE ( * , * ) ' IOERROR= ' , IOVAL 

GOTO 54 
ENDIF 

********** OPEN FILE FOR T-VALUE PLOT DATA IF TVAL=2 ******* 
IF (TVAL.EQ.2) THEN 

WRITE ( * , * ) ’ FILE NAME FOR T-VALUE PLOT DATA OUTPUT? ' 

READ ( * , ' (A) ' ) FTPLOT 
INQUIRE (FILE=FTPLOT, EXIST=EX) 

IF (EX) THEN 

WRITE {*,*) 'PLOT FILE EXISTS - OVERWRITE? (Y/N) ' 
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READ (*, ' (A) * ) DES 
IF (DES.NE. ' Y' ) THEN 
GOTO 56 
ELSE 

OPEN (UNIT=11 , FILE=FTPLOT, STATUS^ ' OLD 1 , IOSTAT=IOVAL) 
ENDIF 
ENDIF 

OPEN (UNIT=1 1 , FILE=FTPLOT, STATUS^ 1 NEW' , IOSTAT=IOVAL) 
ENDIF 

IF (IOVAL.NE.0) THEN 
CLOSE (UNIT=11 ) 

WRITE ( * , * ) ' IOERROR= ' , IOVAL 

GOTO 56 
ENDIF 


CALL SYSTEM TIME AND DATE 

THIS MUST BE CALLED ONLY ONCE SO THAT THE TIME AND DATE WILL 
BE THE SAME ON ALL OUTPUT INFORMATION FOR ONE RUN 
CALL DATTM ( I MONTH , IDAY , I YEAR, IHOUR, IMINUTE, ISECOND) 

SET IPGCTR=COUNTER FOR SEQUENTIAL PAGE NUMBERS ON ALL OUTPUTS 
IPGCTR=0 


★ ★★★★★★★★★ 


********* check basic time increment ************************** 

********* BASIC TIME INCREMENT (HRS) -DD(1, 11) ***************** 

BINC=DD { 1 , 11) 

********* test for BINC= o (CAUSES ENDLESS TIME LOOP) 

IF (BINC.EQ.0) THEN 
CLOSE ( IDEVNO) 

WRITE { * , * ) ' ERROR-BASIC TIME INCREMENT=0 ' 

GOTO 999 
ENDIF 


************** ZERO INITIAL VARIABLES ************************* 

TN=0 

TN1 = 0 

NINC=0 


★ **★★**★★***★★★★★★*★★★★★*★**★★*★★*★★**★★★★★*★★**★★*★★★★★★★**★*★★★ 
* PRECALCULATION SET UP ROUTINE * 

FOR ALL CONTAMINANTS ONE AT A TIME AT CAV PRED=lE-20, CALC INIT 
DEV EFF AND LOAD IT INTO MAT CC-ALSO CALC CAVPRD(CAV PREDICTED) 
OUTPUTS TO PRECALC SET UP ROUTINE: 

TNl= INCREMENT INITIAL TIME (HRS) 

BINC=BASIC INCREMENT SIZE (HRS) (REF. =DD (1 , 11) ) PASS IN???? 

LIN=NO. OF CONT IN MAT CC AND NN 
DD, NROW, NCOL=NAME & SIZE OF MAT DD 
CC , NROWl , NCOLl =NAME & SIZE OF MAT CC 
CD I , NROW2 , NCOL2 =NAME & SIZE OF MAT CDI 
LIN2 =NO. DEVICES IN MAT DD 
NN=NAME OF MAT NN 

INPUTS FROM PRECALC SETUP ROUTINE-SUBROUTINE PCSET: 

PUT IN MAT CC 

CAVPRD=PRED CABIN AV CONC (MG/CU M) : =CC(I,2) 

CEQLIB=EQUILIBRIUM CABIN CONT CONC (MG/CU M):=CC(I,3) 

CFINAL=FINAL CABIN CONT CONC (MG/CU M):=CC(I,4) 

PUTS REM EFF FROM DD COL 20 IN CC (I, 7-10-13-16 ETC) 

PUTS M.REM IN CC ( I , 5- 8- 1 1- 14 . . ) 
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CALL PCSET ( TNI , LIN , DD , NROW , NCOL , CC , NROWl , NCOLl , 

+CDI , NROW2 , NCOL2 , LIN2 , NN , PRTSW1 , IMSGDN) 

************* 0F PRECALCULATION SETUP ROUTINE ******** 

********************************************************** 

* BEGINNING OF CALCULATION FOR EACH TIME INCREMENT 

********************************************************** 

100 CONTINUE 

********* INCREASE INCREMENT COUNTER ******************* 
NINC=NINC+1 

********* SET up TIME INCREMENT SIZE FOR INCREMENT **** 

IF (NINC.EQ.l) TN=BINC/24 
IF (NINC.EQ.2) TN=BINC/2 
IF (NINC.EQ.3) TN=BINC 
IF (NINC.GT.3) TN=TN+BINC 

CHECK FOR INCREMENT FINAL TIME > MISSION TIME 
IF (TN.GT.TMIS) TN=TMIS 

********* CHECK FOR CHANGES IN BASIC TIME INCREMENT **** 

BINCNEW=BINC 

DO 105 K = 1 , LINl 

IF ( (TT(K, 1) .GE.TN1) .AND. (TT(K, 1) .LT.TN) ) THEN 
IF ( ( TT ( K , 4 ) . EQ . 1 ) . AND . ( TT ( K , 6 ) . EQ . 1 1 ) ) THEN 
BINCNEW=TT (K , 7 ) 

ENDIF 

ENDIF 

105 CONTINUE 

IF (BINC .NE. BINCNEW) THEN 
TN=TN-BINC+BINCNEW 
BINC=BINCNEW 
DD { 1 , 11) =BINC 
ENDIF 

STORE PREVIOUS INCREMENT CABIN VOLUME 
PREVCVOL=DD (1,9) 

********* r E ad TIME DEPENDENT DATA *********************** 
CALL TIME DEPENDENT DATA SUBROUT INE-RINCDD 
CALL RINCDD ( I , TN , TNI , DD , NROW , NCOL , LIN2 , 

+CC , NROWl , NCOLl , CDI , NROW2 , NCOL2 , L IN , TT , NTTROW , NTTCOL , LINl ) 

CHECK TO SEE IF CABIN VOLUME HAS CHANGED, AND IF SO 
UPDATE INITIAL CABIN CONCENTRATION FOR NEW VOLUME 
IF ( PREVCVOL . NE . DD (1,9) ) THEN 
DO 200 1=1, LIN 

CC (1,1) =CC (1,6) /DD (1,9) 

200 CONTINUE 
ENDIF 

****** LIST INCREMENT NO. AND TIMES TO CONSOLE *********** 
IF ( (FNAME.NE. 'CON' ) .OR. ( PRTSW8 .NE. 1 ) ) THEN 
OPEN (IMSGDN, FILE= 1 CON 1 , IOSTAT=IOVAL) 

WRITE ( IMSGDN, 6 5 ) NINC , TNI , TN 

065 FORMAT (1X,'INCR NO.= ',15,' BEGIN & END TIMES (hours) 

+ F9.3,2X,F9.3) 
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CLOSE ( IMSGDN) 

ENDIF 

* * * * STORE ORIGINAL Q DEVICE IN DD COL 7 (TAKEN FROM DD COL 2)**** 
DO 110 J=1,LIN2 
DD ( J , 7 ) =DD ( J , 2 ) 

110 CONTINUE 

*********** CHECK FOR REGENERATION IN TIME INCREMENT *********** 

CALL REGENERATION SUBROUTINE REGEN 

CALL REGEN (TN, TNI , DD , NROW, NCOL , CC , NROWl , NCOL1 , 

+CDI , NROW2 , NCOL2 , LIN , LIN2 , IMSGDN) 

********** CALL MAIN CALCULATION SUBROUTINE ******************** 
OUTPUTS TO MAIN CALC SUBROUTINE-MCALC : 

I=CONTAMINANT NO. 

TN, TNI = INCREMENT END & BEGINNING TIME (HRS) 

DD, NROW, NCOL =NAME & DIM OF MAT DD 
CC, NROWl , NCOL 1=NAME & DIM OF MAT CC 
CD I , NROW2 , NCOL 2 =NAME & DIM OF MAT CDI 
NN=NAME OF MAT NN 

LIN=NUMBER OF CONTAMINANTS IN MAT NN & CDI 
LIN2 =NO . DEVICES IN MAT DD 
INPUTS FROM MAIN CALC ROUTINE -MCALC : 

TO MAT CC 

PUTS CAVCLC , CEQLIV , &CFINAL IN CC (1,2-3 &4) 

PUTS REM EFF FROM DD COL20 IN CC (1,7-10-13 ETC) 

PUTS M.REM FOR EACH DEV FROM DD COL2-1 IN CC ( 1 , 6-9-12 ETC) 

PUTS SUM MASS REM FOR EACH DEV IN CC (I, 8-11-14 ETC) 

CALL MCALC (I, TN, TNI, DD, NROW, NCOL, 

+CC , NROWl , NCOLl , CDI , NROW2 , NCOL2 , NN, LIN, LIN2 , 

+PRTSW2 , PRTSW3 , PRTSW4 , IMSGDN) 

************ CALCULATE LIOH USED IN INCREMENT ******************* 

CALL LIOH REMOVAL SUBROUTINE SLIOH 

CALL SLIOH (TN,TN1, DD, NROW, NCOL, CC, NROWl, NCOLl , 

+CDI , NROW2 , NCOL 2 , LIN, LIN2 ) 

************ restore device flow ******************************** 

RESTORE ORIGINAL DEVICE FLOW RATE FROM DD COL 7 TO DD COL 2 
DO 120 J=1 , LIN2 
DD ( J , 2 ) = DD ( J , 7 ) 

120 CONTINUE 

********* PRINTOUT OF DATA FOR EACH TIME INCREMENT ************** 
IF PRTSW5=1 THEN PRINT MAT DD+MAT CC INFO FOR THIS CONTAMINANT 
IF { PRTSW5 . EQ. 1 ) THEN 

OPEN ( IMSGDN, FILE=' CON' , IOSTAT=IOVAL) 

WRITE ( IMSGDN, *)' PRINTOUT FOR MAT CC & DD AT END OF TIME INCR' 
WRITE ( IMSGDN , * ) ' INFO FROM MAT CC ' 

CLOSE (IMSGDN) 

CALL CROUT (NN,CC, NROWl, NCOLl, 1, NCOLl, LIN, 1, LIN, IMSGDN) 

CALL CROUT (NN,CC, NROWl, NCOLl, 1 , NCOLl , LIN, 1 , LIN, IMSGDN, NINC, 

+ FNAME , IDEVNO , I OVAL ) 

OPEN ( IMSGDN, FILE= 'CON' , IOSTAT=IOVAL ) 

WRITE (IMSGDN, *)' INFO FROM MAT DD ' 

CLOSE (IMSGDN) 

CALL RROUT ( DD , NROW , NCOL , 1 , NCOL , LIN2 , IMSGDN ) 

CALL RROUT ( DD , NROW , NCOL , l,NCOL,LIN2, IMSGDN , FNAME , IDEVNO, I OVAL) 
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END IF 
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*********** REGULAR PRINTOUT OF DATA FOR EACH INCREMENT ******** 
IF ( PRTSW8 . EQ . 1 ) THEN 

CALL DATOUT ( TN , TNI , L IN , DD , NROW , NCOL , CC , NROWl , NCOL 1 , 

+ CDI , NROW 2 , NCOL 2 , LIN2 , NN, PRTSW6 , PRTSW8 , PRTSW9 , 

+ IDEVNO , NINC , IMONTH, IDAY, I YEAR, IHOUR, IMINUTE, FNAME , IOVAL, IPGCTR, 
+ TVAL, FCPLOT, IDEVN1 , IDEVN3 , IDEVN2 ) 

IF ( IOVAL. NE.0) THEN 
CLOSE (IDEVNO) 

WRITE (*,*)' PROGRAM DATA OUTPUT ERROR IN INCREMENT = 1 , NINC 
GOTO 999 
END IF 
ENDIF 

********* CONTROLS PLOT DATA OUTPUT IF PRTSW8=0 **************** 
IF (PRTSW8.EQ.0) THEN 

IF ( (PRTSW9.GT.0) .OR. (TVAL.EQ.2) ) THEN 

CALL DATOUT (TN , TNI , LIN , DD , NROW, NCOL , CC , NROWl , NCOL 1 , 

+ CDI , NROW2 , NCOL 2 , LIN2 , NN , PRTSW6 , PRTSW8 , PRTSW9 , 

+ IDEVNO, NINC, IMONTH, IDAY, IYEAR, IHOUR, IMINUTE , FNAME , IOVAL, 

+ IPGCTR, TVAL, FCPLOT, IDEVN1, IDEVN3 , IDEVN2 ) 

ENDIF 

ENDIF 


130 


********************* CHECK FOR END OF MISSION 




IF(TN.GE.TMIS) THEN 
END TIME LOOP 
CONTINUE 
ELSE 

********* UPDATE FOR NEXT TIME INCREMENT AND REPEAT ************ 
SET TFINAL FOR THIS INCR = TINIT FOR NEXT INCR 
TNl=TN 

SET CFINAL FOR INCR=CINIT FOR NEXT INCR-ALL CONTAMINANTS 
DO 130 1=1, LIN 
CC ( I , 1 ) =CC (1,4) 

CONTINUE 
GOTO 100 
ENDIF 


*+*★ + ★******★★ + ★★********** + *★**★***********★★★★**★*****★*★★*★*★★ 

* END OF CALCULATION FOR EACH TIME INTERVAL * 

*★★**★** + ****★**★★**★★*★ + ★★*★■*■******★*★★★★*******★**********★**** 


****** print FINAL ANSWERS AT END OF MISSION IF DESIRED ********* 
900 CONTINUE 

IF (PRTSW8 .EQ. 0) THEN 

CALL DATOUT ( TN , TNI , LIN , DD , NROW , NCOL , CC , NROWl , NCOL1 , 

+ CDI , NROW 2 , NCOL2 , LIN2 , NN, PRTSW6 , PRTSW8 , PRTSW9 , 

+ IDEVNO , -1, IMONTH, IDAY, IYEAR, IHOUR, IMINUTE , FNAME , IOVAL, IPGCTR, 
+ TVAL, FCPLOT, IDEVN1 , IDEVN3 , IDEVN2 ) 

IF ( IOVAL. NE.0) THEN 
CLOSE (IDEVNO) 

WRITE (*,*)' PROGRAM DATA OUTPUT ERROR - FINAL PRINTOUT' 

GOTO 999 
ENDIF 
ENDIF 

************ CLOSE ALL PROGRAM OUTPUT FILES ******************* 
CLOSE (IDEVNO) 
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CLOSE ( IDEVN1 ) 
CLOSE ( IDEVN3 ) 
CLOSE (IDEVN2) 


999 CONTINUE 

WRITE ( * , * ) 'DO YOU WISH TO RUN ANOTHER CASE? (Y/N) ' 

READ ( * , ' (A) ' ) DES 
IF (DES.EQ.'Y') THEN 
GO TO Oil 
ENDIF 

***************************************************************** 

* END OF MAIN PROGRAM * 

***************************************************************** 


END 


NUMBER 

OF 

WARNINGS 

IN 

NUMBER 

OF 

ERRORS 

IN 

NUMBER 

OF 

WARNINGS 

IN 

NUMBER 

OF 

ERRORS 

IN 


PROGRAM UNIT: 
PROGRAM UNIT: 
COMPILATION : 
COMPILATION : 


0 

0 

0 

0 
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RM/ FORTRAN Compiler (V2.42) 

source File: C:\FORTRAN\TCC\ACHBD. Options: /C 80 /L /BIJY 03/15/94 10:38:54 


1 C 

2 C 

3 C 

4 C 

5 C 

6 C 
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8 

9 C 

10 C 
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13 C 
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15 C 
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48 

49 

50 

51 

52 C 

53 C 

54 C 

55 

56 

57 


* **************************************************** 

* FILE :ACHBD. FOR * 

* SUBROUTINE FOR REM EFF-AXIAL FLOW CHARCOAL BED * 

* DOESNT ALLOW FOR DESORPTION * 

* ******************★**★★**★★★★**★★*★★★***********★*★* 


SUBROUTINE ACHBD (TN , TNI , CIIN , TCABIN, COEXIS , BEDQ , EMAX, BEDL , 
+BEDDIA , DENCH , TRTTYP , DCONT , VMOL , MW , VCONC , SOL , SMR , EFF , RH ) 

OUTPUT : 

EFF=BED REMOVAL EFF (DEC) 

INPUTS : 

TN, TNI = INCREMENT INITIAL AND FINAL TIMES (HR) 

Cl IN=BED INLET CONT CONC (MG/CU M) 

TC AB IN =C AB IN TEMP (DEG K) 

COEX I S=COEXI STANCE FACTOR 
BEDQ=BED FLOW RATE (CU M/HR) 

EMAX = MAX I MUM BED EFF (DEC) 

BEDL=BED LENGTH (M) 

BEDDIA=BED DIAMETER (M) 

DENCH=DENS ITY OF CHARCOAL IN BED (KG/CU M) 

TRTTYP=BED TREATMENT TYPE(1=CI CHAR,2=PHOS ACID, OTHER 
# =NONE ) 

DC ONT = CONT LIQUID DENSITY (GM/CC) 

VMOL = CONT MOLAR VOL (GM/CC) 

MW=CONT MOLECULAR WGT 

VCONC=CONT VAPOR CONCENTRATION AT TCABIN (MG/CU M) 

SOL=HENRY ' S LAW CONSTANT FOR WATER SOLUBILITY 
(ATM/MOL FRACTION) 

SMR=SUM OF CONT MASS STORED IN BED (MG) -FROM LAST INCR 

REAL LPREV, LAVNl , LUTIL, LIMM, LAVAV , LADS, MW 
INTEGER FACID, FCI 

SET CIN=CIIN (THIS PREVENTS CIN FROM BEING PASSED BACK UP 
TO OTHER SUBROUTINES IF IT IS SET TO IE-20) 

CIN=CIIN 

BED TREATMENT LOGIC 

FAC I D= FLAG IF BED IS TREATED WITH PHOSPHORIC ACID (Y=l 
N=0 ) 

FCI=FLAG FOR Cl CHAR IN BED (REMOVES FORMALDAH YDE ) 

IF (NINT (TRTTYP) .EQ. 2) THEN 
FACID=1 
FCI = 0 

ELSEIF (NINT (TRTTYP) .EQ. 1) THEN 
FACID=0 
FCI = 1 
ELSE 

FACID=0 
FCI = 0 
END IF 

TEST FOR NO BED FLOW (BEDQ=<0 ) OR 
TN-TN1<=0 ; BEDL , BEDDIA, DENCH=0 

IF ( (BEDQ.LE.0) .OR. (TN-TN1 . LE . 0 ) .OR. (BEDL.LE. 0) .OR. ( BEDDIA . LE . 0 ) 
+. OR. (DENCH. LE.0) ) THEN 
EFF=0 
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GOTO 199 
ENDIF 

TEST FOR Cl CHARCOAL AND FORMALDEHYDE ( FCI = 1 AND MW=30.03 
IF ( ( MW . EQ .30.03) .AND. (FCI.EQ.l) ) THEN 

CALL C ICH ( EFF , EMAX , BEDL , BEDDI A , DENCH , SMR , BEDQ ) 

GOTO 199 
ENDIF 

TEST FOR AMMONIA AND H3P04 ACID ON CHAR (FACID=1 AND 
MW=17 .0) 

IF ( (MW. EQ. 17.0) .AND. (FACID.EQ.l) ) THEN 

CALL ACIDCH ( EFF , EMAX , BEDL , BEDDIA , DENCH , SMR ) 

GOTO 199 
ENDIF 

TEST FOR MOL VOL=0 {NO CHAR REMOVAL) 

IF (VMOL.EQ.O) THEN 
EFF=0 
GOTO 199 
ENDIF 

CHARCOAL REMOVAL EFFICIENCY CALCULATION 
SUPERFICIAL BED VEL(FT/MIN) 

BEDVEL=BEDQ* .06960 /BEDDIA* *2 

TEST FOR CIN TOO SMALL IN AVAL CALC 

IF (CIN.LT. IE-20) CIN=lE-20 

AVAL= (TCABIN/VMOL) *LOG10 (VCONC/CIN) 

ADS ZONE LENGTH FOR 90% REMOVAL (M) 

LADS = AVAL* .000275* (BEDVEL/1 . 3) ** .8 
GET QI(CC LIQ CONT/GM CHAR) 

CALL FQI (AVAL, QI ( FACID, SOL, RH) 

LENGTH OF BED PREVIOUSLY USED BY CONT AT THIS C INLET (M) 
LPREV=SMR* 1 . 27 3E-6 *COEXIS/ (DCONT* DENCH* BEDDI A* *2 *QI ) 

RATE OF BED USAGE (M BED/ MG CONT) 

LIMM=1 . 2 73E- 6 *COEXIS / (DCONT*DENCH*BEDDIA**2*QI ) 

LENGTH OF BED AVAILABLE FOR ADS ZONE AT BEGINNING OF 
INCREMENT (M) 

LAVN1=BEDL-LPREV 
IF (LAVN1.LT.0) LAVN1=0 
FIX HERE IF DESORPTION IS DESIRED 
IF ( LAVN 1 / LADS . GT .20) THEN 
EFFAV=EMAX 
ELSE 

INIT INCR EFF BASED ON C IN AND BED L AVAIL AT BEG OF 
INCR (DEC) 

EFAVN1=EMAX* ( 1-EXP (-2 . 3025851*LAVNl/LADS) ) 

LOOP FOR EFFICIENCY 

EFFAV=EFAVN1 

DO 399 J=l,10,l 

LENGTH OF BED UTILIZ IN INCR (M) 

LUTIL=CIN*BEDQ*EFFAV* (TN-TN1) *LIMM 
IF (LUTIL . GT . LAVN1 ) THEN 
GOTO 299 
ELSE 

AVERAGE BED LENGTH AVAIL (M) 

LAVAV= LAVN l-LUTIL/2 
IF ( (LAVAV/LADS) .GE.20) THEN 
EFFAV=EMAX 
GOTO 299 
ELSE 
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118 C 

119 

120 
121 

122 399 

123 299 

124 C 

125 

126 

127 

128 

129 C 

130 

131 C 

132 199 

133 

134 

135 


AV EFF BASED ON AV BED L AVAIL (DEC) 
EFFAV=EMAX* (1-EXP (-2 . 3025851*LAVAV/LADS) ) 
ENDIF 
END IF 
CONTINUE 
ENDIF 

MAX EFF BASED ON C IN AND RATE OF BED USAGE (DEC) 
EFFMAX=LAVNl / (CIN*BEDQ* (TN-TNl ) *LIMM) 

IF ( EFFAV . GT . EFFMAX ) EFFAV=EFFMAX 
IF ( EFFAV. LT.O) EFFAV=0 
IF ( EFFAV . GT . EMAX ) EFFAV=EMAX 
EFF=ACTUAL EFF OUTPUT FROM SUBROUTINE 
EFF=EFFAV 

REMOVE THIS CHECK IF DESORPTION IS ADDED 

IF (EFF. LT.O) EFF=0 

IF ( EFF. GT. EMAX) EFF=EMAX 

RETURN 

END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 


136 C 

137 C 

138 C 

139 C 

140 

141 C 

142 C 

143 C 

144 C 

145 C 

146 C 

147 C 

148 C 

149 C 

150 C 

151 C 

152 

153 C 

154 

155 

156 

157 

158 

159 

160 C 

161 C 

162 

163 

164 


★ ★*★★**★★*★★******★★*★★★*★★*★ + ★******★*★**★*★★*★★★*★* 

* SUBROUTINE ACIDCH - CALCULATES REMOVAL EFF * 

* BED WITH NH3 AND 1.22 MILLIMOLE H3P04 ON CHAR * 

***************************************************** 

SUBROUTINE ACIDCH (EFF, EMAX , BEDL, BEDDIA, DENCH , SMR) 
OUTPUTS 

EFF=OUTPUT REMOVAL EFF (DEC) 

INPUTS 

EMAX = MAX I MUM BED REMOVAL EFF (DEC) 

BEDL=BED LENGTH (M) 

BEDDIA=BED DIAMETER (M) 

DENCH =CHARCOAL DENSITY (KG/CU M) 

SMR=SUM OF MASS OF CONT REMOVED AT BEG OF INCR (MG) 

FOR AMMONIA CAPACITY AT SMAC 
CHAR USED (KG) 

CHRUSD=1 . 6E-4 *SMR 
CHAR BED WGT(KG) 

BEDWGT=BEDL*BEDDIA* *2 * . 7 85 *DENCH 
IF (CHRUSD.LT. 0. 8 *BEDWGT) THEN 
EFF=EMAX 
ELSE 

EFF=EMAX*SIN( ( BEDWGT-CHRUSD) *1.57/ (BEDWGT*0.2) ) 
ENDIF 

PREVENTS NEGATIVE EFF FOR REACTION 

IF (EFF. LT.O) EFF=0 

IF (EFF. GT. EMAX) EFF=EMAX 

RETURN 

END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 


165 C 

166 C 

167 C 

168 C 

169 


* SUBROUTINE CICH - CALCULATES REMOVAL EFF * 

* FOR FORMALDELYDE AND Cl CHAR BED * 

* ★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★A-*** 

SUBROUTINE C ICH ( EFF , EMAX , BEDL , BEDDIA , DENCH , SMR , BEDQ ) 
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170 C 

171 C 

172 C 

173 C 

174 C 

175 C 

176 C 

177 C 

178 C 

179 C 

180 

181 C 

182 

183 

184 

185 

186 

187 

188 C 

189 

190 

191 

192 

193 C 

194 C 

195 

196 

197 

NUMBER 

NUMBER 

198 C 

199 C 

200 C 

201 C 

202 C 

203 

204 

205 C 

206 C 

207 C 

208 C 

209 C 

210 C 

211 C 

212 C 

213 

214 C 

215 

216 C 

217 

218 C 

219 

220 C 

221 
222 

223 C 

224 C 

225 C 


OUTPUTS 

EFF=OUTPUT REMOVAL EFF (DEC) 

INPUTS 

EMAX =MAX IMUM BED REMOVAL EFF (DEC) 

BEDL=BED LENGTH (M) 

BEDDIA=BED DIAMETER (M) 

DENCH=CHARCOAL DENSITY (KG/CU M) 

SMR=SUM OF MASS OF CONT REMOVED AT BEG OF INCR (MG) 
BEDQ=BED FLOW RATE (CU M/HR) 

BEDWGT=BEDL*BEDDIA**2* . 785*DENCH 
PERCENT OF BED WEIGHT CONSUMED (DEC) 

PBWGT=SMR/ (BEDWGT* 1E6 ) 

IF (PBWGT.LT. . 0012 ) THEN 
EFF=1-PBWGT*83 . 3 
ELSE 

EFF= . 9 *COS ( PBWGT* 1 . 5 7 / . 0 5 ) 

ENDIF 

BED RESIDENCE TIME (SEC) 

BREST=BEDL*BEDDIA* *2 *3600/ (BEDQ* 1 .273 ) 

IF ( BREST . LT . 0 . 2 5 ) THEN 
EFF= EFF* BREST/ . 2 5 
ENDIF 

PREVENTS NEGATIVE EFF FOR REACTION 

IF (EFF.LT.O) EFF=0 

IF (EFF. GT. EMAX) EFF=EMAX 

RETURN 

END 

OF WARNINGS IN PROGRAM UNIT: 0 
OF ERRORS IN PROGRAM UNIT: 0 

★ *★★★★★★★*★★★**★★*★***★*★★★★★★★★**★*★**■★**★★★★★★**★★★* 

* SUBROUTINE FQI - FINDS QI, THE ACID TREATED * 

* CHARCOAL CAPACITY FOR A CONTAMINANT AT 0 TO 75% RH * 

* (CC LIQ CONT/GM CHAR) * 

***************************************************** 

SUBROUTINE FQI (AVAL, QI , FACID, SOL, RH) 

INTEGER FAC ID 
OUTPUTS 

QI=CHARCOAL CAPACITY (CC LIQ CONT/GM CHAR) 

INPUTS 

AVAL= A VALUE OF CONTAMINANT 

FACID= FLAG FOR ACID TREATED CHAR IN BED (Y=l N=0) 

SOL=CONT AM INANT HENRY'S LAW CONSTANT (ATM/MOL FRACTION) 
RH=RELATIVE HUMIDITY (%) 

IF (RH.LT.0) RH=0 

CARBON CAPACITY DATA NOT AVAILABLE ABOVE 75% RH 
IF ( RH . GT .75) RH=75 

IF { AVAL . LT . 0 ) AVAL=0 
A VALUE .GT. 8 AND .LT. 200 

IF ((AVAL .GT. 8) .AND. (AVAL . LT . 200 ) ) THEN 
SOLUBLE CONTAMINANTS 

IF(SOL.GT.O. AND.SOL.LT. 5) THEN 
QI=2 . 1*EXP (-0 . 31*AVAL) 

INSOLUBLE CONTAMINANTS (IF HENRY'S LAW CONSTANT FOR A CONTAM. 
IS NOT AVAILABLE A 0 VALUE IS ASSIGNED AND CONTAMINANT IS 
CONSIDERED WATER INSOLUBLE) 
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226 

227 

228 

229 

230 

231 

232 

233 

234 C 

235 

236 

237 

238 

239 

240 

241 

242 C 

243 

244 

245 

246 

247 


ELSE 

IF (RH . LE . 50 ) THEN 

QI= (0 . 0 00096 *RH* *2 - 0 . 0 1 88 *RH + 2 .11) *EXP (-0.3 1 * AVAL ) 
ELSEIF (RH.GT.50) THEN 

QI= (0 . 000096* RH **2-0. 0188*RH+2 . 11} * 

+ EXP(-AVAL* (0.25+0. 0012*RH) ) 

END IF 
END IF 

A VALUE .LE. 8 

ELSEIF (AVAL . LE . 8 ) THEN 

IF ( SOL. GT.O .AND.SOL.LT. 5) THEN 
QI=0 . 5-AVAL*0 . 0405 
ELSE 

QI=- 0 . 0000128* RH **2-0. 00264* RH+0. 5+ (0. 00000112*RH**2+ 

+ 0. 000208* RH -0.0405) *AVAL 

END IF 

AVAL .GE. 200 
ELSE 

QI=lE-20 

ENDIF 

RETURN 

END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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RM/FORTRAN Compiler (V2.42) 

Source File: C:\RMFORT\TCC\ALIOH.F Options: /C 80 /L /BY 05/21/92 12:54:48 


1 C 

2 C 

3 C 

4 C 

5 C 

6 

7 

8 C 

9 C 

10 C 

11 C 

12 C 

13 C 

14 C 

15 C 

16 C 

17 C 

18 C 

19 C 

20 C 

21 C 

22 C 

23 C 

24 

25 

26 

27 

28 C 

29 

30 C 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 C 


★ ★★★★★★★★★a*******************************************’*********** 

* SUBROUTINE ALIOH - AXIAL FLOW BED * 

* CALCULATES REMOVAL EFFICIENCY FOR LI OH * 

★ *★*★★★*★**★★*★★★★★★*★★**★★★*★★★★★★***★*★*★★★★*★★★*★★*★******★★** 


SUBROUTINE ALIOH (TN, TNI , EMAX , BEDL , DFNLI , BEDDI A , RWUTLI , SWUTLI , 
+REMFCT, EFF) 


INPUTS : 

TN= FINAL INCREMENT TIME (HRS) 

TN1=INITIAL INCREMENT TIME (HRS) 

EMAX = MAX I MUM POSSIBLE REMOVAL EFFICIENCY (DEC) 

BEDL=BED LENGTH (M) 

DENLI=LIOH DENSITY (KG/CU M) 

BEDDIA=BED DIAMETER (M) 

RWUTLI =RATE OF LIOH USAGE FOR ALL CONTAMINANTS FROM LAST INCR (KG/HR) 
SWUTLI =SUM OF WEIGHT OF LIOH UTILIZED FROM LAST INCR (KG) 

REMFCT=L IOH REMOVAL FACTOR (LB LIOH /LB CONTAMINANT) 

OUTPUTS : 

EFF= REMOVAL EFFICIENCY (DEC) 


IF CONT DOESNT REACT WITH LIOH OR BEDL<=0 OR BED DIA <=0 OR 


DENLI <=0 THEN REM EFF =0 

IF ( (REMFCT.LE.O) .OR. (BEDL.LE.0) .OR. (BEDDI A. LE . 0 ) .OR. ( DENLI. LE.0) ) 
+ THEN 

EFF=0 

ELSE 

BED WEIGHT (KG) 

BEDWGT=BEDL* (BEDDIA) **2* . 785*DENLI 

TOTAL WEIGHT OF LIOH UTILIZED AT AVERAGE TIME IN INCREMENT (KG) 
TWUTL I = S WUTL I +RWUTL I * (TN-TNl ) /2 
IF ( TWUTL I / BEDWGT . LE .0.8) THEN 
EFF=EMAX 
ELSE 

EFF=EMAX * S IN ( ( BEDWGT-TWUTLI ) *1.57/ ( BEDWGT* 0 . 2 ) ) 

END IF 

IF (BEDL.LT. 0.0191) THEN 
EFF=EFF*BEDL/0 .0191 
END IF 
END IF 

IF(EFF.LT.O) EFF=0 

IF ( EFF. GT. EMAX) EFF=EMAX 

RETURN 

END 


*★★■*■*★***** 


END OF SUBROUTINE ALIOH **************************** 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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RM/ FORTRAN Compiler (V2.42) 

Source File: C:\RMFORT\TCC\CAFILL. Options: /C 80 /L /BY 05/21/92 12:54:53 

2 C * SUBROUTINE CAFILL 

3 C * SUBROUTINE TO FILL ADJUSTABLE SIZE CHAR ARRAY WITH BLANKS 

5 SUBROUTINE CAFILL (NN, NROW, NCOL) 

6 INTEGER NROW, NCOL 

7 CHARACTER NN ( NROW, NCOL) *30 


★ ★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★■A:******* 


★ ★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★•A-**** 


9 C NN=ARRAY NAME- ARRAY HAS 30 CHARACTERS 

10 C NROWrrNUMBER OF ROWS IN ARRAY (INTEGER) 

11 C NCOL=NUMBER OF COLUMNS IN ARRAY ( INTEGER) 12 

13 DO 110 1=1, NROW 

14 DO 100 J=1 , NCOL 

15 NN ( I , J ) = 1 1 

16 100 CONTINUE 

17 110 CONTINUE 

18 RETURN 

^ g q ************* END OF SUBROUTINE CAFILL ************************** 

2 0 END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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RM/ FORTRAN Compiler (V2.42) 

Source File: C:\RMFORT\TCC\CALCM.F Options: /C 80 /L /BY 05/21/92 12:54:56 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 


Q **************************************************************** 

C * SUBROUTINE CALCM * 

C * SUBROUTINE TO CALCULATE SUM OF MASSES REMOVED BY ALL DEVICES * 

C * USES CAV CABIN, REM EFF DD(J,20), & M.GEN DD(J,19) TO CALC * 

C * DEVICE CIN & COUT DD(J,17)& DD { J , 18 ) , M . REM CABIN+DEV DD(J,21)* 

Q *★*★*★*★***********★*★★★★★***★*★*★★★★★************************** 

SUBROUTINE CALCM ( DD , NROW, NCOL , CAV , SMGEN, SMREM , TN , LIN2 ) 

INTEGER NROW, NCOL, LIN2 
REAL DD (NROW, NCOL) 

C 

C SUBROUTINES REQUIRED: 

C NONE 

C INPUTS: 

C CAV=CABIN CONT AVERAGE CONCENTRATION (MG/CU M) 

C DD , NROW , NCOL=NAME AND SIZE OF MAT DD 
C TN= INCREMENT FINAL TIME (HRS) 

C LIN2=NO . OF DEVICES IN MAT DD 

C REM EFF AND DEVICE+CABIN M.GEN MUST BE LOADED INTO MAT DD 
C BEFORE USING THIS SUBROUTINE 

C OUTPUTS: 

C M.REM (MG/HR) FOR ALL DEVICES + CABIN CALCULATED AT CAV, 

C ARE STORED IN MAT DD COL 21 

C SMREM=SUM OF MASS OF CONT REM IN DEVICES (MG/HR) 

C SMGEN=SUM OF MASS GENERATED BY ALL DEVICES + CABIN (MG/HR) 

C 

C LOAD DEVICES 1 AND 2 WITH CIN AND COUT+M. REMOVED FOR DEV 2 
C DEVICE 1=CABIN 

C LOAD MAT DD WITH CAV CABIN 

DD ( 1 , 17) =CAV 
DD(1, 18) =CAV 
DD(1, 21) =0 
C DEVICE 2 

DD (2,17) =CAV 
DD (2 , 18) =CAV 

DD(2, 21) =DD{2 ,2 ) *DD(2 , 18) 

C FOR DEV 3-15 CALC CIN COUT AND M . REMOVED BY DEVICE 
C IF DEVICE FLOW RATE=0 THEN SET CIN, COUT, & M.GEN=0 

DO 100 J=3 , LIN2 
IF (DD{ J, 2 ) .EQ. 0) THEN 
DD(J, 17) =0 
DD(J, 18) =0 
DD ( J, 21) =0 
GOTO 100 
END IF 

C IF UPSTREAM DEVICE=1 OR 2 THEN SET INLET=CABIN CONC+DEV M.GEN/Q 
I F ( ( DD ( J , 4 ) . EQ . 1 ) . OR . ( DD ( J , 4 ) . EQ . 2 ) ) THEN 
DD ( J , 17 ) =DD (1,18) +DD ( J , 19 ) /DD ( J, 2 ) 

ELSE 

C DETERMINE FLOWS, CIN AND COUT FOR DEVICES WITH RELATIVE ADDRESSES 

IF ( DD ( J , 4 ) . EQ . 0 ) THEN 
QNOl=0 
CNOl=0 
ELSE 

QN01=DD (NINT ( DD ( J , 4 ) ) ,2) 

CN01 = DD (NINT (DD ( J, 4 ) ) , 18) 

ENDIF 
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59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 C 

74 C 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 C 

86 

87 

88 

89 C 

90 

91 C 

92 

93 C 

94 

95 

96 C 

97 C 

98 

99 
100 
101 

102 C 

103 

104 

105 

106 

107 C 

108 

109 

110 C 

111 


IF (DD(J, 5) .EQ.O) THEN 
QNO2=0 
CNO2=0 
ELSE 

QN02=DD(NINT(DD(J, 5) ) ,2) 

CN02 = DD (NINT (DD (J, 5) ) , 18) 

END IF 

IF (DD(J, 6) .EQ.O) THEN 
QNO3=0 
CNO3=0 
ELSE 

QN03=DD (NINT (DD ( J, 6 ) ) ,2) 

CN03=DD (NINT ( DD ( J , 6 ) ) , 18) 

ENDIF 

IF ALL UPSTREAM DEVICE FLOWS=0 
IF (QN01+QN02+QN03 .EQ.O) THEN 
DD ( J, 17 ) =0 
DD ( J, 18 ) =0 
DD ( J, 21 ) =0 

OPEN ( 2 , FILE= ' CON ' , IOSTAT=IOVAL) 

WRITE (*,*)' FLOW HALTED -UPSTREAM DEV TURNED OFF- 
+ INC END TIME; DEV= 1 , TN, DD ( J, 1 ) 

CLOSE (2) 

GO TO 100 
ELSE 

CALCULATE CIN 

DD ( J , 17) = (QNOl *CN01+QN02 *CN02 +QN03 *CN03 ) / 

+ ( QNOl +QN02 +QN03 ) +DD(J, 19) /DD(J, 2) 

ENDIF 

END OF DETERMINE FLOWS , CIN, COUT OF DEV WITH REL ADDR. 

ENDIF 

CALCULATE COUT 

DD ( J , 18) =DD { J , 17) * ( 1 -DD ( J, 20 ) ) 

CALCULATE SUM OF MASS REMOVED (CIN*Q*REM EFF) 

DD { J, 21 ) =DD ( J, 17)* DD ( J , 2 ) * DD ( J , 2 0 ) 

100 CONTINUE 

END OF LOADING OF MAT DD WITH DATA AND CALCULATING CIN COUT, M. REM 
SUM TOTAL MASS OF CONT REMOVED BY ALL DEVICES (2-15) (MG/HR) 

SMREM=0 

DO 101 J=2 , LIN2 
SMREM=SMREM+DD ( J , 2 1 ) 

101 CONTINUE 

SUM MASS OF CONT GENERATED IN ALL DEVICES+CABIN (1-15) (MG/HR) 
SMGEN=0 

DO 102 J=1 , LIN2 
SMGEN=SMGEN+DD ( J , 1 9 ) 

102 CONTINUE 

CALC M.REM CABIN AND PUT IN DD(1,21) 

DD (1,21) =SMGEN-SMREM 
RETURN 

************ END OF SUBROUTINE CALCM *************************** 
END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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RM/ FORTRAN Compiler (V2.42) 

Source File: C:\FORTRAN\TCC\CATBNR Options: /C 80 /L /BIJY 


03/15/94 10:39:12 


1 

2 C 

3 C 

4 C 

5 C 

6 C 

7 

8 C 

9 

10 C 

11 C 

12 C 

13 C 

14 C 

15 C 

16 C 

17 C 

18 C 

19 C 

20 
21 
22 

23 C 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 


* ★★★★★★★★a-*****************:************************************** 

* SUBROUTINE CATBNR * 

* CALCULATES EFFICIENCY OF CAT BURNER * 

***************************************************************** 

SUBROUTINE CATBNR ( P , EMAX , OXNEW , OXID , MW, EFF ) 

REAL MW 
INPUTS : 

P= TOTAL POISON (CONTAMINANTS IN NHB CATEGORIES 6, 7 AND 12: 
CHLOROCARBONS , CHLOROFLUOROCARBONS AND SULFIDES) REMOVED 
BY CAT BURNER (MG) 

EMAX = MAX I MUM BED EFFICIENCY (DEC) 

OXID=DEGREE OF OXIDIZATION OF CHEMICAL ( 1= FULLY, 0=NONE) 
MW=MOLECULAR WEIGHT 
OUTPUTS : 

EFF=REMOVAL EFF (DEC) 

IF(OXID.LT.O) OXID=0 
IF(OXID.GT.l) OXID=l 
EFF= EMAX *0X1 D 

EFFICIENCY FOR METHANE (MW=16.04) IS A FUNCTION OF P 
IF (MW. EQ. 16 . 04 ) THEN 
OXNEW=OXID 
IF (P.LE. 5500) THEN 

EFF=0 .97506*10** (-0. 00010507 *P) *EMAX*OXID 
ELSEIF (P.GT.5500) THEN 

EFF= (31.453-1. 151* 1E-3*P+1. 9046* 1E-8*P* *2 
+ -1.0389*1E-13*P**3 ) *0 . 01*EMAX*OXID 

ENDIF 
ELSE 

EFF=EMAX*OXID 
END IF 

IF(EFF.LE.O) EFF=0 

IF ( EFF. GT. EMAX) EFF=EMAX 

RETURN 

END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 


52 



RM/ FORTRAN Compiler (V2.42) 

Source File: C:\FORTRAN\TCC\CNRSUB Options: /C 80 /L /BIJY 03/15/94 10:39:14 


1 C 

2 C 

3 C 

4 C 

5 C 

6 


************************ 


***************************************** 


SUBROUTINE CNRSUB 

FOR 1 CONT AT A TIME AT CAV PRED, CALL EFF SUBROUTINES 
FOR DEVICES AND PUT EFFICIENCY IN MAT DD COL 20 


★ ★**************************************** ifc ,^^^^ iki ^^^ + ^ 5t 


******** 


* 

* 


* 

★ * * 


7 C NOTE: DEVICE NUMBERS IN THIS SUBROUTINE REFER TO DEVICE TYPES, 

8 C NOT THEIR RELATIVE POSITION IN MAT DD 

9 

10 SUBROUTINE CNRSUB ( I , TN , TNI , DD, NROW, NCOL, CC , NROW1 , NCOL1 , 

11 +CDI , NROW2 , NCOL2 , LIN2 , KK) 

12 INTEGER NROW, NCOL , NROW1 , NCOL 1 , NROW2 , NCOL2 , CATEG, D1 , D2 , D3 , KK 

13 REAL DD (NROW, NCOL) 

14 REAL CC (NROW1 , NCOL1 ) 

15 REAL CDI (NROW2 , NCOL2 ) 

16 REAL RGTM1 , RGTM2 , RGTM3 , P, POISN, OLDP 

17 C 

18 C SUBROUTINES REQUIRED: 

19 C ACHBD- REMOVAL EFF OF AXIAL CHARCOAL BED 

20 C RCHBD- REMOVAL EFF OF RADIAL CHARCOAL BED 

21 C AL I OH -REMOVAL EFF OF AXIAL LIOH BED 

22 C COOXID- REMOVAL EFF OF CO OXIDIZER 

23 C CATBNR- REMOVAL EFF OF CAT BURNER 

24 C CONDHX- REMOVAL EFF OF CONDENSING HX 

25 C 

26 C INPUTS: 

27 C FROM PCSET PREDCT , AND CONVRG 

28 C I =CONT NO. 

29 C TN,TNl=CONT INCREMENT FINAL, INIT TIME (HRS) 

30 C CAVPRD=CABIN AVERAGE CONT CONC (MG/CU M) 

31 C DD , NROW , NCOL=NAME & DIMENSIONS OF MAT DD 

32 C CC , NROW1 , NCOL1 =NAME & DIM OF MAT CC 

33 C CDI , NROW2 , NCOL2=NAME & DIM OF MAT CDI 

34 C LIN2 =NUMBER OF DEVICES IN MAT DD 

35 C DD (1,14) =PERCENT RELATIVE HUMIDITY 

36 C FROM EFFICIENCY SUBROUTINES 

37 C EFF =REMOVAL EFF (DEC) 

38 C OUTPUTS: 

39 C TO EFFICIENCY SUBROUTINES 

40 C DEVICE AND CONTAMINANT INFORMATION AS REQUIRED 

41 C DD ( J , 22 ) =DEVICE INLET CONCENTRATION 

42 C TO PCSET, PREDCT, & CONVRG 

43 C PUTS REMOVAL EFF FOR EACH DEVICE IN MAT DD COL 20 

44 C 

45 C EFFICIENCY FOR DEVICES 1 AND 2 

46 DD ( 1 , 2 0 ) =0 

47 DD (2,20) =DD (2,8) 

48 


49 C SET COUNTER FOR READING CHAR SMR IN MAT CC FOR DEVICE TYPE 3 

50 K=9 

51 C 

52 C BEGIN LOOP FOR DEVICES 3 TO 15 OF MAT DD 

53 DO 100 J=3 , LIN2 

54 C INDEX COUNTER 

55 K=K+3 

56 C SET REM EFF=0 AND GO TO END OF J LOOP IF DEVICE FLOW = 0 

57 IF(DD(J,2) .LE.1E-10) THEN 
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58 

59 

60 

61 C 
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DD(J, 20) =0 
GOTO 80 
END IF 

DECISIONS FOR VARIOUS DEVICES 
IF (NINT (DD ( J, 3 ) ) . EQ . 3 ) THEN 

GO TO SUBROUTINE FOR DEVICE TYPE 3 -AXIAL CHARCOAL BED 
CALL ACHBD ( TN , TNI , DD ( J , 2 2 ) , DD ( 1 , 10 ) , DD ( 1 , 13 ) , DD ( J, 2 ) , DD ( J , 8 ) , 

+ DD( J, 9) , DD ( J, 10 ) , DD ( J, 12 ) , DD( J, 13 ) , CDI (I, 2) ,CDI (1,3) , 

+ CDI (1,4) ,CDI(I, 5) ,CDI (1,6) ,CC(I,K) , EFF, DD ( 1 , 14 ) ) 

STORE EFF IN MAT DD 
DD ( J, 20 ) =EFF 

ELSEIF (NINT (DD ( J, 3 ) ) . EQ. 4 ) THEN 

GO TO SUBROUTINE FOR DEVICE TYPE 4 -RADIAL CHARCOAL BED 
CALL RCHBD(TN, TNI , DD ( J, 22 ) , DD (1 , 10) ,DD(1, 13 ) , DD ( J, 2 ) , DD ( J, 8 ) , 

+ DD ( J, 9) , DD ( J, 10 ) , DD (J, 11) , DD ( J, 12 ) , DD ( J, 13 ) , CDI (1,2) , CDI (1,3) , 
+ CDI (1,4) ,CDI (I, 5) ,CDI (1,6) , CC ( I , K) , EFF, DD ( 1 , 14 ) ) 

STORE EFF IN MAT DD 
DD ( J , 20) =EFF 

ELSEIF (NINT (DD ( J , 3 ) ) . EQ . 5 ) THEN 

GO TO SUBROUTINE FOR DEVICE TYPE 5-LIOH BED 
CALL ALIOH (TN,TN1 ,DD(J,8),DD(J,9),DD(J,10),DD(J,12),DD(J,15), 

+ DD ( J , 16) , CDI (1,7) , EFF ) 

STORE EFF IN MAT DD 
DD(J,20)=EFF 

ELSEIF ( NINT ( DD ( J , 3 ) ) . EQ . 6 ) THEN 

GO TO SUBROUTINE FOR DEVICE TYPE 6 -CO OXIDIZER 
CALL COOXID (DD ( J , 2 ) , DD ( J , 8 ) , DD { J , 9 ) , DD ( J , 10 ) , CDI (I, 4) , EFF) 
STORE EFF IN MAT DD 
DD ( J, 20 ) =EFF 

ELSEIF (NINT (DD ( J, 3 ) ) . EQ . 7 ) THEN 

SUM POISONS (CONTAMINANTS IN NHB CATEGORIES 6, 7 AND 12: 
CHLOROCARBONS , CHLOROFLUOROCARBONS AND SULFIDES) REMOVED 
BY CAT BURNER (MG) 

P = 0 . 0 
POISN=0 . 0 

DO 200 L=1 , NROW1 , 1 
CATEG=CDI (L, 8) 

IF ( (CATEG.EQ. 6) .OR. (CATEG.EQ.7) .OR. (CATEG.EQ. 12 ) ) THEN 
P=CC (L, ( J+l ) *3 ) 

POISN=POISN+P 
END IF 

200 CONTINUE 
RGTM1=0 
RGTM2 = 0 
RGTM3=0 
Dl=DD ( J , 4 ) 

D2=DD { J, 5) 

D3=DD) J, 6) 

IF ( (D1.NE.0) .OR. (D2.NE.0) .OR. (D3.NE.0) ) THEN 
IF ( (DD ( D1 , 3 ) . EQ . 3 ) . OR . (DD (D1 , 3 ) . EQ. 4 ) ) THEN 
TRCI=DD(D1, 15) 

TRD=DD (Dl , 16) 

TIR=DD(D1, 14) 

IF (TN1.EQ.0) GOTO 50 
IF (TN1.LT.TIR) GOTO 50 
IF (TRCI.LE.0) GOTO 50 

IF ( AINT ( (TN1-TIR) /TRCI ) . EQ . ( (TN1-TIR) /TRCI ) ) THEN 
RGTMl=l 
GOTO 55 
ENDIF 
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050 RGTM1=0 

055 CONTINUE 

END IF 

IF ( ( DD ( D2 , 3 ) . EQ . 3 ) . OR . ( DD ( D2 , 3 ) . EQ . 4 ) ) THEN 
TRCI =DD ( D2 ,15) 

TRD=DD ( D2 ,16) 

TIR=DD (D2 , 14 ) 

IF (TN1.EQ.0) GOTO 60 
IF (TN1.LT.TIR) GOTO 60 
IF (TRCI.LE.0) GOTO 60 

IF (AINT( (TNl-TIR) /TRCI) .EQ. ( (TN1-TIR) (/TRCI) ) THEN 
RGTM2=1 
GOTO 65 
ENDIF 

060 RGTM2 = 0 

065 CONTINUE 
ENDIF 

IF { ( DD { D3 , 3 ) . EQ . 3 ) . OR . ( DD ( D3 , 3 ) . EQ . 4 ) ) THEN 
TRCI=DD(D3 , 15) 

TRD=DD (D3 , 16) 

TIR=DD (D3 , 14 ) 

IF (TN1.EQ.0) GOTO 70 
IF (TN1.LT.TIR) GOTO 70 
IF (TRCI.LE.0) GOTO 70 

IF (AINT( (TNl-TIR) /TRCI) .EQ. ( (TNl-TIR) /TRCI) ) THEN 
RGTM3=1 
GOTO 75 
ENDIF 

070 RGTM3 =0 
075 CONTINUE 
ENDIF 
ENDIF 

IF ( ( RGTM1 . NE . 0 ) . OR . ( RGTM2 . NE . 0 ) . OR . ( RGTM3 . NE . 0 ) ) THEN 
OLDP=POISN 

REINITIALIZES METHANE OXIDATION EFFICIENCY TO 90% OF PREVIOUS 
IF ( (CDI (1,4) .EQ.16.04) .AND. (KK.EQ.l) ) THEN 
CDI (1,23) =0 .9*CDI (1,23) 

ENDIF 

ENDIF 

PO I SN= PO I SN - OLDP 

GO TO SUBROUTINE FOR DEVICE TYPE 7-CAT BURNER 
CALL C ATBNR ( PO I SN , DD ( J , 8 ) ,DD(J,9) , CDI (I, 23) ,CDI(I,4) , EFF ) 

STORE EFF IN MAT DD 
DD ( J , 2 0 ) =EFF 

ELSEIF (NINT ( DD ( J , 3 ) ) . EQ . 8 ) THEN 

GO TO SUBROUTINE FOR DEVICE TYPE 8 -CONDENSING HX 
CALL CONDHX ( DD ( J , 2 ) , DD ( J, 8 ) , DD ( J , 9 ) , CDI (1,4) ,CDI(I,5) ,CDI(I,6) 
+ DD ( J , 2 2 ) , EFF ) 

STORE EFF IN MAT DD 
DD(J, 20) =EFF 

ELSEIF (NINT(DD(J, 3) ) .EQ. 9) THEN 

SUBROUTINE FOR DEVICE TYPE 9 -DUMMY 
EFF=0 

STORE EFF IN MAT DD 
DD ( J, 20 ) =EFF 
ELSE 

DEVICES >9 OR <1 (OR ANY DEVICE NOT IN ABOVE CASES) 

DD ( J , 20 ) =0 
ENDIF 

080 CONTINUE 
100 CONTINUE 
RETURN 

END 


END OF SUBROUTINE CNRSUB ************************** 



NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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C ***************************************************************** 

C * SUBROUTINE CONDHX * 

C * CALCULATES EFFICIENCY OF CONDENSING HX * 

c ***************************************************************** 

c 

SUBROUTINE CONDHX ( BEDQ , EMAX , MLIQ , MW, VCONC , SOL , CAVPRD , EFF ) 

REAL MW, MLIQ, H, PA, XA, NOUT 

C 

C INPUTS: 

C BEDQ=BED FLOW RATE (M3 /HR) 

C EMAX=MAXIMUM POSSIBLE REMOVAL EFFICIENCY (DEC) 

C MLIQ=WATER FLOW RATE IN HX DUE TO CONDENSING (KG/HR) 

C MW=MOLECULAR WEIGHT OF CONTAMINANT 
C VCONT= VAPOR CONCENTRATION OF CONTAMINANT (MG/M3 ) 

C SOL=HENRY ' S LAW COEFFICIENT (ATM/MOL FRACTION) 

C OUTPUTS: 

C EFF=REMOVAL EFF (DEC) 

C 


C 

c 


c 

c 


c 


CAIN=CAVPRD 

IF (CAIN. LE . IE- 10 ) THEN 
CAIN=0 . IE- 10 
END IF 

IF (SOL.LE. IE-10) THEN 
EFF=0 
ELSE 

IF CONTAMINANT IS AMMONIA - USES EXPERIMENTAL DATA FROM JSC- 08797 
FOR AMMONIA REMOVAL AS A FUNCTION OF C02 CONCENTRATION (8/23/76) 

IF (MW. EQ. 17.0) THEN 

CAOUT= ( (CAIN*BEDQ) - (MLIQ* 18 9 . 5847418 *CAIN* * 0 . 53 4 9 152 56 ) ) /BEDQ 
EFF= ( (CAIN-CAOUT) /CAIN) *EMAX 
ELSE 

CONTAMINANT IS NOT AMMONIA 

CALCULATE CONTAMINANT PARTIAL PRESSURE AND WATER MOLE FRACTION 
PA=CAIN*1 . 0E- 9 * 82 . 06 *27 8 /MW 

XA= (PA/1) / ( (MLIQ* (1000/18) ) / ( BEDQ* ( 1000 /22 . 4 ) ) +SOL/1) 

NOUT=ML IQ * XA *1000/18 

CAOUT= ( (CAIN*BEDQ) - (NOUT*MW*1000 ) ) /BEDQ 
EFF = ( (CAIN-CAOUT) /CAIN) *EMAX 
END IF 
END IF 

IF(EFF.LE.O) EFF=0 

IF (EFF. GT. EMAX) EFF=EMAX 

RETURN 

END 

***************** END OF SUBROUTINE CONDHX *********************** 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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■**;;*r*r;;;;;;;;r;;;*;*^;;;*.**********-***********-^**-**** 

* SUBROUTINE CONVRG 

* MAIN CONVERGENCE LOOP SUBROUTINE 

* USING CAV PRED & BASED ON SUM MASS REM OF LAST INCR, CALC 

* NEW REM EFF, CAV CALC, CEQULIB , C FINAL, & M. REMOVED 

* PUT THEM IN MAT DD - WORKS FOR ONE CONT AT A TIME 
***************************************************************** 

SUBROUTINE CONVRG ( I , TN , TNI , CAVPRD , DD, NROW, NCOL , 

+CC , NROWl , NCOLl , CDI , NROW2 , NCOL2 , CAVCLC , CFINAL , CEQLIB , KK , LIN , LIN2 , 
+NN , PRTSW3 , IMSGDN ) 

INTEGER NROW, NCOL, NROWl , NCOLl ,NROW2 ,NCOL2 , PRTSW3 , KK 
CHARACTER NN( NROWl) *30 
REAL DD (NROW, NCOL) 

REAL CC (NROWl , NCOLl ) 

REAL CDI ( NROW2 , NCOL 2 ) 

REAL XI, X2,Y1,Y2, SLOPE 

SUBROUTINES REQUIRED: 

PRAFIL-ZERO MAT DD COL 17-21 . „ „ A 

CNRSUB -USING CAV PRED CALC- REM EFF FOR ALL DEVICES-PUT IN OD COL 20 
MASBAL-CALC CAV CALC , CFINAL , CEQULIB, M. REM 
INPUTS : 

FROM MCALC 

I =CONTAM INANT NO. 

TN, TNI = INCREMENT END & BEGINNING TIME (HRS) 

DD , NROW , NCOL^NAME & DIM OF MAT DD 
CC, NROWl, NCOLl=NAME & DIM OF MAT CC 
CDI , NROW2 , NCOL2 =NAME & DIM OF MAT CDI 
LIN=NO. OF CONTAMINANTS IN MAT CDI 
LIN2 =NO . DEVICES IN MAT DD 
IMSGDN=DEVICE NO. FOR MESSAGE OUTPUT 

CAVPRD= PREDICTED CABIN CONC FOR INCREMENT (MG/CU M) 

FROM PRAFIL 

ZEROS IN MAT DD COL 17-21 

CALL^REM EFF SUBROUTINE & PUTS REM EFF FOR EACH DEV IN DD COL 20 
FROM MASBAL 

CAVCLC = CALC CABIN CONT CONC (MG/CU M) 

CFINAL= FINAL INCR CABIN CONT CONC (MG /CUM) 

CEQLIB^ EQUILIBRIUM CABIN CONT CONC (MG/CU M) 

OUTPUTS : 

TO MCALC 

CAVCLC =CALC CABIN CONT CONC (MG/CU M) 

CFINAL= FINAL INCR CABIN CONT CONC (MG/CUM) 

KK=COUNTER FOR CONVERGENCE 

CEQLIB= EQUILIBRIUM CABIN CONT CONC (MG/CU M) 

M.REM IS IN MAT DD COL 21 
TO PRAFIL 

NAME & SIZE OF MATRIX + FIRST AND LAST COL TO BE ZEROED 
TO CNRSUB 

I =CONT NO. 

TN, TNl=CONT INCREMENT FINAL, INIT TIME (HRS) 

DEVICE AVERAGE CONT CONC (MG/M3) = DD(J,22) 

DD, NROW, NCOL =NAME & DIMENSIONS OF MAT DD 
C C , NROW 1 , NC OL 1 =NAME & DIM OF MAT CC 
CDI , NROW2 , NCOL2 =NAME & DIM OF MAT CDI 
LIN2 =NUMBER OF DEVICES IN MAT DD 
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C TO MASBAL 

C I 

C TN , TNI 

C CVOL=CABIN VOL (CU M) =DD ( 1 , 9 ) 

C CINIT=INCR INIT CABIN CONT CONC (MG/CU M)=CC(I,1) 

C 

C CONVERGENCE ERROR (DEC) 

CNVERR=DD ( 1 , 12) 

IF (CNVERR.LT. IE-10) THEN 

OPEN (IMSGDN, FILE= ' CON’ , IOSTAT=IOVAL ) 

WRITE ( IMSGDN, *) ' CONV ERROR<1E-10 : PROGRAM TERMINATED' 

CLOSE (IMSGDN) 

STOP 
END IF 
KK=1 

DO 100 KK=1 , 20 

C ZERO MAT DD COL 17-21 

CALL PRAFIL (DD, NROW, NCOL , 17,21) 

C USING CAVPRD FIND REM EFF OF EACH DEV & PUT IN DD COL 20 

CALL CNRSUB ( I , TN , TNI , DD , NROW , NCOL , CC , NROW1 , NCOL 1 , 

+ CDI , NROW2 , NCOL2 , LIN2 , KK ) 

C FIND CAVCLC FOR THESE REMOVAL EFFICIENCIES 

CALL MASBAL ( I , TN, TNI , DD, NROW, NCOL , CC, NROW1 , NCOL1 , 

+ CAVCLC , CDI , NROW2 , NCOL2 , CFINAL, CEQLIB , LIN, LIN2 ) 

C IF PRTSW3=1 THEN PRINT NAME & NO + CONV VALUES 

IF ( PRTSW3 . EQ . 1 ) THEN 
OPEN ( IMSGDN, FILE= ' CON' , IOSTAT=IOVAL) 

WRITE (IMSGDN, *) 'PRINTOUT FOR CONVERGENCE VALUES IN CONVRG 1 
WRITE (IMSGDN, 50) I , NN ( I ) 

050 FORMAT (IX, 'CONT NO. = ' ,I4,2X,A) 

WRITE (IMSGDN, *) ' CAVPRD , CAVCLC = ' , CAVPRD , CAVCLC 
CLOSE ( IMSGDN) 

END IF 

IF CAVCLC=CAVPRD THEN EXIT THE KK LOOP 
IF (CAVCLC. EQ. CAVPRD) GOTO 101 
IF CAVPRD<1E-10 THEN SKIP CONVERGENCE STEP 
IF (CAVPRD.LT. IE-10) GOTO 80 
IF CONVERGENCE IS REACHED EXIT THE KK LOOP 

IF ( ABS ( (CAVPRD- CAVCLC) /CAVPRD) . LT . CNVERR) THEN 
GOTO 101 
END IF 

: CONVERGENCE CALCULATION ROUTINE 

: USE THE BISECTION METHOD FOR THE ITERATION WHERE KK=1 

080 IF (KK.EQ.l) THEN 

: INITIALIZE X2 AND Y2 FOR THE NEXT ITERATION 

X2 =CAVPRD 
Y2=CAVCLC- CAVPRD 
CAVPRD= (CAVPRD+CAVCLC ) /2 
ELSE 

: USE THE NEWTON- RAPHSON METHOD FOR ITERATIONS WHERE KK>1 

Xl = X2 
Yl=Y2 
X2=CAVPRD 
Y2=CAVCLC -CAVPRD 
SLOPE= ( Y2-Y1 ) / (X2-X1) 

CAVPRD=X2-0 . 95*Y2 /SLOPE 
ENDIF 

: SET CAV IN PRED DD(I,22) = CAV IN CALC DD(I,17) 
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119 


DO 90 J=1 , LIN2 

120 


DD ( J , 22 ) =DD ( J, 17) 

121 

090 

CONTINUE 

122 

C 


123 

100 

CONTINUE 

124 

c 

END OF KK LOOP 

125 

101 

CONTINUE 

126 

c 


127 



128 


RETURN 

129 

c 

************ END OF SUBROUTINE CONVRG 

130 


END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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q ***************************************************************** 

C * SUBROUTINE COOXID * 

C * CALCULATES EFFICIENCY OF CO OXIDIZER ( Pt on charcoal) * 

Q ***************************************************************** 

c 

SUBROUTINE COOXID {BEDQ, EMAX, BEDL, BEDDIA, MW, EFF) 

REAL MW 

C 

C INPUTS: 

C BEDQ=BED FLOW RATE (CU M/HR) 

C EMAX=MAXIMUM POSSIBLE REMOVAL EFFICIENCY (DEC) 

C BEDL=BED LENGTH (M) 

C BEDDIA=BED DIAMETER (M) 

C MW=MOLECULAR WEIGHT OF CONTAMINANT 

C OUTPUTS: 

C EFF=REMOVAL EFF (DEC) 

C 

C WORKS ONLY FOR CO MW=28.01 OR H2=2.02; OTHERWISE REM EFF=0 

IF ( (MW.EQ.28.01) .OR. (MW.EQ.2 .02) ) THEN 
EFF=EMAX 

C IF RESIDENCE TIME <0.2 SEC THEN REM EFF DROPS LINEARLY 

C BREST = BED RESIDENCE TIME (SEC) 

BREST= (3.141592654/4 ) * BEDL* BEDDIA* *2 *3 600 /BEDQ 
IF ( BREST.LT. 0.2) THEN 
EFF = EMAX* BREST /0 . 2 
END IF 
ELSE 

C REM EFF FOR OTHER THAN CO OR H2 = ZERO 

EFF = 0 
END IF 

IF(EFF.LE.O) EFF=0 

IF ( EFF. GT. EMAX) EFF=EMAX 

RETURN 

END 


C 


***************** 


END OF SUBROUTINE COOXID *********************** 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 

NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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****************************************************************** 

* SUBROUTINE CRIN * 

* SUBROUTINE TO READ STRING OF LENGTH 30 INTO MAT NN * 

* AND READ REAL DATA INTO MAT XX (ROW, COL) * 

* RETURNS NUMBER OF LINES OF DATA READ FROM FILE * 

****************************************************************** 

NOTE: INPUT STRING MUST HAVE SINGLE QUOTES AROUND IT 
NOTE: INPUT NUMBERS MUST BE SEPARATED BY BLANKS 
SUBROUTINE CRIN (NN , XX , NROW, NCOL , LIN) 

INTEGER NROW , NCOL , IOVAL , LIN 
CHARACTER NN (NROW) *30 , FNAME*24 
REAL XX (NROW, NCOL) 

READ ( * , ' (A) ' ) FNAME 

OPEN ( 1 , FILE=FNAME , STATUS= ' OLD ' , IOSTAT=IOVAL ) 

I F ( IOVAL . NE . 0 ) GOTO 900 
LIN=0 

DO 100 1=1, NROW 

READ ( 1 , * , IOSTAT=IOVAL , END=500 , ERR=900 ) NN ( I ) , (XX ( I , J) , J=1 , NCOL ) 
LIN=LIN+1 
CONTINUE 

WRITE (*, ' (A) ' ) ' DONE WITH FILE INPUT' 

WRITE ( * , * ) ' ' 

CLOSE (1) 

GOTO 990 

WRITE ( * , * ) ' IOERROR= ' , IOVAL 
CLOSE (1) 

WRITE (*,*) 'WHAT IS THE INPUT FILE NAME? ' 

GOTO 10 
RETURN 

*★****★***★**★★ end op SUBROUTINE CRIN ************************ 
END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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RM/ FORTRAN Compiler (V2.42) 

Source File: C:\RMFORT\TCC\CROUT.F Options: /C 80 /L /BY 05/21/92 12:56:01 

1 C ***************************************************************** 

2 C * SUBROUTINE CROUT 

3 c * SUBROUTINE TO WRITE DATA TO CONSOLE, OR PRINTER 

4 C * WRITES STRG OF LENGTH 30 FROM MAT NN & REAL DATA FROM MAT 

5 c * XX (ROW, COL) STARTING WITH COL FSTCOL, AND ENDING WITH LSTCOL 

6 C * AND FROM LINE FSTLIN TO LINE LSTLIN 

7 C ***************************************************************** 

8 SUBROUTINE CROUT (NN, XX, NROW, NCOL, FSTCOL, LSTCOL, LIN, FSTLIN LSTLIN 

9 +IMSGDN, NINC , FNAME , IDEVNO, IOVAL) 

10 INTEGER 

NROW, NCOL, IOVAL, FSTCOL, LSTCOL, LIN, FSTLIN, LSTLIN, NINC, IDEVNO 

11 + IOVAL 

12 CHARACTER FNAME * 2 4 , DES * 1 

13 CHARACTER NN (NROW) *30- 

14 REAL XX (NROW, NCOL) 

15 IF ( FSTCOL. GT. NCOL) FSTCOL=NCOL 

16 IF (LSTCOL. GT. NCOL) LSTCOL=NCOL 

17 IF (FSTCOL. GT. LSTCOL) FSTCOL=LSTCOL 

18 IF (FSTLIN. GT. LIN) FSTLIN=LIN 

19 IF (LSTLIN. GT. LIN) LSTLIN=LIN 

20 IF (FSTLIN. GT. LSTLIN) FSTLIN=LSTLIN 


OR CON OR 


GOTO 


NINC 


0 OPEN ( IMSGDN, FILE=' CON' , IOSTAT=IOVAL) 

WRITE (IMSGDN, * (A) ' ) 1 WRITE TO LPT1 OR CON OR END ' 

CLOSE (IMSGDN) 

READ ( * , 1 (A) ‘ ) FNAME 

QUIT IF FNAME = END 

IF ( FNAME. EQ. 'END' ) GO TO 990 

IF( (FNAME. NE. 'LPT1 ') .AND. (FNAME. NE. 'CON' ) ) GOTO 10 
OPEN ( 1 , FILE=FNAME, IOSTAT=IOVAL) 

I F ( I OVAL . NE . 0 ) GOTO 900 

WRITE (IDEVNO, 55, IOSTAT= IOVAL, ERR=900) NINC 
FORMAT ( 1 INCREMENT NO. = ',17) 

DO 110 I=FSTLIN, LSTLIN 

WRITE (1, 60, IOSTAT=IOVAL, ERR=900) I,NN(I) 

WRITE ( IDEVNO, 60, IOSTAT=IOVAL , ERR=9 00 ) I,NN(I) 

FORMAT ( IX, 'CONT NO. = ',I4,2X,A) 

WRITE ( 1 , 7 0 , IOSTAT= IOVAL , ERR= 900) (XX(I,J), J=FSTCOL , LSTCOL ) 

WRITE (IDEVNO, 70, IOSTAT=IOVAL,ERR=900) (XX (I, J) , J=FSTCOL , LSTCOL ) 
FORMAT ( IX , 7G11 . 4 ) 

CONTINUE 
CLOSE ( 1 ) 

GOTO 990 

OPEN (IMSGDN, FILE= 1 CON’ , IOSTAT=IOVAL) 

WRITE ( IMSGDN, * ) ' IOERROR= ' , IOVAL 
CLOSE ( IMSGDN) 

CLOSE ( 1 ) 

CLOSE (IDEVNO) 

GOTO 10 
RETURN 

********** END q F SUBROUTINE CROUT ****************************** 


★ ★**★★**** 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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RM/FORTRAN Compiler (V2.42) 

source File: C:\RMFORT\TCC\CROUT2. Options: /C 80 /L /BY 


05/21/92 12:56:09 


1 C 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 


C 

C 

c 

c 

c 

c 


***************************************************************** 

* SUBROUTINE CROUT2 * 

* SUBROUTINE TO WRITE DATA TO CONSOLE, OR PRINTER 

* WRITES STRG OF LENGTH 30 FROM MAT NN & REAL DATA FROM MAT 

* XX (ROW, COL) STARTING WITH COL FSTCOL, AND ENDING WITH LSTCOL 

* AND FROM LINE F STL IN TO LINE LSTLIN 
****************************************************************** 
SUBROUTINE CROUT2 (NN, XX , NROW, NCOL , FSTCOL, LSTCOL, LIN, FSTLIN, LSTLIN, 

+IMSGDN) 

INTEGER NROW, NCOL, I OVAL , FSTCOL , LSTCOL , LIN , FSTLIN , LSTLIN 
CHARACTER FNAME*24 , DES* 1 
CHARACTER NN(NROW)*30 
REAL XX (NROW, NCOL) 

IF (FSTCOL. GT. NCOL) FSTCOL=NCOL 
IF (LSTCOL. GT. NCOL) LSTCOL=NCOL 
IF (FSTCOL. GT. LSTCOL) FSTCOL=LSTCOL 
IF (FSTLIN. GT. LIN) FSTLIN=LIN 
IF (LSTLIN. GT. LIN) LSTLIN=LIN 
IF (FSTLIN. GT. LSTLIN) FSTLIN=LSTLIN 


010 


070 

110 


900 


990 


IOSTAT=IOVAL ) 
WRITE TO LPT1 


OR CON OR END 


■CON*)) GOTO 10 


OPEN ( IMSGDN, FILE= ’ CON 
WRITE (IMSGDN, 1 (A ) ‘ ) 

CLOSE (IMSGDN) 

READ ( * , 1 (A) ' ) FNAME 
QUIT IF FNAME = END 
IF (FNAME. EQ. ' END' ) GO TO 990 
IF ( (FNAME .NE. ' LPTl 1 ) .AND. (FNAME. NE. 

OPEN ( 1 , FILE=FNAME , IOSTAT=IOVAL) 

I F ( IOVAL . NE . 0 ) GOTO 900 
DO 110 I=FSTLIN, LSTLIN 

WRITE (1,70, IOSTAT=IOVAL , ERR=900 ) (XX (I , J) , J=FSTCOL , LSTCOL ) 
FORMAT ( IX, 7G1 1.4) 

CONTINUE 
CLOSE (1) 

GOTO 990 

OPEN (IMSGDN, FILE=' CON' , IOSTAT= IOVAL) 


WRITE (IMSGDN, *) 
CLOSE (IMSGDN) 
CLOSE (1) 

GOTO 10 
RETURN 

********** end 

END 


• IOERROR= ' , IOVAL 


OF SUBROUTINE GROUT 


****************************** 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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RM/ FORTRAN Compiler (V2.42) 

Source File: C:\RMFORT\TCC\DATOUT. Options: /C 80 /L /BY 05/21/92 12:56:16 

1 C **************************************************************** 

2 C * SUBROUTINE DATOUT * 

3 C * SUBROUTINE TO PRINT HEADINGS AND DATA TO PRINTER, CON, OR FILE * 

4 C **************************************************************** 

5 C NOTES: (l)FILE MUST BE OPEN BEFORE CALLING THIS SUBROUTINE 

6 C (2) IDEVNO MUST BE 6 FOR FORM FEEDS TO BE PRINTED 

7 


8 


9 


10 


11 


12 


13 


14 


15 

c 

16 

c 

17 

c 

18 

c 

19 

c 

20 

c 

21 

c 

22 


23 

c 

24 

c 

25 

c 

26 

c 

27 

c 

28 

c 

29 

c 

30 

c 

31 

c 

32 

c 

33 

c 

34 

c 

35 

c 

36 

c 

37 

c 

38 

c 

39 


40 

c 

41 

c 

42 


43 

★ 

44 

★ 

45 


46 

c 

47 


48 


49 


50 

c 

51 


52 


53 


54 


55 


56 


57 


58 

c 


SUBROUTINE DATOUT ( TN , TNI , LIN , DD , NROW, NCOL , CC , NROW1 , NCOL1 , 

+CDI , NROW2 , NCOL2 , LIN2 , NN, PRTSW6 , PRTSW8 , PRTSW9 , 

+ IDEVNO , NINC , IMONTH, IDAY, IYEAR , IHOUR, IMINUTE , FNAME , IOVAL, IPGCTR, 
+TVAL , FCPLOT , IDEVN1 , IDEVN3 , IDEVN2 ) 

INTEGER PRTSW6 , PRTSW8 , PRTSW9 , TVAL , NINC 

SUBROUTINES REQUIRED: 

PRCDA= PRINT OUT OF CONTAMINANT DATA 

PRREM1= PRINTOUT OF RATE OF CONTAMINANT REMOVAL (MG/HR) -SHEET1 
PRREM2= PRINTOUT OF RATE OF CONTAMINANT REMOVAL (MG/HR) -SHEET2 
PRMAS1= PR INTOUT OF SUM MASSES REMOVED BY DEVICES (MG) -SHEET1 
PRMAS2 = PRINTOUT OF SUM MASSES REMOVED BY DEVICES (MG) -SHEET2 
PREFF= PRINTOUT OF INCREMENT END REMOVAL EFFICIENCIES 

INPUTS FROM MAIN PROGRAM: 

TN= INCREMENT FINAL TIME (HRS) 

TNI = INCREMENT INITIAL TIME (HRS) 

LIN=NO . OF CONT IN MAT CC A D NN 
DD, NROW, NCOL =NAME & SIZE OF MAT DD 
CC , NROWl , NCOLl=NAME & SIZE F MAT CC 
CDI , NROW2 , NCOL2 =NAME & SIZE OF MAT CDI 
LIN2=NO. DEVICES IN MAT DD 
NN=NAME OF MAT NN 
NINC=TIME INCREMENT NUMBER 

= 0 THEN PRINT HEADINGS & DATA FOR PRECALCULATION SET UP ROUTINE 
=-l THEN PRINT HEADINGS & DATA FOR FINAL ANSWERS 
ELSE PRINT WITH PROPER INCREMENT NUMBER 
IDEVNO=OUTPUT DEVICE NUMBER (SHOULD BE 6) 

IMONTH. . I HOUR=DATE AND TIME VARIABLES 

I PGCTR= PAGE COUNTER FOR SEQUENTIAL PAGE NO.'S ON ALL PAGES 

OUTPUT TO MAIN PROG: 

IOVAL=STATUS OF IOERROR IN SUBROUTINES 

TEST CASE ********** 

LIN=120 

PRINT OUT CONCENTRATION DATA 

CALL PRCDA ( TN, TNI , LIN , CC , NROWl , NCOL1 , CDI , NROW2 , NCOL 2 , NN , 

+ IDEVNO, NINC, IMONTH, IDAY, IYEAR, IHOUR, IMINUTE , FNAME , IOVAL, IPGCTR, 

+ PRTSW8 , PRTSW9 , FCPLOT, IDEVN1 ) 

PRINT OUT NHB 8060.1 GROUP CONTRIBUTION VALUES (T-VALUES) 

IF ( ( TVAL . EQ . 1 ) . OR . ( TVAL . EQ . 2 ) ) THEN 

CALL GROUP ( TN , TNI , L IN , CC , NROWl , NCOL1 , CDI , NROW 2 , NCOL2 , NN , 

+ IDEVNO, NINC, IMONTH, IDAY, IYEAR, IHOUR, IMINUTE , FNAME , IOVAL, 

+ IPGCTR, TVAL , IDEVN3 , PRTSW8 ) 

END IF 

IF (PRTSW6. EQ. 1) THEN 

IF ( ( PRTSW8 . EQ . 1 ) .OR. ( ( PRTSW8 . EQ . 0 ) .AND. (NINC.EQ.-l) ) ) THEN 
PRINTOUT OF RATE OF CONTAMINANT REMOVAL BY DEVICES-SHEETl 
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59 

60 

61 C 

62 

63 C 

64 

65 

66 

67 C 

68 

69 

70 C 

71 

72 C 

73 

74 

75 

76 

77 C 

78 

79 

80 
81 
82 

83 

84 


CALL PRREMl (TN, TNI , LIN, CC , NROW1 , NCOLl , CDI , NROW2 , NCOL2 , NN, 

IDEVNO , NINC , I MONTH , IDAY , I YEAR , IHOUR , IMINUTE , FNAME , IOVAL, IPGCTR} 
DON'T PRINT SHEET 2 UNLESS NUMBER OF DEVICES IN MAT DD > 8 
IF (LIN2.GT.8) THEN 

PRINTOUT OF RATE OF CONTAMINANT REMOVAL BY DEVICES- SHEET2 
CALL PRREM2 (TN , TNI , LIN, CC , NROWl , NCOLl , CDI , NROW2 , NCOL2 , NN, 
IDEVNO, NINC, IMONTH, IDAY, I YEAR, IHOUR, IMINUTE , FNAME , IOVAL, IPGCTR) 
ENDIF 

PRINTOUT OF SUM OF MASS REMOVED BY DEVICES- SHEETl 
CALL PRMAS 1 ( TN , TNI , LIN , CC , NROWl , NCOLl , CDI , NROW2 , NCOL2 , NN , 
IDEVNO, NINC, IMONTH, IDAY, I YEAR, IHOUR, IMINUTE , FNAME , IOVAL, IPGCTR) 
DON'T PRINT SHEET 2 UNLESS NUMBER OF DEVICES IN MAT DD > 8 
IF (LIN2 .GT.8) THEN 

PRINTOUT OF SUM OF MASS REMOVED BY DEVICES- SHEET2 
CALL PRMAS2 (TN , TNI , LIN , CC , NROWl , NCOLl , CDI , NROW2 , NCOL2 , NN , 
IDEVNO, NINC, IMONTH, IDAY, IYEAR, IHOUR, IMINUTE, FNAME, IOVAL, IPGCTR) 


ENDIF 

ENDIF 

PRINTOUT OF INCREMENT END REMOVAL EFFICIENCIES 
CALL PREFF (TN, TNI , LIN, CC , NROWl , NCOLl , CDI , NROW2 , NCOL2 , NN, 
+ IDEVNO, NINC, IMONTH, IDAY, IYEAR, IHOUR , IMINUTE , FNAME , IOVAL, 

+ I PGCTR , PRTSW8 , PRTSW9 , IDEVN2 ) 

ENDIF 


RETURN 

END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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c 

c 




25 


30 


C 


05/21/92 12:56:54 


C 


50 


r********************************^^^^^^^^^^^ 

* SUBROUTINE GROUP * 

* PROGRAM TO PRINT THE GROUP TOXICITY LEVELS AND T LEVEL 

******************************* i ,***i,*i,. ki , ir i :i , iti , i , i 'i, i , i ' i , i ' i ' itiri ' 1l . iri ' iri ' i ' i(i ' 

SUBROUTINE GROUP (TN , TNI , LIN , CC , NROW1 , NCOL1 , CDI , NROW2 , NCOL2 , NN 
+IDEVNO, NINC , IMONTH, IDAY, IYEAR, IHOUR, IMINUTE , FNAME , IOVAL 
+ I PGCTR, TVAL, IDEVNT, PRTSW8) 

DIMENSION GL (16) 

REAL CC ( NROW1 , NCOL1 ) 

REAL CDI (NROW2 , NCOL2 ) 

REAL TLEVL 

INTEGER TVAL, PRTSW8,NINC, IDEVNO , IDEVNT 
DETERMINE THE SUMS FOR EACH GROUP LEVEL 
DO 25 J=1 , 16 
GL (J) =0 . 

DO 30 1=1, LIN 

FRACT = CC (I, 4 ) /CDI (I, 9) 

TSTR = CDI (I , 8 ) 

NHB = IFIX(TSTR) 

GL(NHB) = GL {NHB ) + FRACT 
CONTINUE 

CALCULATE THE TLEVEL OF THE ASSOCIATED GROUP LEVELS 
TLEVL = GL ( 1 ) +GL (2 ) +GL ( 3 ) +GL ( 4 ) +GL ( 5) +GL (9) +GL { 10 ) +GL ( 11 ) + 

+ GL ( 13 ) +GL { 14 ) +GL ( 16 ) 

IF ( (PRTSW8.EQ. 1) .OR. ( (PRTSW8.EQ.0) .AND. (NINC.EQ.-l) ) ) THEN 
WRITE (IDEVNO, *, IOSTAT=IOVAL, ERR=900) ' 1 

WRITE ( IDEVNO, * , IOSTAT=IOVAL, ERR=900) ' ' 

WRITE ( IDEVNO, * , IOSTAT=IOVAL , ERR=900) 

+' GROUP T-VALUES AS SPECIFIED IN NHB 8060. IB APPENDIX D 1 
WRITE (IDEVNO, *, IOSTAT=IOVAL, ERR=900) ' ' 

WRITE (IDEVNO, * , IOSTAT=IOVAL, ERR=900) 

+' -01- -02- -03- -04- -05- -06- -07- -08- -09- 

+ -10- -11- -12- -13- -14- -15- -16- 1 

WRITE ( IDEVNO, 1 ( /IX, 16 (F6 . 2 , IX) / ) ' , IOSTAT=IOVAL , ERR=900 ) GL ( 1 ) 

+ GL(2) ,GL(3) ,GL(4) ,GL(5) ,GL(6) ,GL(7) ,GL(8) ,GL(9) ,GL(10) , 

+ GL (11) , GL ( 12 ) , GL ( 13 ) ,GL(14) ,GL(15) ,GL(16) 

WRITE (IDEVNO, *, I OSTAT= IOVAL, ERR=900) 1 OVERALL T- VALUE * 

OALLT = GL ( 1 ) +GL ( 2 ) +GL ( 3 ) +GL ( 4 ) +GL ( 5 ) +GL ( 9 ) +GL (10) +GL (11)+ 

+ GL (13) +GL ( 14 ) +GL ( 16 ) 

WRITE ( IDEVNO, ' (T4 , F7 . 2 ) 1 , IOSTAT=IOVAL, ERR=900 ) OALLT 
END IF 

IF (NINC.NE.-l) THEN 

************* WRITE t-VALUE DATA TO A PLOT FILE *************** 

IF (TVAL.EQ.2) THEN 

WRITE (IDEVNT, 50, IOSTAT=IOVAL,ERR=900) TNI, TN, TLEVL 
FORMAT (T2,2(F8.2,1X) ,F7.2) 

END IF 
END IF 
GOTO 999 

********************** 


900 

999 


WRITE ( * 
RETURN 
END 


■ ) ' IO ERROR IN 


END OF 
GROUP = 


SUBROUTINE 
' , IOVAL 


********************** 
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NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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RM/ FORTRAN Compiler (V2.42) 

Source File: C:\FORTRAN\TCC\HEADGS Options: /C 80 /L /BIJY 03/17/94 09:59:34 

1 C FILE :HEADGS. FOR 

2 Q ***★***★*****★★★******★**★★*★**★★+**★*★*****★★*★★*★★*★******★*** 

3 C * SUBROUTINE HDG1 * 

4 C * PROGRAM TO PRINT HEADING- DATE , TIME, FILE NAME, & PAGE NO. * 

5 C ★******★★★★*★★★★★************★*★★★★***★★*★★**********★*★***★**** 

6 

7 C NOTE: FILE MUST BE OPEN BEFORE STARTING THIS SUBROUTINE 

8 

9 SUBROUTINE HDG1 ( I MONTH , IDAY, IYEAR, IHOUR, IMINUTE, 

10 +FNAME, IPGNO, IDEVNO) 

11 C IMONTH. . . ISECOND=TIME AND DATE NAMES 

12 C FNAME=FILE NAME 

13 C I PGNO=PAGE NUMBER 

14 C I DEVNO= DEVICE NUMBER FOR OUTPUT 

15 

16 CHARACTER FNAME*24 

17 

18 

19 C WRITE HEADING 

20 IF (IPGNO .EQ. 1) THEN 

21 WRITE (IDEVNO, 5, IOSTAT=IOVAL,ERR=900) 

22 005 FORMAT (2X, 1 PROGRAM VERSION 8.1 Alpha ', 5X ,' March 15, 1994’,/) 

23 ENDIF 

24 WRITE ( IDEVNO ,10, IOSTAT=IOVAL , ERR=9 00 ) IMONTH , IDAY, IYEAR, 

25 + IHOUR, IMINUTE, FNAME, IPGNO 

26 010 FORMAT(2X, 12, 1 /’ , 12, ' / ' , I4,5X, 12, ' : ' ,12, 5X,A,2X, 'PAGE ',14) 

27 

28 GO TO 999 

29 900 WRITE (*,*)' IO ERROR IN HDG1= ' , IOVAL 

30 999 RETURN 

31 END 

NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 

32 

33 C **************************************************************** 

34 C * SUBROUTINE DATTM * 

35 C * SUBROUTINE TO READ SYSTEM DATE AND TIME FOR IBM PC OR AT * 

36 C **************************************************************** 

37 

38 SUBROUTINE DATTM ( IMONTH , IDAY, IYEAR, IHOUR, IMINUTE, ISECOND) 

39 

40 C REQUIRED FOR IBM PROF FORTRAN 

41 INTEGER *2 IMONTH , IDAY, IYEAR, IHOUR, IMINUTE, ISECOND, IHUNSEC 

42 

43 CALL GETDATf IYEAR, IMONTH, IDAY) 

44 CALL GETTIM( IHOUR, IMINUTE, ISECOND, IHUNSEC) 

45 

4 6 RETURN 

47 END 

NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
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48 

49 C 

50 C 

51 C 

52 C 

53 C 

54 

55 C 

56 

57 

58 

59 C 

60 C 

61 C 

62 C 

63 C 

64 

65 C 

66 

67 

68 

69 

70 

71 

72 

73 


* * 
* 

* 

★ 

★ * 


*************************************** 


*★*******★★*★★**★***★★* 


SUBROUTINE HDG2 

PROGRAM TO PRINT HEADING-TIME INCR+INCR INIT AND FINAL TIME * 
PRINTS TIME INCREMENT NUMBER 

************************************************************** 


NOTE: FILE MUST BE OPEN BEFORE STARTING THIS SUBROUTINE 
SUBROUTINE HDG2 ( INCRNO , TNI , TN , IDEVNO ) 

INPUTS : 

INCRNO=TIME INCREMENT NUMBER 
TNI = INCREMENT INITIAL TIME (HRS) 

TN= INCREMENT FINAL TIME (HRS) 

IDEVNO=DEVICE NUMBER FOR OUTPUT 

WRITE HEADING 

WRITE ( IDEVNO, 10, IOSTAT=IOVAL , ERR=900 ) INCRNO, TNI , TN 
010 FORMAT ( IX, ' TIME INCR 1 , 15 , 2X ,' INITIAL TIME (HRS)= ',F8.2,2X, 
+• FINAL TIME (HRS)= ',F8.2) 

GO TO 999 

900 WRITE (*,*)' IO ERROR IN HDG2= 1 , IOVAL 
999 RETURN 
END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 


74 

75 C 

76 C 

77 C 

78 C 

79 C 

80 

81 C 

82 

83 

84 

85 C 

86 C 

87 C 

88 C 

89 C 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 
101 
102 
103 


**************************************************************** 

* SUBROUTINE HDG3 

* PROGRAM TO PRINT HEADING-TIME INCR+INCR INIT AND FINAL TIME * 

* PRINTS PCALC OR FINAL INSTEAD OF TIME INCREMENT NUMBER 
**************************************************************** 


NOTE: FILE MUST BE OPEN BEFORE STARTING THIS SUBROUTINE 
SUBROUTINE HDG3 ( IFLAG, TNI , TN, IDEVNO) 

INPUTS : 

IFLAG=FLAG FOR TIME INCREMENT (1=PCALC, 2=FINAL ) 

TNI = INCREMENT INITIAL TIME (HRS) 

TN= INCREMENT FINAL TIME (HRS) 

I DE VNO= DEV ICE NUMBER FOR OUTPUT 

CHARACTER INAME* 5 

IF (IFLAG. EQ. 1) THEN 
INAME = ' PCALC ' 

ELSEIF ( IFLAG . EQ. 2 ) THEN 
I NAME = ' FINAL 1 
ELSE 

INAME =' ERROR' 

ENDIF 

C WRITE HEADING 

WRITE (IDEVNO, 10, IOSTAT= IOVAL, ERR=900 ) INAME, TNl,TN 
010 FORMAT (IX, 'TIME INCR ', A, 2X ,' INITIAL TIME (HRS)= ',F8.2,2X, 
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104 +• FINAL TIME (HRS ) = ' , F8 . 2 ) 

105 

106 GO TO 999 

107 900 WRITE (*,*)* 10 ERROR IN HDG3 = ' , IOVAL 

108 999 RETURN 

109 END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 


110 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 
12 3 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 


C 

C 

C 

C 

C 


*********************************** + * i ,i, i ,. ki , i ,i ci ,i ri ' i ' i , i ' i ' iriri , i ' i ' i ' 1 ' iriri ' iri ' ir 


* SUBROUTINE HDG4 

* PROGRAM TO PRINT HEADING-CONT NO., 

* MAC, EXCEEDS MAC 
************************************ 


NAME, FINAL CABIN CONC 
************************ 


★ 

★ 


★ 

★ ★ ★ ★ 


C NOTE: FILE MUST BE OPEN BEFORE STARTING THIS SUBROUTINE 
SUBROUTINE HDG4 ( IDEVNO) 


C INPUTS : 

C IDEVNO=DEVICE NUMBER FOR OUTPUT 
C WRITE HEADING 

WRITE (IDEVNO, 10 , IOSTAT=IOVAL, ERR=900) 

010 FORMAT (IX, 1 CONT ' , 14X, 'NAME 1 ,14X, 'FINAL CABIN' ,5X, 'MAC 1 5X 
+ ' EXCEEDS ' ) ' 

WRITE ( IDEVNO, 20 , IOSTAT=IOERR, ERR=900) 

020 FORMAT ( IX , 1 NO . 1 , 32X , 1 CONC (MG/M3 ) 1 , 12X, 1 MAC *) 

GO TO 999 

900 WRITE ( * , * ) 1 IO ERROR IN HDG4= 1 , IOVAL 
999 RETURN 
END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 


135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 


C 

C 

c 

c 

c 


* SUBROUTINE HDG5 * 

* PROGRAM TO PRINT HEADING-TOTAL CONT REMOVED BY EACH DEV (MG) * 

* PRINTS SHEET 1 -NO, NAME , CABIN , LEAK, &DEV3 .. DEV8 


NOTE: FILE MUST BE OPEN BEFORE STARTING THIS SUBROUTINE 
SUBROUTINE HDG5 (IDEVNO) 

C INPUTS : 

C IDEVNO=DEVICE NUMBER FOR OUTPUT 


C WRITE HEADING 

WRITE (IDEVNO, 10 , IOSTAT=IOVAL, ERR=900 ) 

010 FORMAT ( IX , 2 4X , 1 TOTAL CONTAMINANT MASS REMOVED BY EACH 
+ SHEET 1 1 ) 

WRITE (IDEVNO, 20, IOSTAT=IOVAL, ERR=900) 

020 FORMAT (IX, ' NO. ', 14X ,' NAME 1 , 1 6X ,* CABIN 1 , 9X, 1 LEAK 1 , 8X, 
+8X, ' DEV4 1 , 8X, ' DEV5 ' ,8X, ' DEV6 1 ,8X, ' DEV7 ' ,8X, ' DEV8 1 ) 


DEVICE (MG) 


DEV 3 ' , 
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156 

157 GO TO 999 

158 900 WRITE (*,*)' 10 ERROR IN HDG5= ', IOVAL 

159 999 RETURN 

160 END 

NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 


161 

162 C 

163 C 

164 C 

165 C 

166 C 

167 

168 C 

169 

170 

171 

172 C 

173 C 

174 

175 C 

176 

177 

178 

179 

180 
181 
182 

183 

184 

185 

186 


**************************************************************** 

* PROGRAM°TO PRINT HEADING-TOTAL CONT REMOVED BY EACH DEV {MG) * 

* PRINTS SHEET 2-NO, NAME, & DEV9 .. DEV15 *************** * 

***************************************************** 


NOTE: FILE MUST BE OPEN BEFORE STARTING THIS SUBROUTINE 
SUBROUTINE HDG6 ( IDEVNO) 


INPUTS: 

IDEVNO=DEVICE NUMBER FOR OUTPUT 


WRITE HEADING 

ud T TP i T DFVNO 10 , IOSTAT=IOVAL , ERR=900 ) 

010 FORMAT ( IX, 2 4X, 'TOTAL CONTAMINANT MASS REMOVED BY EACH DEVICE (MG) 

SHEET 2 ' ) 


WRITE (IDEVNO, 20, IOSTAT=IOERR,ERR=900) iTW , M _ i _ v 

020 FORMAT ( IX , ' NO. ' , 14X, ’ NAME ' , 16X, ' DEV9 ' , 8X, DEV10 , 7X , 
+7X, ' DEV12 ' , 7X , ' DEVI 3 ' , 7X, ' DEV14 ' , 7X, ' DEVI 5 ' ) 


DEVI 1 1 


GO TO 999 

900 WRITE (*,* ) 'IO ERROR IN HDG6= 1 , IOVAL 
999 RETURN 
END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 


187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 
207 


**************************************************************** 

* SUBROUTINE HDG7 * 

* PROGRAM TO PRINT HEADING-DEVICE REM EFF AT END OF TIME INCR * 



NOTE: FILE MUST BE OPEN BEFORE STARTING THIS SUBROUTINE 
SUBROUTINE HDG7 (IDEVNO) 


C INPUTS : 

C IDEVNO=DEVICE NUMBER FOR OUTPUT 


C WRITE HEADING 

WRITE ( IDEVNO, 10 , IOSTAT=IOVAL, ERR=900) 

010 FORMAT ( IX, 24 X, ' DEVICE REMOVAL EFFICIENCY 


AT END OF 


TIME INCREMENT 


+ (DEC) ' ) 

WRITE { IDEVNO, 20 , IOSTAT=IOERR, ERR=900) 

020 FORMAT ( IX, ' NO. ' , 14X, ' NAME ' , 15X, ' #2 ' , 4X, ' #3 ' , 
+ 4X, '#4' ,4X, '#5' ,4X, '#6' ,4X, '#7' ,4X, '#8' ,4X, '#9 


3X, '#10 
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208 

209 

+ 3X, ' #11 ' , 3X, ' #12 ' , 3X, 

1 #13 ' , 3X, ' #14 ' , 3X, ' #15 ' ) 

210 

GO TO 999 


211 

900 WRITE (*,*)' IO ERROR 

IN HDG7 = 1 , IOVAL 

212 

999 RETURN 


213 

END 


NUMBER 

OF WARNINGS IN PROGRAM 

UNIT: 0 

NUMBER 

OF ERRORS IN PROGRAM UNIT: 0 


214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 


C 

C 

C 

C 

C 

C 


C 

C 


* SUBROUTINE HDG8 * 

* PROGRAM TO PRINT HEADING-RATE OF CONT REMOVAL- EACH DEV (MG) * 

PRINTS SHEET 1 -NO, NAME, CABIN, LEAK, &DEV3 . .DEV8 * 


NOTE: FILE MUST BE OPEN BEFORE STARTING THIS SUBROUTINE 
SUBROUTINE HDG8 (IDEVNO) 

INPUTS: 

IDEVNO=DEVICE NUMBER FOR OUTPUT 


-EACH DEVICE (MG/HR) 


: WRITE HEADING 

WRITE (IDEVNO, 10, IOSTAT= IOVAL , ERR=900) 

010 FORMAT ( IX, 2 4X ,' RATE OF CONTAMINANT REMOVAL- 
+ SHEET 1 • ) 

WRITE ( IDEVNO, 20 , IOSTAT= IOVAL , ERR=900) 

020 FORMAT (IX, ' NO. 1 , 14X, 'NAME' ,16X, 'CABIN' ,9X, 'LEAK' ,8X, ' DEV 3 
+8X, ' DEV 4 , 8X , DEV 5 ' , 8X , ‘ DEV 6 ' , 8X , ' DEV7 ' , 8X, ' DEV8 ' ) 


GO TO 999 

900 WRITE ( * , * ) ' IO ERROR IN HDG8= ', IOVAL 
999 RETURN 
END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 


240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 


C 

C 

C 

C 

C 

C 


* SUBROUTINE HDG9 

PROGRAM TO PRINT HEADING-RATE OF CONT REMOVAL-EACH DEV (MG/HR 
PRINTS SHEET 2-NO, NAME, & DEV9..DEV15 


* 

★ 

★ 

* 

* 


NOTE: FILE MUST BE OPEN BEFORE STARTING THIS SUBROUTINE 
SUBROUTINE HDG9 (IDEVNO) 


C INPUTS : 

C IDEVNO=DEVICE NUMBER FOR OUTPUT 
C WRITE HEADING 

WRITE (IDEVNO, 10, IOSTAT=IOVAL , ERR=900) 

010 FORMAT ( IX, 2 4X, ' RATE OF CONTAMINANT REMOVAL-EACH DEVICE (MG/HR) 
+ SHEET 2 1 ) 

WRITE (IDEVNO, 20 , IOSTAT= IOERR , ERR=900 ) 

020 FORMAT ( IX , ' NO . ' , 14X , ' NAME ' , 1 6X , ' DEV9 • , 8X , • DEV10 ' , 7X , 1 DEVI 1 ' , 
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260 +7X , ' DEVI 2 1 , 7X , 1 DEVI 3 ' , 7X, 'DEVl4',7X, 1 DEV15 1 ) 

261 

'y c o GO TO 9 99 

263 900 WRITE (*,*)' 10 ERROR IN HDG9= ' , IOVAL 

264 999 RETURN 

265 END 

266 

NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 


* SUBROUTINE LDIGEN 


. SUBROUTINE TO LOAD INTERNAL GENERATION FROM MAT CDI COL 1 L ; 



SUBROUTINE LDIGEN ( I , DD , NROW , NCOL , CDI , NROW2 , NCOL2 , LIN2 ) 

INTEGER NROW, NCOL , NROW2 , NCOL2 , LIN2 
REAL DD (NROW, NCOL ) 

REAL CDI (NROW2 , NCOL2 ) 


RM/FORTRAN Compiler (V2.42) 05/21/92 12:56:50 

Source File: C: \™FORT\TCC\LDIGEN. Options: 

1 C 

2 C 

3 C 

4 C 

5 C 

6 

7 

8 
9 

10 

11 C INPUTS: 

12 C I =C ONTAM I NANT NUMBER 

13 C DD, NROW, NCOL =NAME AND DIMENSIONS OF MAT DD 

14 C CDI , NROW2 , NCOL2 =NAME AND DIMENSIONS OF MAT CDI 

15 C LIN2=NUMBER OF DEVICES IN MAT DD 

16 C OUTPUT^ internal GENERATI on FROM MAT CDI INTO MAT DD COL 19 

18 C 

19 DD ( 1 , 19 ) =CDI (1,1) 

20 DD (2 , 19 ) =0 

21 DO 10 J=3 , LIN2 

22 DD (J, 19) =CDI (I,J+7) 

23 010 CONTINUE 

25 q ^**** ******** END OF SUBROUTINE LDIGEN ********************* 

26 END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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RM/ FORTRAN Compiler (V2.42, 

Source File: C:\RMFORT\TCC\LODEFF. Options: /C 80 /L /BY 


1 
2 

3 

4 

5 

6 

7 

8 
9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 C 

25 


05/21/92 12:57:05 


*.**** 0 *** 0 **^,**********,,****,,.,***.,,*,**,,*,;;*^;;^;^, 

* SUBROUTINE LODEFF * 

SUBROUTINE TO LOAD LAST INCR EFF FROM MAT CC INTO MAT DD COL 20* 

* USES ADJUSTABLE SIZE ARRAYS * 

SUBROUTINE LODEFF ( I , DD, NROW, NCOL , CC , NR0W1 , NC0L1 LIN2 ) 

INTEGER NROW, NCOL, NROWl , NCOL 1 

REAL DD (NROW, NCOL) 

REAL CC (NROWl , NC0L1 ) 

INPUTS : 

I =CONTAM I NANT LINE NUMBER IN MAT CC 
DD, NROW, NCOL=NAME & DIMENSIONS OF MAT DD 
CC, NROWl ,NC0L1=NAME &. DIMENSIONS OF MAT CC 
LIN2=N0. OF DEVICES IN MAT DD 


DD(1, 20) =0 
K=7 

DO 100 J=2 , LIN2 
DD ( J , 2 0 ) =CC ( I , K) 

K=K+3 
100 CONTINUE 
RETURN 

************* END OF SUBROUTINE LODEFF ************************* 
END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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RM/FORTRAN Compiler (V2.42) 
source File: C:\RMFORT\TCC\MASBAL. 


05/21/92 12:57:51 


1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 


: *Sii^^?ii^*.““i*«*i**“ i ****** ; 

* MASS BALANCE SUBROUTINE -MASBAL r&T PUT ATES * 

* rnR 1 PONT AT A TIME AT A GIVEN DEVICE EFFICIENCY CALCULAT S 

! . 5 ?.. • 

NOTE: ( Zh 0 R b E v R S, WFSFS? S T «ch RO d“ice DD in?o 17 ' 21 

MAT DD COL 20 (DONE BY LODEFF OR CNRSUB) 

SUBROUTINE MASBAL ( I , TN , TNI , DD , NROW , NCOL , CC , NROWl , NCOL 1 , 

+CAVCLC, CDI ,NROW2 ,NCOL2 , CFINAL, CEQLIB, LIN, LIN2) 

INTEGER NROW , NCOL , NROWl , NCOLl , NROW2 , NCOL2 
REAL DD (NROW, NCOL) 

REAL CC (NROWl , NCOLl ) 

REAL CDI (NROW2 , NCOL2 ) 

saw a - - - » 

PCAVCF-USING CEQ & CINIT, CALC CFINAL & CAVERAGE 


C INPUTS: 

C FROM PCSET, PREDCT, AND CONVRG 

C TN^INCREMENT T END ’ TIME ( HRS ) ; TNI=INCR BEGINNING TIME HRS 

C DD, NROW, NCOL =NAME AND SIZE OF MAT DD 

C CC , NROWl , NCOLl =NAME AND SIZE OF MAT CC 

C CDI , NROW2 , NCOL2=NAME AND SIZE OF MAT CDI 
C LIN=NO. OF CONT IN MAT CDI 
C LIN2 =NO. OF DEVICES IN MAT DD 

£ FR fMRFM L SUM OF MASS REM FOR ALL DEVICES (MG/HR ) -TOTAL OF DD COL 21 
? iSlt Sm OF £11 G^N IN ALL DEVICES INCL CABIN (MG/HR) -DD COL 19 

c FR "t LOADS^CABIN M.GEN (MG/HR) FROM MAT CDI INTODDU^l^ 

C IT LOADS M.GEN DEVICES FROM MAT CDI COL 2-15,19 INTO DD COL 9 

C FROM PCAVCF 

C CAVCLC =C ALC INCR CABIN CONT CONC (MG/CU M) 

C CFINAL=FINAL INCR CABIN CONT CONC (MG/CU M) 

C OUTPUTS: 

C TO PCSET, PREDCT, AND CONVRG 

C CAVCLC CALCULATED CABIN AVERAGE CONC (MG/CU M) 

C CEOLIB=CABIN EQUILIBRIUM CONCENTRATION (MG/CU M) 

C C F INAL = INCREMENT FINAL CABIN CONCENTRATION (MG/CU M) 

C PUTS M.REM FOR CABIN + DEVICES IN MAT DD COL 21 
C TO CALCM 

C CAV=CABIN CONT CONC (MG/CU M) 

C OTliE^DEV^E^DD^-lS 0 ^) MUST =0 AT THIS POINT (SEE PRAFIL) 

C TO LDIGEN 

C I = CONTAMINANT NUMBER 

C DD, NROW, NCOL =NAME AND DIMENSIONS OF MAT DD 
C CDI NROW2 , NCOL2 =NAME AND DIMENSIONS OF MAT CDI 
C LIN2=NUMBER OF DEVICES IN MAT DD 
C TO PCAVCF 
P TN TNI 

57 C CINIT= INITIAL INCR CABIN CONT CONC (MG/CU M)=CC(I,1) 

58 C CEQLIB=CABIN EQULIB CONC (MG/CU M) 
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59 C SQEFFN=SUM OF Q*REM EFF NET FOR ALL DEVICES (CU M/HR) 

60 C CVOL=CABIN VOL (CU M)=DD(1,9) 

61 C SMNTC=SUM OF MASS NET TO CABIN (MG/HR) 


63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 


C CABIN VOL (CU M) 
CVOL=DD ( 1 , 9) 

C C INITIAL (MG/CU M) 
CINIT=CC ( I , 1) 


C 

C 

C 

C 

C 

C 

C 


C 


C 


C 


EVALUATE SUM Q*REM EFF NET USING M.GEN IN DEVICES=0 (DD COL 19) 

SET CABIN AVERAGE CONCENTRATION = TO ARBITRARY VALUE OF 100 
AND INTERNAL GENERATION IN DEVICES =0 (NOT YET LOADED) 

SET CABIN M . GEN=ARBITRARY VALUE OF 50 (DD(1,19)) 

SMGEN=SUM M.GEN IN ALL DEVICES +CABIN (MG/HR) 

SMREM=SUM M. REMOVED BY ALL DEVICES (MG/HR) 

CAV=CABIN AVERAGE CONCENTRATION (MG/CU M) 

CAV=100 
DD (1, 19) =50 

CALL CALCM ( DD, NROW, NCOL , CAV , SMGEN , SMREM , TN , LIN2 ) 

SQEFFN=SUM OF Q* REMOVAL EFF NET (MG/HR) 

SQEFFN= SMREM /CAV 

LOAD INTERNAL GENERATION FOR ALL DEV+CABIN FROM CDI INTO DD COL 19 
CALL LDIGEN ( I , DD, NROW, NCOL, CDI , NROW2 , NCOL2 , LIN2 ) 

EVALUATE SUM OF M . NET TO CABIN=M . GEN CABIN+SUM M.GEN ALL DEVICES - 
SUM M.REM ALL DEVICES 

SMNTC=SUM M . NET TO CABIN=AMT GEN WHICH GETS TO THE CABIN DIRECTLY 
SET C CABIN AV= 0 """ ' 

CAV=0 

GET SUM MASS GEN CABIN+ INTERNAL DEVICES AND SUM MASS REMOVED ALL 
DEVICES FROM SUBROUTINE-SINCE CABIN C = 0 NO CABIN CONT WILL BE REM 
CALL CALCM ( DD, NROW, NCOL , CAV, SMGEN, SMREM, TN, LIN2 ) 

SMNTC=SMGEN- SMREM 

NOTE : SMNTC IS ALSO PUT IN DD (1,21) BY CALCM 
GET CALCULATED CABIN EQUILIBRIUM CONCENTRATION (CAVCLC) (MG/CU M) 

IF ( SQEFFN . LT . IE- 1 0 ) THEN 
IF(CVOL.EQ.O) THEN 
CFINAL=1E10 
ELSE 

CF INAL =C INI T+ (TN-TN1) *SMGEN/CVOL 
END IF 

CAVCLC = ( C INI T+CF INAL ) /2 
CEQLIB=1E10 
ELSE 

CEQLIB=SMNTC /SQEFFN 

CALCULATE CAVCLC AND CFINAL FROM SUBROUTINE 
CALL PCAVCF (TN, TNI ,CEQLIB, SQEFFN, CVOL, SMNTC, CINIT 
+ CAVCLC, CFINAL) 

END IF 

USING CAV CALC EVALUATE M.REM FOR CABIN + DEVICES AND PUT 
IN MAT DD COL 21 
CAV=CAVCLC 

CALL CALCM ( DD , NROW, NCOL , CAV , SMGEN , SMREM , TN, LIN2 ) 
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119 

120 C 

121 C 

122 


RETURN 

END OF SUBROUTINE MASBAL 

************** END OF SUBROUTINE MASBAL 
END 


************************** 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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RM/FORTRAN Compiler (V2.42) 

Source File: C:\FORTRAN\TCC\MCALC. Options: /C 80 /L /BIJY 03/15/94 10:40 

***************************** + ***i,*i,i,** i , i , i , i ' i , i , iri , ir i ti , i , i , iti ' i ' i ' + ******** 


1 

2 

3 

4 

5 

6 
7 


C 

C 

C 

C 

C 

C 

C 


SUBROUTINE MCALC 

MAIN CALCULATION LOOP SUBROUTINE FOR 1 TIME INCREMENT 
FOR ALL CONTAMINANTS ONE AT A TIME 
BASED ON SUM MASS REM LAST INCR, FOR EACH CONT 
CALCULATE NEW REMOVAL EFF, CAV CALC CABIN, 

CEQUILIB, CFINAL , & M. REMOVED ALL DEVICES- PUT IN MAT CC 


★ 

★ 

★ 

★ 

* 

* 


8 Q *********************************************************^******* 

9 SUBROUTINE MCALC ( I , TN, TNI , DD , NROW, NCOL , 

+CC , NROW1 , NCOL1 , CDI , NROW2 , NCOL2 , NN , LIN , LIN2 , 

+PRTSW2 , PRTSW3 , PRTSW4 , IMSGDN) 

INTEGER NROW , NCOL , NROW1 , NCOLl , NROW2 , NCOL2 , PRTSW2 , PRTSW3 PRTSW4 KK 
CHARACTER NN(NROW1)*30 ' ' 

REAL DD (NROW, NCOL) 

REAL CC (NROW1, NCOLl) 

REAL CDI (NROW2 , NCOL2 ) 

SUBROUTINES REQUIRED: 

PREDCT= PREDICT CAV BASED ON M. GEN OF THIS INCR & REM EFF OF LST INC 
CONVRG=CALC CAV CALC , CEQ, CFINAL, M. REM, REM EFF 
CROUT=PRINT TEST VALUES OF MAT CC 
RROUT= PRINT TEST VALUES OF MAT DD 

INPUTS : 

FROM MAIN PROG 
I=CONTAMINANT NO. 

TN, TNI = INCREMENT END & BEGINNING TIME (HRS) 

DD,NROW, NCOL=NAME & DIM OF MAT DD 
CC,NROWl , NCOLl =NAME & DIM OF MAT CC 
CDI , NROW2 , NCOL2 =NAME & DIM OF MAT CDI 
NN=NAME OF MAT NN 

LIN=NUMBER OF CONTAMINANTS IN MAT NN & CDI 
LIN2=NO. DEVICES IN MAT DD 
FROM PREDCT 

CAVPRD=PRED CABIN CONT CONC (MG/CU M) 

FROM CONVRG 

CAVCLC=CALC CABIN CONT CONC (MG/CU M) 

CFINAL=FINAL INCR CABIN CONT CONC (MG/CUM) 

KK=COUNTER FOR CONVERGENCE 

CEQLIB=EQUILIBRIUM CABIN CONT CONC (MG/CU M) 

IMSGDN=DEVICE NO FOR MESSAGE AND TEXT PRINTOUT OUTPUT 
OUTPUTS : 

TO PREDCT 
I =CONTAMINANT NO. 

TN , TNI = INCREMENT END & BEGINNING TIME (HRS) 

DD, NROW, NCOL=NAME & DIM OF MAT DD 
CC,NROWl , NCOLl =NAME & DIM OF MAT CC 
CDI , NROW2 , NCOL2=NAME & DIM OF MAT CDI 
LIN=NO . OF CONTAMINANTS IN MAT CDI 
LIN2 =NO . DEVICES IN MAT DD 
TO CONVRG 

C I = C ONTAM I N ANT NO. 

C TN , TNI INCREMENT END & BEGINNING TIME (HRS) 

C DD, NROW, NCOL =NAME & DIM OF MAT DD 

54 C CC , NROW1 , NCOLl=NAME & DIM OF MAT CC 

55 C CDI , NROW2 , NCOL2=NAME & DIM OF MAT CDI 

56 C LIN=NO . OF CONTAMINANTS IN MAT CDI 

57 C LIN2 =NO . DEVICES IN MAT DD 


10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 


36 C 

37 C 

38 C 

39 C 


40 

41 

42 

43 

44 

45 

46 

47 

48 C 

49 C 

50 C 

51 

52 

53 


: 19 
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c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 C 

74 C 

75 

76 

77 

78 C 

79 C 

80 
81 

82 C 

83 C 

84 

85 C 

86 

87 C 

88 C 

89 

90 

91 C 

92 

93 C 

94 

95 C 

96 C 

97 

98 

99 

100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 
117 


CAVPRD= PREDICTED CABIN CONC FOR INCREMENT (MG/CU M) 

TO MAT CC „ 

PUTS CAVCLC , CEQLIV, &CFINAL IN CC(I,2-3 &4) 

PUTS REM EFF FROM DD COL20 IN CC(I, 7-10-13 ETC) 

PUTS M. REM FOR EACH DEV FROM DD COL21 IN CC (I, 6-9-12 
PUTS SUM MASS REM FOR EACH DEV IN CC(I, 8-11-14 ETC) 


ETC) 


BEGIN LOOP FOR EACH CONTAMINANT FOR EACH TIME INCREMENT 

D CALC°CAV 1 PRED CABIN FOR CONT BASED ON REM EFF OF LAST INCREMENT 
AND GENERATION RATES OF THIS INCREMENT 
CALL PREDCT ( I , TN , TNI , CAVPRD , DD , NROW, NCOL , 

+CC , NROWl , NCOLl , CDI ,NROW2 ,NCOL2 , CAVCLC , CFINAL , CEQLIB , LIN, LIN2 ,NN) 

CONVERGE UNTIL CCALC=CPRED 
CALL CONVRG ( I , TN , TNI , CAVPRD , DD , NROW , NCOL , 

+CC , NROWl , NCOLl , CDI , NROW2 , NCOL2 , CAVCLC , CFINAL , CEQLIB , KK , LIN , 

+LIN2 ,NN, PRTSW3 , IMSGDN) 

IF KK>2 0 THEN BEGIN 1/20 TIME INCREMENT CONVERGENCE ROUTINE 
IF (KK.GT.20) THEN 
KK=1 

BEGIN 1/20 INCREMENT CONVERGENCE ROUTINE 
NEW INCREMENT INITIAL TIME (HRS) 

TN1NEW=TN1 

NEW TIME INCREMENT (HRS) 

BINEW= (TN-TN1 ) /20 

BEGIN LOOP FOR 1/20 INCREMENT SIZE TIME INCREMENT 
NEW INCREMENT FINAL TIME (HRS) 

200 TNNEW=TNlNEW+BINEW 

! ZERO MAT DD COL 17-21 

CALL PRAFIL ( DD , NROW , NCOL ,17,21) 

; LOAD EFFICIENCY FROM LAST INCREMENT INTO MAT DD COL 20 
CALL LODEFF ( I, DD, NROW, NCOL, CC , NROWl, NCOLl, LIN2 ) 

; CALC CAV PRED CABIN FOR CONT BASED ON REM EFF OF LAST INCREMENT 
; AND GENERATION RATES OF THIS INCREMENT 

CALL PREDCT ( I , TNNEW , TN1NEW , CAVPRD , DD , NROW , NCOL , 

+ CC , NROWl , NCOLl , CDI , NROW2 , NCOL2 , CAVCLC , CFINAL , CEQLIB , LIN, 

+ LIN2 , NN) 

C 

C CONVERGE UNTIL CCALC=CPRED 

CALL CONVRG ( I , TNNEW , TN1NEW, CAVPRD , DD , NROW , NCOL , 

+ CC , NROWl , NCOLl , CDI , NROW2 , NCOL2 , CAVCLC , CFINAL , CEQLIB , KK , LIN , 

+ LIN2,NN, PRTSW3, IMSGDN) 

IF KK>20 THEN PRINT CONVERGENCE WARNING 
IF (KK.GT.20) THEN 

OPEN (IMSGDN, FILE=' CON' , IOSTAT=IOVAL ) 

WRITE ( IMSGDN ,* ) 'WARNING: CALCULATION DID NOT CONVERGE FOR 
WRITE (IMSGDN, *) ' FULL AND 1/20 INCREMENT ROUTINES' 

WRITE (IMSGDN, 50) I , NN ( I ) , TN1NEW, TNNEW 
050 FORMAT ( IX, ' CONT NO. = ' , 14 , 2X, A, / , IX , 

+ ‘FOR INCREMENT INIT & FINAL TIMES= ’,F8.2,F8.2) 

CLOSE (IMSGDN) 

ENDIF 

: FILL MAT CC WITH RESULTS 


C 

c 
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118 C PUT CAVCLC , CEQLIB , AND CFINAL IN CC 

119 CC ( I , 2 ) =CAVCLC 

120 CC (1,3) ^CEQLIB 

121 CC (I, 4 ) =CFINAL 

122 C PUT REM EFF FROM LAST ITER DD COL 20 IN CC ( I , 7- 1 0 - 1 3 ETC ) 

123 K=7 

124 DO 102 J=2 , LIN2 

125 CC ( I , K ) =DD ( J , 2 0 ) 

126 K=K+3 

127 102 CONTINUE 

128 

129 C TAKE CABIN M. REM (MG /HR) FROM DD ( 1 , 2 1 ) & PUT IN MAT CC ( I , 5 ) 

130 CC(I, 5) =DD(1,21) 

131 C TAKE M.REM FROM DD COL 21 & PUT IN CC ( I , 8- 1 1 - 14 ETC ) 

132 K=8 

133 DO 103 J=2 , LIN2 

134 CC ( I , K ) =DD ( J, 2 1 ) 

135 K=K+3 

136 103 CONTINUE 

137 

138 C CALCULATE SUM OF MASS REMOVED IN DEVICES + CABIN TO DATE AND 

139 C PUT IN CC<I, 6-9-12ETC) 

140 K=5 

141 DO 104 J=1,LIN2 

142 CC ( I, K+l) =CC ( I, K+l ) +CC (I, K) * (TNNEW-TN1NEW) 

143 K=K+3 

144 104 CONTINUE 

145 C 

14 6 C IF PRTSW2 = 1 THEN PRINT MAT CC INFO FOR THIS CONTAMINANT 

147 IF ( PRTSW2 . EQ . 1 ) THEN 

148 OPEN ( IMSGDN, FILE= ' CON ' , IOSTAT=IOVAL) 

149 WRITE ( IMSGDN, *) 'NINC,TN, TNI ' , NINC , TN, TNI 

150 WRITE ( IMSGDN, *)' PRINTOUT FOR ONE CONT INSIDE 1/20 INCR 

151 + LOOP OF MCALC ' 

152 WRITE ( IMSGDN, * ) ' INFO FROM MAT CC ’ 

153 CLOSE (IMSGDN) 

154 CALL CROUT (NN, CC,NROWl ,NCOLl , 1 ,NCOLl ,LIN, I , I, IMSGDN) 

155 END IF 

156 C 

157 

158 C REPEAT LOOP FOR 1/10 INCREMENT IF END OF 1/20 INCREMENT TIME 

159 C (TNNEW) IS < THAN END OF LARGER TIME INCR (TN) 

160 C ELSE IF TNNEW>=TN , END 1/20 ENCR CONV & PRINT ANSWERS+REPEAT 

161 C FOR ANOTHER CONTAMINANT 

162 

163 IF (TNNEW. LT.TN) THEN 

164 C RESET FOR ANOTHER 1/20 TIME INCREMENT 

165 TN1NEW= TNNEW 

166 CC (1,1) =CC (1,4) 

167 GO TO 200 

168 ELSE 

169 C END 1/20 INCR CONV ROUTINE-REPEAT FOR ANOTHER CONT 

170 GOTO 100 

171 END I F 

172 

173 C END OF CONVERGENCE ROUTINE 

174 END IF 

175 

176 C CALC SUM MASS REMOVED & FILL MAT CC WITH RESULTS 

177 C PUT CAVCLC, CEQLIB, AND CFINAL IN CC 
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178 

179 

180 
181 
182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 


CC ( I , 2 ) =CAVCLC 
CC (1,3) =CEQLIB 
CC (1,4) =CFINAL 

C PUT REM EFF FROM LAST ITER DD COL 20 IN CC { 1 , 7 - 1 0- 13 ETC ) 

K=7 

DO 302 J=2 , LIN2 
CC ( I , K) =DD(J,20) 

K=K+3 

302 CONTINUE 

C TAKE CABIN M . REM (MG/HR) FROM DD(1,21) & PUT IN MAT CC(I,5) 

CC(I, 5) =DD(1,21) 

C TAKE M. REM FROM DD COL 21 & PUT IN CC ( I , 8-11-14ETC ) 

K=8 

DO 303 J=2 , LIN2 
CC { I , K ) =DD ( J , 2 1 ) 

K=K+3 

303 CONTINUE 

C CALCULATE SUM OF MASS REMOVED IN DEVICES + CABIN TO DATE AND PUT 
C IN CC ( I, 6- 9 -12 ETC) 

K=5 

DO 304 J=1 , LIN2 

CC ( I , K+l ) =CC ( I , K+l ) +CC ( I , K) * (TN-TNl ) 

K=K+3 

304 CONTINUE 
C 

C IF PRTSW4 = 1 THEN PRINT MAT DD+MAT CC INFO FOR THIS CONTAMINANT 
IF ( PRTSW4 . EQ. 1 ) THEN 
OPEN { IMSGDN, FILE= ' CON* , IOSTAT=IOVAL) 

WRITE ( IMSGDN, *) 'PRINTOUT FOR ONE CONT AT END OF MCALC ' 

WRITE ( IMSGDN, * ) 1 INFO FROM MAT CC ' 

CLOSE (IMSGDN) 

CALL CROUT (NN, CC , NROW1 , NCOLl , l,NCOLl,LIN, 1,1, IMSGDN) 

OPEN ( IMSGDN, FILE= ' CON ' , IOSTAT=IOVAL) 

WRITE ( IMSGDN, *)' INFO FROM MAT DD' 

CLOSE (IMSGDN) 

CALL RROUT ( DD , NROW , NCOL , 1 , NCOL , LIN2 , IMSGDN ) 

ENDIF 


C 

C 

c 


END OF I LOOP FOR EACH CONTAMINANT 


100 CONTINUE 
RETURN 

c ****** END op SUBROUTINE MCALC ********************************** 
END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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RM/ FORTRAN Compiler (V2.42) 

Source File: C:\RMFORT\TCC\PCAVCF. Options: /C 80 /L /BY 05/21/92 12 -57 -51 

1 C **************************************************************** 

* SUBROUTINE PCAVCF 

* SUBROUTINE TO PREDICT INCREMENT CALCULATED AVERAGE EFF 

* (C AV CALC), FINAL EFF (CFINAL) & CABIN CONTAMINANT 

* CONCENTRATION 

*************************** ************************************** 


2 C 

3 C 

4 C 

5 C 

6 C 

7 


8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 


+ CAVCLC ' T S I CFINAL) VCF (TN ' TN1 ' SCE Q LIB ' SQEFFN, CVOL, SMNTC, CINIT, 

DOUBLE PRECISION EXPON, CEQLIB , CFINAL 
CEQLIB=DBLE (SCEQLIB) 

SUBROUTINES REQUIRED: NONE 


INPUTS : 

TN, TN1=INITIAL & FINAL INCREMENT TIME 
=EQUILIBRIUM CABIN CONC 
NET (MG/HR) 


(HRS) 

(MG /CUM) 


SQEFFN=SUM Q* REMOVAL EFF 
CVOL=CABIN VOLUME (CU M) 

SMNTC = SUM MASS CONT NET TO CABIN (MG/HR) 

CINIT= INITIAL INCREMENT CONT CONC (MG/CU 
OUTPUTS: 

SCAVCLC (CAVCLC) =CALC AVERAGE CABIN CONC (MG/CU 
SCFINAL (CFINAL) = FINAL INCREMENT CONC (MG/CU M) 


M) 


M) 


26 

IF ( CVOL . LE . 0 ) THEN 

27 

CAVCLC =CEQLIB 

28 

CFINAL=CEQLIB 

29 

GOTO 9 9 

30 

END IF 


31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 


C 


099 


CALCULATION FOR CFINAL 
EXPON= (TN-TN1 ) *SQEFFN/CVOL 
IF(ABS(EXPON) .GT. 50) THEN 
CAVCLC =CEQL IB 
CFINAL=CEQLIB 
GOTO 99 
END IF 

IF (ABS ( EXPON) . LT. IE-6 ) THEN 

CFINAL=CINIT+ SMNTC * (TN-TN1 ) /CVOL 
CAVCLC= (CINIT+CFINAL) /2 
CEQLIB=1E10 
GOTO 99 
ELSE 

CFINAL=CINIT+ ( SMNTC/SQEFFN-CINIT) 

END IF 

CALCULATION FOR C AVERAGE CALC 

IF ( (CINIT. EQ. CFINAL) .OR. (CFINAL. EQ. CEQLIB) ) THEN 
CAVCLC =C F I NAL 
GOTO 99 
END IF 

IF ( (CEQLIB -CINIT) / (CEQLIB-CFINAL ) .LT.1E-6) 

CAVCLC = (CINIT+CFINAL) / 2 
ELSE 

CAVCLC=CEQLIB- (CFINAL-CINIT) /LOG ( ( CEQLIB-CINIT ) / 
+ (CEQLIB -CFINAL)) 

END IF 
CONTINUE 

SCFINAL=REAL (CFINAL) 


(1 -EXP (-EXPON) ) 


THEN 
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59 

60 C 

61 

62 C 


RETURN 

*********** 

END 


END OF SUBROUTINE PCAVCF ************************** 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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RM/ FORTRAN Compiler (V2.42) 

Source File: C:\RMFORT\TCC\PCSET.F Options: /C 80 /L /BY 05/21/92 12:58:31 

1 C ***************************************************************** 

2 C * SUBROUTINE PCSET * 

3 C * SUBROUTINE FOR PRECALCULATION SETUP ROUTINE * 

4 c * FOR ALL CONTAMINANTS ONE AT A TIME * 

l C * CALL EFF SUBROUTINES FOR DEVICES; GET CAV CABIN PRED, * 

6 C * CEQUILIB , CFINAL, & M. REMOVED ALL DEVICES-PUT IN MAT CC * 

7 C ***************************************************************** 

8 SUBROUTINE PCSET ( TNI , LIN, DD , NROW, NCOL , CC , NROW1 , NCOL1 , 

9 +CDI , NROW2 , NCOL2 , LIN2 , NN, PRTSW1 , IMSGDN) 

10 INTEGER NROW, NCOL, NROW1 , NCOL 1 ,NROW2 ,NCOL2 , PRTSW1 , KK 

11 CHARACTER NN(NROWl) *30 

12 REAL DD (NROW, NCOL) 

13 REAL CC (NROW1 , NCOLl ) 

14 REAL CDI (NROW2 , NCOL2 ) 

15 KK=0 

16 C SUBROUTINES REQUIRED: 

17 C PRAFIL-ZERO MAT DD COL 17-21 

18 C CNRSUB-USING CAV=lE-20, FIND REMOVAL EFF AND PUT IN MAT DD COL 20 

19 C MASBAL-CALC CAV CALC PRED, CEQ, CFINAL, M.REM 

20 C CROUT-TEST PRINTOUT OF CONT INFO 


INPUTS : 

FROM MAIN CALC LOOP 

TN 1 = INCREMENT INITIAL TIME (HRS) 

LIN=NO . OF CONT IN MAT CC AND NN 
DD, NROW, NCOL=NAME & SIZE OF MAT DD 
CC,NROWl , NCOLl =NAME & SIZE OF MAT CC 
CDI , NROW2 , NCOL2=NAME & SIZE OF MAT CDI 
LIN2 =NO . DEVICES IN MAT DD 
NN=NAME OF MAT NN 

PRTSW1 =PRINTSWITCH WHICH CONTROLS TEST PRINTOUT 
IMSGDN=DEVICE NUMBER FOR MESSAGES AND TEST PRINTOUT 
FROM PRAFIL 

PUTS ZEROS IN MAT DD COL 17-21 
FROM CNRSUB 

CNRSUB PUTS REM EFF (DEC) FOR EACH DEVICE IN MAT D 
FROM MASBAL (PREDICTED VALUES) 

CAVCLC=AVERAGE CABIN CONC (MG/CU M) 

CF INAL = FINAL INCREMENT CONT CONC (MG/CU M) 
CEQLIB=EQUILIBRIUM CONT CONC (MG/CU M) 

M.REM FOR ALL DEVICES PLACED IN COL 21 OF MAT DD 
OUTPUTS 

TO MAIN PROGRAM 
PUT IN MAT CC 

CAVPRD=PRED CABIN AV CONC (MG/CU M) : =CC(I,2) 
CEQLIB=EQUILIBRIUM CABIN CONT CONC (MG/CU M) : =C 
CFINAL= FINAL CABIN CONT CONC (MG/CU M):=CC(I,'4) 


FOR EACH DEVICE IN MAT DD COL 20 


PUTS REM EFF FROM DD COL 20 IN CC (I , 7-10-13-16 
PUTS M.REM IN CC ( I , 5- 8- 1 1 - 1 4 . . . ) 

TO PRAFIL 

NAME AND SIZE OF MAT DD+FIRST & LAST CO 
TO CNRSUB 

I = CONT NO. 

TN, TNl=CONT INCREMENT FINAL, INIT TIME (H 
CAVPRD=CABIN AVERAGE CONT CONC (MG/CU M) 

DD, NROW, NCOL =NAME & DIMENSIONS OF MAT DD 
CC.NROWl , NCOLl =NAME & DIM OF MAT CC 
CDI , NROW2 , NCOL2=NAME & DIM OF MAT CDI 


LAST COLUMN TO ZERO 


(HRS) 
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59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 


CAV CABIN CALC 


LIN2 =NO . ACTIVE DEVICES IN MAT DD 
TO MASBAL 

I=CONT NO. 

TN , TNI =CONT INCREMENT FINAL, INIT TIME (HRS) 

DD , NROW, NCOL=NAME & DIMENSIONS OF MAT DD 
CALCLC=CALC CABIN AV CONC (MG/CU M) 

CC , NROWl , NCOL1 =NAME & DIM OF MAT CC 
CD I , NROW2 , NCOL2=NAME & DIM OF MAT CDI 
CFINAL=CABIN FINAL CONCENTRATION (MG/CU M) 

CEQLIB=CABIN EQUILIBRIUM CONCENTRATION (MG/CU M) 

LIN=NO . OF CONTAMINANTS IN MAT CDI 
LIN2=NO. ACTIVE DEVICES IN MAT DD 
TO RROUT 

MATRIX NAME, #ROWS,#COLS, FIRST & LAST COL TO PRINT, #LINES TO PRINT 

BASIC TIME INCREMENT (HRS) 

BINC=DD (1,11) 

SET FINAL INCREMENT TIME (HRS) 

TN=0. l*BINC/24 

BEGIN LOOP FOR EACH CONTAMINANT - ONE AT A TIME 

CALCULATE REM EFF FOR EACH DEVICE, GET M.REM, 

CEQULIB , C FINAL- PUT IN MAT CC 
DO 100 1=1, LIN 

ZERO MAT DD COL 17 TO 21 

CALL PRAFIL ( DD , NROW , NCOL , 17,21) 

: SET CAVPRD = MINIMUM VALUE TO ALLOW COMPUTATION 

CAVPRD=lE-20 
DD ( J , 2 2 ) =CAVPRD 

CALC REMOVAL EFFICIENCYS (THROUGH EFF CALLING SUBROUTINE) 

' THIS STORES REM EFF IN MAT DD COL 20 FOR EACH DEVICE 

CALL CNRSUB(I,TN, TNI, DD, NROW, NCOL, CC, NROWl, NCOL1, 

+ CDI , NROW2 , NCOL2 , LIN2 , KK ) 

; CALL MASS BALANCE-GET CAVPRD ( =CAVCLC IN MASBAL) , CEQ, CF INAL, M.REM 

CALL MASBAL ( I , TN , TNI , DD , NROW, NCOL , CC , NROWl , NCOL 1 , 

+ CAVCLC , CDI , NROW2 , NCOL2 , CFINAL , CEQLIB , LIN , LIN2 ) 

CAVPRD=CAVCLC 

; PUT CAVPRD, CEQLIB, & CFINAL IN MAT CC 

CC ( I , 2 ) =CAVPRD 
CC (1,3) =CEQLIB 
CC ( I , 4 ) =C FINAL 

: GET REM EFF FROM DD COL 20 AND PUT IN CC (I, 7-10-13 ETC) 

K=7 

DO 101 J=2 , LIN2 
CC ( I , K) =DD ( J, 20 ) 

K=K+3 

101 CONTINUE 

TAKE M. REMOVED FROM MAT DD COL 21 AND PUT IN MAT CC(I,5-8-...) 

CABIN REMOVAL RATE 
CC (1,5) =DD (1,21) 

: DEVICE 2-15 REMOVAL RATE 

K=8 

DO 102 J=2 , LIN2 
CC ( I , K) =DD ( J, 21 ) 

K=K+3 

102 CONTINUE 

' IF PRTSWl=l THEN PRINT MAT DD+MAT CC INFO FOR THIS CONTAMINANT 
IF (PRTSW1.EQ. 1) THEN 

OPEN ( IMSGDN , FILE= 1 CON 1 , IOSTAT=IOVAL) 

WRITE ( IMSGDN, *) 'PRINTOUT FOR ONE CONT AT END OF PCSET' 

WRITE (IMSGDN, * ) 'INFO FROM MAT CC ' 


C 

C 
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119 

120 
121 
122 

123 

124 

125 

126 C 

127 

128 C 

129 

130 C 

131 

NUMBER 

NUMBER 

NUMBER 

NUMBER 


CLOSE ( IMSGDN ) 

CALL CROUT (NN, CC , NROW1 , NCOL1 , 1 , NCOL1 , LIN, I , I 
OPEN ( IMSGDN, FILE= ' CON' , IOSTAT=IOVAL ) 

WRITE ( IMSGDN, *) 'INFO FROM MAT DD' 

CLOSE (IMSGDN) 

CALL RROUT ( DD, NROW, NCOL , 1 , NCOL, LIN2 , IMSGDN) 
END IF 

100 CONTINUE 
RETURN 

********** END OF SUBROUTINE PCSET ****************** 
END 

OF WARNINGS IN PROGRAM UNIT: 0 
OF ERRORS IN PROGRAM UNIT-: 0 
OF WARNINGS IN COMPILATION : 0 
OF ERRORS IN COMPILATION : 0 


, IMSGDN) 
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Options: /C 80 /L /BY 05/21/92 12:58:44 

****************************************************************** 


RM/FORTRAN Compiler (V2.42) 

Source File: C:\RMFORT\TCC\PRAFIL. 

1 C 

* SUBROUTINE PRAFIL 

* SUBROUTINE TO FILL ADJUSTABLE SIZE REAL ARRAY WITH ZEROS * 

* PARTIAL FILL-FROM COL FSTCOL TO COL LSTCOL * 

* ***************************************************************** 

SUBROUTINE PRAFIL ( X , NROW , NCOL , FSTCOL , LSTCOL ) 

INTEGER NCOL , NROW , FSTCOL , LSTCOL 
REAL X (NROW, NCOL) 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 


100 

110 


INPUTS : 

X , NROW, NCOL=NAME AND DIMENSIONS OF MATRIX X 

FSTCOL, LSTCOL=FIRST AND LAST COLUMN TO FILL WITH ZEROS 

IF ( FSTCOL. GT. NCOL) FSTCOL=NCOL 
IF ( LSTCOL. GT. NCOL) LSTCOL=NCOL 
IF ( FSTCOL . GT . LSTCOL ) FSTCOL=LSTCOL 
DO 110 1=1, NROW 
DO 100 J=FSTCOL, LSTCOL 
X(I,J)=0.0 
CONTINUE 
CONTINUE 
RETURN 

************* END OF SUBROUTINE PRAFIL 
END 


*★*★*★*★**★**★★★★*★***★* 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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RM/ FORTRAN Compiler (V2.42) 

Source File: C:\RMFORT\TCC\PRFANS 
1 C FILE PRFANS 


Options: /C 80 /L /BY 05/21/92 12:58 


58 


2 

3 

4 

5 

6 

7 C 

8 C 

9 

10 

11 

12 

13 C 

14 C 

15 

16 C 

17 C 

18 C 

19 C 


20 
21 
22 

23 

24 

25 

26 C 

27 C 

28 C 

29 

30 

31 

32 

33 

34 

35 C 

36 

37 C 

38 

39 

40 C 

41 C 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 


020 


040 


********************** + ********** i ,* i '*i r i, ir i riri , i , i , i ' i ,„ i , i ' i ' i , i ' i ' iriri ' i ' iriii ' iti ' ir 

* SUBROUTINE PRC DA 

PROGRAM TO PRINT ANSWERS FOR CONCENTRATION DATA * 

************ + ***************** + i , i < i ,* i ,* i , + i ' ir i :iri ,i, i!i ' ili ' i ' iri ' itiri ' i ' ir * iriii ' i! * i ' 

NOTES: (l)FILE MUST BE OPEN BEFORE STARTING THIS SUBROUTINE 
( 2 ) IDEVNO MUST BE 6 FOR FORM FEEDS TO BE PRINTED 

PRCDA < ™ * TN1 ' LIN , CC , NROW1 , NCOL1 , CDI , NROW2 , NCOL2 , NN , 
IDEVNO, NINC, IMONTH, IDAY, I YEAR , IHOUR, IMINUTE , FNAME IOVAL IPGCTR 
PRTSW8 , PRTSW9 , FCPLOT, IDEVN1 ) 

SUBROUTINES REQUIRED: 

HDG1 , HDG2 , HDG3 , HDG4 

TN, TNl=FINAL AND INITIAL INCREMENT TIME (HRS) 

LIN=TOTAL NUMBER OF CONTAMINANTS 
CC , NROW1 , NCOLl=NAME & SIZE OF MAT CC 
CDI , NROW2 , NCOL2 =NAME & SIZE OF MAT CDI 
NN=NAME OF MAT NN 
IDEVNO=DEVICE NUMBER FOR OUTPUT 
NINC=TIME INCREMENT NUMBER 
=0 THEN PRINT HDG3 WITH PCALC 
=-l THEN PRINT HDG3 WITH FINAL 
ELSE PRINT HDG2 WITH INCREMENT NUMBER 
IMONTH .. IMINUTE=TIME AND DATE INFO 
FNAME=FILE NAME OUTPUT DATA IS STORED ON 

IPGCTR=COUNTER FOR SEQUENTIAL PAGE NUMBERS ON ALL OUTPUT 

REAL CC (NROW1 , NCOL1 ) 

REAL CDI (NROW2 , NCOL2 ) 

CHARACTER CNAME*3 0 , FNAME* 2 4 , ECHR* 1 , FCPLOT*24 

CHARACTER NN(NROW1)*30 

INTEGER PRTSW8 , PRTSW9 , NINC 

ECHR= EXCEEDS MAC CHARACTER (Y OR N) 

IF ( (PRTSW8.EQ. 1) .OR. ( (PRTSW8.EQ.0) .AND. (NINC.EQ -1))) THEN 
INCREMENT PAGE COUNTER BY ONE 
I PGCTR= I PGCTR+ 1 

START FIRST PAGE 
PRINT FORM FEED 

WRITE { IDEVNO, 20 , IOSTAT=IOVAL , ERR=900) 

FORMAT ( • 1 • ) 

PRINT BLANK LINE FOLLOWED BY HEADINGS 1 , 3 &4 
WRITE ( IDEVNO, 40 , IOSTAT=IOVAL , ERR=900 ) 

FORMAT (IX, ’ 1 ) 

t^LHDGI (IMONTH, IDAY, I YEAR, IHOUR, IMINUTE, FNAME, IPGCTR, IDEVNO) 
lr (NINC.EQ.0) THEN 

CALL HDG3 (1, TNI, TN, IDEVNO) 

ELSEIF (NINC. EQ. -1) THEN 

CALL HDG3 (2 , TNI, TN, IDEVNO) 

ELSE 

CALL HDG2 (NINC, TNI, TN, IDEVNO) 

ENDIF 

CALL HDG4 (IDEVNO) 

PRINT ANOTHER BLANK LINE 

WRITE ( IDEVNO, 40 , IOSTAT=IOVAL , ERR=900 ) 


89 


59 

60 
61 
62 

63 

64 

65 

66 
67 


BEGIN LOOP FOR EACH CONTAMINANT 1 TO LIN 
DO 100 1=1. LIN 

CNAME=CONTAMINANT NAME 
CNAME=NN ( I ) 

FCONC= FINAL CONT CONCENTRATION (MG CU M) 

FCONC =CC (I 4 ) 

RMAC = MAX I MUM ALLOWABLE C ON C ENTRAT 1 1 ON (MG/CU M) 
RMAC=CDI (1,9) 


68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

HDG1 ( 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 
117 


010 


IF CABIN CONC>MAC PRINT 'Y' OTHERWISE PRINT ‘N 1 
IF (FCONC .GT. RMAC) THEN 
ECHR= ' Y ' 

ELSE 

ECHR= ' N ‘ 

ENDIF 

PRINT 56 LINES OF DATA AND THEN START NEW PAGE 

WRITE (IDEVNO, 10, IOSTAT=IOVAL , ERR=900 ) I , CNAME , FCONC , RMAC , ECHR 
FORMAT (IX, 14 , IX, A, IX, Gil . 4 , 1X,G11. 4, 5X, A) 


C CHECK FOR 56 LINES- IF SO, INCREMENT PAGE NUMBER+ START NEW PAGE 

IF ( INT (REAL ( I ) / 5 6 ) .-EQ . REAL ( I ) / 56 ) THEN 
I PGCTR= I PGCTR+ 1 • 

C START SUBSEQUENT PAGES 

C PRINT FORM FEED 

WRITE (IDEVNO, 50 , IOSTAT=IOVAL, ERR=900 ) 

050 FORMAT ( ' 1 ' ) „ , , . 

C PRINT BLANK LINE FOLLOWED BY HEADINGS 1 , 3 &4 

WRITE (IDEVNO, 30 , IOSTAT=IOVAL, ERR=900) 

030 FORMAT ( IX, *' ) 

CALL 

IMONTH, IDAY, I YEAR, IHOUR, IMINUTE, FNAME, IPGCTR, IDEVNO) 

IF(NINC.EQ.O) THEN 
CALL HDG3 (1, TNI, TN, IDEVNO) 

ELSEIF (NINC . EQ . -1 ) THEN 
CALL HDG3 (2, TNI, TN, IDEVNO) 

ELSE 

CALL HDG2 (NINC, TNI, TN, IDEVNO) 

ENDIF 

CALL HDG4 (IDEVNO) 

C PRINT ANOTHER BLANK LINE 

WRITE (IDEVNO, 30 , IOSTAT=IOVAL, ERR=900 ) 

ENDIF 
100 CONTINUE 
ENDIF 


. ********* WRITE CONCENTRATION DATA TO A PLOT FILE ************* 

IF (NINC.NE.-l) THEN 

IF ( ( PRTSW9 . EQ. 1 ) .OR. (PRTSW9 . EQ. 3 ) ) THEN 

DO 120 1=1, LIN, 300 
IS = I 

IE=I+299 

IF (IE.GT.LIN) IE=LIN 

WRITE ( IDEVNl , 110 , IOSTAT=IOVAL, ERR=900 ) TNI , TN, 

+ (CC ( J , 4 ) , J =IS , IE) 

HO FORMAT (T2 , 2 (F8 . 3 , IX) , 300 (Gil . 4 , : , IX) ) 

120 CONTINUE 
ENDIF 
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118 END IF 

119 GO TO 999 

120 900 WRITE (*,*)' 10 ERROR IN PRCDA= 1 , I OVAL 

121 999 RETURN 

122 END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 


123 

124 C 
12 5 C 

126 C 

127 C 

128 C 

129 

130 C 

131 C 

132 

133 

134 

135 

136 C 

137 C 

138 

139 C 

140 C 

141 C 

142 C 

143 C 

144 C 

145 C 

146 C 

147 C 

148 C 

149 C 

150 C 

151 C 

152 

153 

154 

155 

156 

157 

158 C 

159 

160 

161 C 

162 C 

163 

164 

165 C 

166 

167 

168 

169 

170 

171 

172 


*********** ******************* ********************************** 
SUBROUTINE PRREM1 * 

* SHEET^L T ° PRINT ANSWERS-RATE 0F CONTAMINAMT REMOVAL (MG/HR) * 
******************** *-******************************************* 


NOTES: ( 1 ) FILE MUST BE OPEN BEFORE STARTING THIS SUBROUTINE 
( 2 ) IDEVNO MUST BE 6 FOR FORM FEEDS TO BE PRINTED 

TT ^„, Jp T ^ BR0UTINE PRREM 1 (TN, TNI , LIN, CC ,NROWl ,NCOLl , CDI ,NROW2 NCOL2 NN 
+- IDEVNO, NINC, I MONTH, I DAY, I YEAR, I HOUR, I MINUTE , FNAME , I OVAL, IPGCTR) ’ 

SUBROUTINES REQUIRED: 

HDG1 , HDG2 , HDG3 , HDG8 


TN, TN1=FINAL AND INITIAL INCREMENT TIME (HRS) 

LIN=TOTAL NUMBER OF CONTAMINANTS 
CC , NROW1 , NCOL 1 =NAME & SIZE OF MAT CC 
CDI , NROW2 , NCOL2=NAME & SIZE OF MAT CDI 
NN=NAME OF MAT NN 
IDEVNO=DEVICE NUMBER FOR OUTPUT 
NINC=TIME INCREMENT NUMBER 
=0 THEN PRINT HDG3 WITH PCALC 
=-l THEN PRINT HDG3 WITH FINAL 
ELSE PRINT HDG2 WITH INCREMENT NUMBER 
IMONTH . . IMINUTE=TIME AND DATE INFO 
FNAME=FILE NAME OUTPUT DATA IS STORED ON 
I PGCTR=:COUNTER FOR SEQUENTIAL PAGE NUMBERS ON ALL OUTPUT 

REAL CC (NROW1 , NCOL1 ) 

REAL CDI ( NROW2 , NCOL2 ) 

CHARACTER CNAME * 3 0 , FNAME * 2 4 
CHARACTER NN(NROW1)*30 


INCREMENT PAGE COUNTER BY ONE 
IPGCTR=IPGCTR+1 


START FIRST PAGE 
PRINT FORM FEED 

WRITE ( IDEVNO, 20 , IOSTAT=IOVAL , ERR=900 ) 

020 FORMAT ( 1 1 ' ) 

PRINT BLANK LINE FOLLOWED BY HEADINGS 1, 3 &8 
WRITE ( IDEVNO, 4 0 , IOSTAT=IOVAL , ERR=900 ) ' 

040 FORMAT (IX, ' 1 ) 

CALL HDG1 ( IMONTH , IDAY, I YEAR, IHOUR, IMINUTE, FNAME 
IF (NINC . EQ . 0 ) THEN 

CALL HDG3 (1, TNI, TN, IDEVNO) 

ELSEIF (NINC. EQ. -1 ) THEN 

CALL HDG3 (2 , TNI , TN, IDEVNO) 


IPGCTR, IDEVNO) 
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173 

174 

175 

176 

177 C 

178 
17 9 

180 C 

181 
182 

183 C 

184 

185 

186 

187 

188 

189 C 

190 

191 

192 C 

193 C 

194 

195 

196 C 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 C 

209 

210 
211 
212 

213 

214 

215 

216 
217 


CALL HDG2 (NINC , TNI , TN , IDEVNO) 

END IF 

CALL HDG8 (IDEVNO) 

PRINT ANOTHER BLANK LINE 

WRITE ( IDEVNO , 40 , IOSTAT=IOVAL , ERR=900 ) 


BEGIN LOOP FOR EACH CONTAMINANT 1 TO LIN 
DO 100 1=1, LIN 


PRINT 56 LINES OF DATA AND THEN START NEW PAGE 
WRITE (IDEVNO, 10, IOSTAT=IOVAL , ERR=900 ) I,NN(I), 

+CC (1,5) ,CC(I,8) ,CC( 1,11) ,CC( 1,14) ,CC( 1,17) ,CC( 1,20) ,CC( I, 23) 

+CC (1,26) 

010 FORMAT(lX, 14 , IX, A, 8 (IX, Gil .4) ) 


CHECK FOR 56 LINES- IF SO, INCREMENT PAGE NUMBER+START NEW PAG 
I F ( INT ( REAL ( I ) / 5 6 ) . EQ . REAL ( I ) / 56 ) THEN 
IPGCTR=IPGCTR+1 
START SUBSEQUENT PAGES 
PRINT FORM FEED 

WRITE (IDEVNO, 50 , IOSTAT=IOVAL , ERR=900 ) 


050 

030 

CALL 


FORMAT Cl') , , . Q 

PRINT BLANK LINE FOLLOWED BY HEADINGS 1 , 3 &8 

WRITE (IDEVNO, 30 , IOSTAT=IOVAL, ERR=900) 

FORMAT ( IX 1 1 ) 

HDGl (IMONTH, IDAY, IYEAR, IHOUR, IMINUTE , FNAME , IPGCTR, IDEVNO) 
IF (NINC . EQ. 0) THEN 
CALL KDG3 (1, TNI, TN, IDEVNO) 

ELSEIF (NINC . EQ. -1 ) THEN 
CALL HDG3 (2, TNI, TN, IDEVNO) 

ELSE 

CALL HDG2 (NINC, TNI, TN, IDEVNO) 

END IF 

CALL HDG8 (IDEVNO) 

PRINT ANOTHER BLANK LINE 

WRITE ( IDEVNO, 30 , IOSTAT=IOVAL, ERR=900 ) 


END IF 
100 CONTINUE 

GO TO 999 

900 WRITE (*,*)' IO ERROR IN PRREM1= 1 , IOVAL 
999 RETURN 
END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 


218 

219 C 

220 C 

221 C 

222 C 

223 C 

224 

225 C 

226 C 

227 

228 


**************************** 
* SUBROUTINE PRREM2 


************************************ 

★ 


PROGRAM TO PRINT ANSWERS-RATE OF CONTAMINANT REMOVAL (MG/HR) * 


* SHEET 2 

********************* 


*********************** 


******************** 


MOTES * ( 1 ) FILE MUST BE OPEN BEFORE STARTING THIS SUBROUTINE 
(2 (IDEVNO MUST BE 6 FOR FORM FEEDS TO BE PRINTED 

SUBROUTINE PRREM2 (TN, TNI , LIN, CC , NROWl , NCOLl , CDI , NROW2 , NCOL2 , NN, 
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229 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

274 

275 

276 

277 

278 

279 

280 
281 
282 

283 

284 

285 

286 

287 

288 
289 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 


020 


C 


040 


+ IDEVNO , NINC , IMONTH, IDAY, I YEAR, IHOUR, IMINUTE, FNAME, IOVAL, IPGCTR) 230 
SUBROUTINES REQUIRED: 

HDG1 , HDG2 , HDG3 , HDG9 

TN, TNI = FINAL AND INITIAL INCREMENT TIME (HRS) 

LIN=TOTAL NUMBER OF CONTAMINANTS 
CC , NROW1 , NCOLl=NAME & SIZE OF MAT CC 
CDI , NROW2 , NCOL2 =NAME & SIZE OF MAT CDI 
NN=NAME OF MAT NN 
IDEVNO=DEVICE NUMBER FOR OUTPUT 
NINC=TIME INCREMENT NUMBER 

=0 THEN PRINT HDG3 WITH PCALC 
=-l THEN PRINT HDG3 WITH FINAL 
ELSE PRINT HDG2 WITH INCREMENT NUMBER 
IMONTH .. IMINUTE=TIME AND DATE INFO 
FNAME=FILE NAME OUTPUT DATA IS STORED ON 

I PGCTR=COUNTER FOR SEQUENTIAL PAGE NUMBERS ON ALL OUTPUT 

REAL CC (NROW1 , NCOL1 ) 

REAL CD I ( NROW2 , NCOL2 ) 

CHARACTER CNAME *30, FNAME * 2 4 
CHARACTER NN(NROW1)*30 

INCREMENT PAGE COUNTER BY ONE 
I PGCTR= I PGCTR+ 1 

START FIRST PAGE 

DON'T PRINT FORM FEED UNLESS NO. CONT >20 
I F ( LIN . GT .20) THEN 
PRINT FORM FEED 

WRITE ( IDEVNO , 20 , IOSTAT=IOVAL , ERR=900 ) 

FORMAT ( ' 1 ' ) 

END IF 

PRINT BLANK LINE FOLLOWED BY HEADINGS 1, 

WRITE ( IDEVNO, 40, IOSTAT= IOVAL , ERR=9 00 ) 

FORMAT (IX, ' • ) 

CALL HDG1 (IMONTH, IDAY, I YEAR, IHOUR, IMINUTE , FNAME , IPGCTR, IDEVNO) 
IF (NINC. EQ. 0) THEN 
CALL HDG3 (1, TNI ,TN, IDEVNO) 

ELSEIF (NINC. EQ. -1) THEN 
CALL HDG3 (2 , TNI, TN, IDEVNO) 

ELSE 

CALL HDG2 (NINC, TNI , TN, IDEVNO) 

END IF 

CALL HDG9 (IDEVNO) 

PRINT ANOTHER BLANK LINE 

WRITE ( IDEVNO, 40 , IOSTAT=IOVAL , ERR=900 ) 


C 


3 &9 


BEGIN LOOP FOR 
DO 100 1=1, LIN 


EACH CONTAMINANT 1 TO LIN 


: PRINT 56 LINES OF DATA AND THEN START 

WRITE (IDEVNO, 10, IOSTAT=IOVAL , ERR=900 ) 
+CC (1,29) , CC ( 1 , 3 2 ) , CC { 1 , 3 5 ) , CC ( 1 , 3 8 ) , CC (1,41) 
010 FORMAT (IX, 14 , IX, A, 7 (IX, Gil .4) ) 


NEW PAGE 
I , NN ( I ) , 

, CC ( 1 , 4 4 ) 


CC(I, 47) 


C CHECK FOR 56 LINES- IF SO, INCREMENT PAGE NUMBER+START NEW PAGE 
IF ( INT (REAL ( I ) /56) . EQ. REAL ( I ) / 56 ) THEN 
IPGCTR=IPGCTR+1 

C START SUBSEQUENT PAGES 
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290 

C 

291 


292 

050 

293 

C 

294 


295 

030 

296 


297 


298 


299 


300 


301 


302 


303 


304 


305 

C 

306 


307 


308 


309 

100 

310 


311 


312 

900 

313 

999 

314 



PRINT FORM FEED 

WRITE (IDEVNO, 50 , IOSTAT= IOVAL , ERR=900 ) 

FORMAT Cl') n n 

PRINT BLANK LINE FOLLOWED BY HEADINGS 1, 3 &8 

WRITE (IDEVNO, 30 , IOSTAT=IOVAL, ERR=900 ) 

FORMAT (IX 1 1 ) 

CALL HDG1 (IMONTH, IDAY, I YEAR, IHOUR, IMINUTE , FNAME , IPGCTR, IDEVNO) 
IF(NINC.EQ.O) THEN 
CALL HDG3 (1, TNI, TN, IDEVNO) 

ELSEIF (NINC . EQ . - 1 ) THEN 
CALL HDG3 (2, TNI, TN, IDEVNO) 

ELSE 

CALL HDG2 (NINC, TNI, TN, IDEVNO) 

END IF 

CALL HDG9 (IDEVNO) 

PRINT ANOTHER BLANK LINE 

WRITE ( IDEVNO, 30 , IOSTAT=IOVAL, ERR=900 ) 

END IF 
CONTINUE 

GO TO 999 

WRITE ( * , * ) 1 10 ERROR IN -PRREM2= ' , IOVAL 
RETURN 
END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 


315 

316 C 

317 C 

318 C 

319 C 

320 C 

321 

322 C 

323 C 

324 

325 

326 

327 

328 C 

329 C 

330 

331 C 

332 C 

333 C 

334 C 

335 C 

336 C 

337 C 

338 C 

339 C 

340 C 

341 C 

342 C 

343 C 

344 

345 


**************************************************************** 

* SUBROUTINE PRMAS1 * 

* PROGRAM TO PRINT ANSWERS-SUM OF CONT REMOVED BY DEVICE (MG) ^ 

* SHEET 1 

**************************************************************** 


NOTES- ( 1 ) FILE MUST BE OPEN BEFORE STARTING THIS SUBROUTINE 
(2) IDEVNO MUST BE 6 FOR FORM FEEDS TO BE PRINTED 

SUBROUTINE PRMAS1 (TN, TNI , LIN, CC ,NROWl ,NCOLl , CDI , NROW2 , NCOL2 , NN, 
+ IDEVNO, NINC, IMONTH, IDAY, IYEAR, IHOUR, IMINUTE, FNAME, IOVAL, IPGCTR) 


SUBROUTINES REQUIRED: 

HDG1 , HDG2 , HDG3 , HDG5 

TN, TNI =F INAL AND INITIAL INCREMENT TIME (HRS) 

LIN=TOTAL NUMBER OF CONTAMINANTS 
CC , NROWl , NCOLl=NAME & SIZE OF MAT CC 
CDI , NROW2 , NCOL2 =NAME & SIZE OF MAT CDI 
NN=NAME OF MAT NN 
IDEVNO=DEVICE NUMBER FOR OUTPUT 
NINC=TIME INCREMENT NUMBER 

=0 THEN PRINT HDG3 WITH PCALC 
=-l THEN PRINT HDG3 WITH FINAL 
ELSE PRINT HDG2 WITH INCREMENT NUMBER 
IMONTH .. IMINUTE=TIME AND DATE INFO 
FNAME=FILE NAME OUTPUT DATA IS STORED ON 
IPGCTR=COUNTER FOR SEQUENTIAL PAGE NUMBERS ON ALL OUTPUT 

REAL CC (NROWl ,NCOLl) 
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346 

347 

348 

349 

350 C 

351 

352 

353 C 

354 C 

355 

356 

357 C 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 C 

370 

371 

372 C 

373 

374 

375 C 

376 

377 

378 

379 

380 

381 C 

382 

383 

384 C 

385 C 

386 

387 

388 C 

389 

390 

391 

392 

393 

394 

395 

396 

397 

398 

399 

400 C 

401 

402 

403 

404 

405 


REAL CD I ( NROW2 , NCOL2 ) 

CHARACTER CNAME*3 0 , FNAME*24 
CHARACTER NN(NROWl)*30 

INCREMENT PAGE COUNTER BY ONE 
I PGCTR= I PGCTR+ 1 

START FIRST PAGE 
PRINT FORM FEED 

WRITE ( IDEVNO ,20, IOSTAT=IOVAL , ERR=900 ) 

020 FORMAT ( 1 1 * ) 

PRINT BLANK LINE FOLLOWED BY HEADINGS 1 , 3 &5 
WRITE ( IDEVNO, 40 , IOSTAT=IOVAL , ERR=900 ) 

040 FORMAT (IX, ' • ) 

CALL HDG1 (IMONTH, IDAY, IYEAR , IHOUR, IMINUTE , FNAME , IPGCTR, IDEVNO) 
IF(NINC.EQ.O) THEN • 

CALL HDG3 ( 1 , TNI , TN , IDEVNO ) 

ELSEIF (NINC. EQ. -1) THEN 

CALL HDG3 (2 , TNI ,TN, IDEVNO) 

ELSE 

CALL HDG2 (NINC, TNI, TN, IDEVNO) 

END IF 

CALL HDG5 (IDEVNO) 

PRINT ANOTHER BLANK LINE 

WRITE (IDEVNO, 40, IOSTAT=IOVAL , ERR=900 ) 

BEGIN LOOP FOR EACH CONTAMINANT 1 TO LIN 
DO 100 1=1, LIN 

PRINT 56 LINES OF DATA AND THEN START NEW PAGE 
WRITE ( IDEVNO, 10, IOSTAT=IOVAL , ERR=900 ) I,NN(I) , 

+CC(I,6) ,CC( 1,9) ,CC (1,12) ,CC( 1,15) , CC ( 1,18) ,CC(I, 21) , CC (1,24) , 

+CC (1,27) 

010 FORMAT (IX, 14, IX, A, 8 (IX, Gil. 4) ) 

CHECK FOR 56 LINES-IF SO, INCREMENT PAGE NUMBER+START NEW PAGE 
IF ( INT ( REAL ( I ) / 5 6 ) . EQ . REAL ( I ) /56) THEN 
IPGCTR=IPGCTR+1 
START SUBSEQUENT PAGES 
PRINT FORM FEED 

WRITE ( IDEVNO , 50, IOSTAT=IOVAL, ERR=900 ) 

050 FORMAT ( ' 1 ' ) 

PRINT BLANK LINE FOLLOWED BY HEADINGS 1, 3 &8 
WRITE (IDEVNO, 30 , IOSTAT=IOVAL , ERR=900 ) 

030 FORMAT (IX, ' ' ) 

CALL HDG1 (IMONTH, IDAY, IYEAR, IHOUR, IMINUTE , FNAME , IPGCTR, IDEVNO) 
IF (NINC . EQ. 0 ) THEN 
CALL HDG3 ( 1 , TNI, TN, IDEVNO) 

ELSEIF (NINC. EQ, -1) THEN 
CALL HDG3 (2, TNI, TN, IDEVNO) 

ELSE 

CALL HDG2 (NINC, TNI, TN, IDEVNO) 

END IF 

CALL HDG5( IDEVNO) 

PRINT ANOTHER BLANK LINE 

WRITE (IDEVNO, 30 , IOSTAT=IOVAL, ERR=900) 

ENDIF 
100 CONTINUE 
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406 

407 

408 

409 

NUMBER 

NUMBER 

410 

411 C 

412 C 

413 C 

414 C 

415 C 

416 

417 C 

418 C 

419 

420 

421 

422 

423 C 

424 C 

425 

426 C 

427 C 

428 C 

429 C 

430 C 

431 C 

432 C 

433 C 

434 C 

435 C 

436 C 

437 C 

438 C 

439 

440 

441 

442 

443 

444 

445 C 

446 

447 

448 C 

449 C 

450 

451 C 

452 

453 

454 

455 C 

456 

457 

458 

459 

460 

461 


GO TO 999 

900 WRITE (*,*)' 10 ERROR IN PRMAS1= ' , IOVAL 
999 RETURN 
END 

OF WARNINGS IN PROGRAM UNIT: 0 
OF ERRORS IN PROGRAM UNIT: 0 


★ ★*★★***********★**★★★★★*★★*★**★★*★★★*****★*********★★★*★*★***** 

* SUBROUTINE PRMAS2 * 

* PROGRAM TO PRINT ANSWERS- SUM OF CONT REMOVED BY DEVICE (MG) * 

* SHEET 2 * 

* ★★★★★★★★★★★★★★★★★★★★★★★a*************************************** 


NOTES: ( 1 ) FILE MUST BE OPEN BEFORE STARTING THIS SUBROUTINE 
(2 ) IDEVNO MUST BE 6 FOR FORM FEEDS TO BE PRINTED 

SUBROUTINE PRMAS2 ( TN , TNI , LIN , CC , NROW1 , NCOLl , CDI , NROW2 , NCOL2 , NN , 
+IDEVNO, NINC , IMONTH, IDAY , IYEAR, IHOUR, I MINUTE , FNAME , IOVAL, IPGCTR) 

SUBROUTINES REQUIRED: 

HDG1 , HDG2 , HDG3 , HDG6 

TN , TNl=FINAL AND INITIAL INCREMENT TIME (HRS) 

LIN=TOTAL NUMBER OF CONTAMINANTS 
CC,NROWl , NCOLl =NAME & SIZE OF MAT CC 
CDI , NROW2 , NCOL2=NAME & SIZE OF MAT CDI 
NN=NAME OF MAT NN 
I DEVNO= DEVICE NUMBER FOR OUTPUT 
NINC=TIME INCREMENT NUMBER 

=0 THEN PRINT HDG3 WITH PCALC 
=-l THEN PRINT HDG3 WITH FINAL 
ELSE PRINT HDG2 WITH INCREMENT NUMBER 
IMONTH . . IMINUTE=TIME AND DATE INFO 
FNAME=FILE NAME OUTPUT DATA IS STORED ON 
I PGCTR=COUNTER FOR SEQUENTIAL PAGE NUMBERS ON ALL OUTPUT 

REAL CC (NROW1 , NCOLl) 

REAL CDI (NROW2 , NCOL2 ) 

CHARACTER CNAME *30, FNAME * 2 4 
CHARACTER NN(NROWl)*30 

INCREMENT PAGE COUNTER BY ONE 
I PGCTR= I PGCTR+ 1 

START FIRST PAGE 

DON'T PRINT FORM FEED UNLESS NO. CONT >20 
IF (LIN. GT. 20) THEN 
PRINT FORM FEED 

WRITE ( IDEVNO ,20, IOSTAT= IOVAL , ERR=90 0 ) 

020 FORMAT ( ' 1 ' ) 

END IF 

PRINT BLANK LINE FOLLOWED BY HEADINGS 1 , 3 &9 
WRITE ( IDEVNO ,40, IOSTAT= IOVAL , ERR=9 00 ) 

040 FORMAT (IX, ' ' ) 

CALL HDG1 (IMONTH, IDAY, IYEAR, IHOUR, IMINUTE , FNAME , IPGCTR, IDEVNO) 
IF(NINC.EQ.O) THEN 
CALL HDG3 (1 ,TN1 ,TN, IDEVNO) 

ELSEIF (NINC. EQ. -1) THEN 
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462 

463 

464 

465 

466 

467 

468 

469 

470 

471 

472 

473 

474 

475 

476 

477 

478 

479 

480 

481 

482 

483 

484 

485 

486 

487 

488 

489 

490 

491 

492 

493 

494 

495 

496 

497 

498 

499 

500 

501 

502 

503 

504 

505 

506 


CALL HDG3 (2 , TNI , TN, IDEVNO) 

ELSE 

CALL HDG2 (NINC , TNI , TN, IDEVNO) 

ENDIF 

CALL HDG6( IDEVNO) 

C PRINT ANOTHER BLANK LINE 

WRITE ( IDEVNO, 40, IOSTAT=IOVAL, ERR=900 ) 

C BEGIN LOOP FOR EACH CONTAMINANT 1 TO LIN 

DO 100 1=1 , LIN 

C PRINT 56 LINES OF DATA AND THEN START NEW PAGE 

WRITE (IDEVNO, 10, IOSTAT=IOVAL , ERR=900 ) I,NN(I) , 

+CC (1,30) , CC(I,33) ,CC( 1,36) , CC ( 1,39) , CC (1,42) , CC (1,45) , CC ( I, 48) 
010 FORMAT (IX, 14 , IX , A, 7 ( IX , Gil . 4 ) ) 


C 


c 

c 


050 

C 


030 


C 


CHECK FOR 56 LINES- IF SO, INCREMENT PAGE NUMBER+START NEW PAGE 
I F ( INT ( REAL ( I ) / 56 ) . EQ . REAL ( I ) /56) THEN 
IPGCTR=IPGCTR+1 
START SUBSEQUENT PAGES 
PRINT FORM FEED 

WRITE ( IDEVNO, 50 , IOSTAT=IOVAL , ERR=900 ) 

FORMAT ( 1 1 ’ ) 

PRINT BLANK LINE FOLLOWED BY HEADINGS 1, 3 &8 
WRITE ( IDEVNO, 30 , IOSTAT=IOVAL, ERR=900 ) 

FORMAT (IX, ' 1 ) 

CALL HDG1 (IMONTH, I DAY , I YEAR, IHOUR, IMINUTE, FNAME , IPGCTR, IDEVNO) 
IF (NINC. EQ. 0) THEN 
CALL HDG3 (1, TNI, TN, IDEVNO) 

ELSEIF (NINC . EQ . - 1 ) THEN 
CALL HDG3 (2 , TNI, TN, IDEVNO) 

ELSE 

CALL HDG2 (NINC, TNI , TN, IDEVNO) 

ENDIF 

CALL HDG6 (IDEVNO) 

PRINT ANOTHER BLANK LINE 

WRITE (IDEVNO, 30 , IOSTAT=IOVAL, ERR=900) 


ENDIF 

100 CONTINUE 


GO TO 999 

900 WRITE (*,*)' IO ERROR IN PRMAS2= 1 , IOVAL 
999 RETURN 
END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 

507 

508 C ************************************************************** 

509 C * SUBROUTINE PREFF 

510 C * PROGRAM TO PRINT ANSWERS-END OF INCREMENT REMOVALL EFF (DEC) 

511 C * 

512 C ************************************************************** 

513 

514 C NOTES: (1) FILE MUST BE OPEN BEFORE STARTING THIS SUBROUTINE 

515 C (2) IDEVNO MUST BE 6 FOR FORM FEEDS TO BE PRINTED 

516 

517 SUBROUTINE PREFF ( TN , TNI , LIN , CC , NROW1 , NCOL1 , CDI , NROW2 , NCOL2 , NN , 



518 +IDEVNO, NINC , IMONTH, IDAY, IYEAR , IHOUR, IMINUTE , FNAME , IOVAL, IPGCTR, 

519 +PRTSW 8 , PRTSW9, IDEVN2) 

520 

521 C SUBROUTINES REQUIRED: 

522 C HDG1 , HDG2 , HDG3 , HDG7 

523 

524 C TN,TN1=FINAL AND INITIAL INCREMENT TIME (HRS) 

525 C LIN=TOTAL NUMBER OF CONTAMINANTS 

526 C CC , NROWl , NCOLl=NAME & SIZE OF MAT CC 

527 C CD I , NROW2 , NCOL2 =NAME & SIZE OF MAT CDI 

528 C NN=NAME OF MAT NN 

529 C IDEVNO=DEVICE NUMBER FOR OUTPUT 

530 C NINC=TIME INCREMENT NUMBER 

531 C =0 THEN PRINT HDG3 WITH PCALC 

532 C =-l THEN PRINT HDG3 WITH FINAL 

533 C ELSE PRINT HDG2 WITH INCREMENT NUMBER 

534 C IMONTH. . IMINUTE=TIME AND DATE INFO 

535 C FNAME=FILE NAME OUTPUT DATA IS STORED ON 

536 C I PGCTR=COUNTER FOR SEQUENTIAL PAGE NUMBERS ON ALL OUTPUT 

537 

538 REAL CC (NROWl , NCOL1 ) 

539 REAL CDI (NROW2 , NCOL2 ) 

540 CHARACTER CNAME*30 , FNAME* 24 

541 CHARACTER NN ( NROWl) *30 

542 INTEGER PRTSW 8 , PRTSW9 , IDEVN2 , IDEVNO , I , J , K , H , NINC 

543 IF { (PRTSW 8 .EQ. 1) .OR. ( (PRTSW8.EQ.0) .AND. (NINC.EQ. -1) ) ) THEN 

544 C INCREMENT PAGE COUNTER BY ONE 

545 I PGCTR= I PGCTR+ 1 

546 

547 C START FIRST PAGE 

548 C PRINT FORM FEED 

549 WRITE (IDEVNO, 20, IOSTAT=IOVAL,ERR=900) 

550 020 FORMAT { ' 1 ’ ) 

551 C PRINT BLANK LINE FOLLOWED BY HEADINGS 1 , 3 &9 

552 WRITE ( IDEVNO, 40 , IOSTAT=IOVAL , ERR=900 ) 

553 040 FORMAT ( IX , ' ' ) 

554 CALL HDG1 (IMONTH, IDAY, IYEAR, IHOUR, IMINUTE, FNAME, IPGCTR, IDEVNO) 

555 IF (NINC.EQ. 0) THEN 

556 CALL HDG3 (1 ,TN1 ,TN, IDEVNO) 

557 ELSE IF (NINC . EQ . - 1 ) THEN 

558 CALL HDG3 (2 , TNI, TN, IDEVNO) 

559 ELSE 

560 CALL HDG2 (NINC, TNI , TN, IDEVNO) 

561 END IF 

562 CALL HDG7 ( IDEVNO) 

563 C PRINT ANOTHER BLANK LINE 

564 WRITE (IDEVNO, 40, IOSTAT=IOVAL , ERR=900 ) 

565 

566 C BEGIN LOOP FOR EACH CONTAMINANT 1 TO LIN 

567 DO 100 1=1, LIN 

568 

569 C PRINT 56 LINES OF DATA AND THEN START NEW PAGE 

570 WRITE (IDEVNO, 10, IOSTAT= IOVAL, ERR=9 00) I,NN(I), 

571 +CC(I,7) ,CC< 1,10) ,CC (1,13) , CC (1,16) ,CC (1,19) ,CC (1,22) , CC ( 1,25) , 

572 +CC(I,28) ,CC( 1,31) ,CC( 1,34) , CC ( 1,37) ,CC< 1,40) , CC ( 1,43) , CC ( I, 46) 

573 010 FORMAT (IX, 14 , IX, A, 14 (IX, F5 . 3 ) ) 

574 

575 c CHECK FOR 56 LINES-IF SO, INCREMENT PAGE NUMBER+START NEW PAGE 

576 IF ( INT (REAL ( I ) / 5 6 ) . EQ. REAL ( I ) /56 ) THEN 

577 IPGCTR=IPGCTR+1 
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578 

579 

580 

581 

582 

583 

584 

585 

586 

587 

588 

589 

590 

591 

592 

593 

594 

595 

596 

597 

598 

599 

600 
601 
602 

603 

604 

605 

606 

607 

608 

609 

610 
611 
612 

613 

614 

615 

616 

617 

618 

619 

620 


C START SUBSEQUENT PAGES 

C PRINT FORM FEED 

WRITE ( IDEVNO , 50 , IOSTAT=IOVAL , ERR=9 00 ) 

050 FORMAT ('1') 

C PRINT BLANK LINE FOLLOWED BY HEADINGS 1 , 3 &8 

WRITE (IDEVNO, 3 0 , IOSTAT= IOVAL , ERR=900 ) 

030 FORMAT (IX, ' ' ) 

CALL HDG1 (IMONTH, IDAY, I YEAR , IHOUR, IMINUTE , FNAME , IPGCTR, IDEVNO) 
IF(NINC.EQ.O) THEN 
CALL HDG3 (1, TNI, TN, IDEVNO) 

ELSEIF (NINC . EQ . - 1 ) THEN 
CALL HDG3 (2, TNI, TN, IDEVNO) 

ELSE 

CALL HDG2 (NINC, TNI, TN, IDEVNO) 

END IF 

CALL HDG7 ( IDEVNO) 

C PRINT ANOTHER BLANK LINE 

WRITE (IDEVNO, 30 , IOSTAT=IOVAL , ERR=900 ) 

END IF 

100 CONTINUE 
END IF 

c *************** WRITE DATA TO A FILE FOR PLOTTING ************* 

IF (NINC.NE.-l) THEN 

IF ( (PRTSW9.EQ.2) .OR. (PRTSW9.EQ.3) ) THEN 
DO 70 1=1 , LIN, 300 
IS=I 

IE=I+299 

IF (IE.GT.LIN) IE=LIN 
K=7 

DO 60 H=2 , 15 

WRITE (IDEVN2, 55, IOSTAT=IOVAL,ERR=900) TNI , TN, H , 

+ (CC (J,K) , J = IS, IE) 

055 FORMAT (T2 , 2 ( F8 . 2 , IX) , 12 , IX , 3 00 (F5 . 3 , : , IX) ) 

K=K+3 
060 CONTINUE 
070 CONTINUE 
END IF 
END IF 
GO TO 999 

900 WRITE (*,*) 1 IO ERROR IN PREFF= 1 , IOVAL 
999 RETURN 
END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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RM/FORTRAN Compiler (V2.42) 

Source File: C:\RMFORT\TCC\PREDCT. Options: /C 80 /L /BY 05/21/92 12:58:48 

Q ***************************************************************** 

2 C * SUBROUTINE PREDCT * 

3 C * BASED ON REMOVAL EFF & SUM MASS REMOVED OF LAST INCREMENT, * 

4 C * AND M.GEN OF THIS INCREMENT, PREDICT CAV PRED * 

5 C * (CEQUILIB , CFINAL, M.REM ARE ALSO CALC, BUT NOT NEEDED) 

6 C * WORKS FOR ONE CONT AT A TIME * 

J Q *****★*******★****** + *★***★*★*★**★*★*********★******************* 

8 

9 SUBROUTINE PREDCT ( I , TN , TN 1 , CAVPRD , DD , NROW , NCOL , 

10 +CC , NROW1 , NCOLl , CDI , NROW2 , NCOL2 , CAVCLC , CFINAL , CEQLIB , LIN, LIN2 , NN) 

11 

12 INTEGER NROW, NCOL , NROWl , NCOLl , NROW2 , NCOL2 

13 CHARACTER NN ( NROWl) *30 

14 REAL DD (NROW, NCOL) 

15 REAL CC (NROWl, NCOLl) 

16 REAL CDI ( NROW2 , NCOL2 ) 

17 C 

18 C SUBROUTINES REQUIRED: 

19 C PRAFIL-ZERO MAT DD COL 17-21 

20 C LODEFF-LOAD REM EFF FOR LAST INCR FROM MAT CC INTO MAT DD COL 20 

21 C MASBAL-CALC CAV PRED BASED ON REM EFF OF LAST INC & M.GEN OF THIS INC 

22 C 

23 C INPUTS: 

24 C FROM MCALC 

25 C I=CONTAMINANT NO. 

26 C TN, TNI = INCREMENT END & BEGINNING TIME (HRS) 

27 C DD, NROW, NCOL =NAME & DIM OF MAT DD 

28 C CC, NROWl, NCOLl =NAME & DIM OF MAT CC 

29 C CDI , NROW2 , NCOL2=NAME & DIM OF MAT CDI 

30 C LIN=NO. OF CONTAMINANTS IN MAT CDI 

31 C LIN2 =NO. DEVICES IN MAT DD 

32 C NN=NAME OF MAT NN 

33 C FROM LODEFF 

34 C TAKES REM EFF FOR LAST INCR (IN MAT CC) AND PUTS IT IN 

35 C MAT DD COL 20 (FOR ALL DEVICES) 

36 C FROM MASBAL 

37 C CAVCLC =CALC CABIN CONT CONC (MG/CU M) 

38 C OUTPUTS: 

39 C TO LODEFF 

40 C I=CONTAMINANT LINE NUMBER IN MAT CC 

41 C DD , NROW , NCOL=NAME & DIMENSIONS OF MAT DD 

42 C CC, NROWl , NCOLl =NAME & DIMENSIONS OF MAT CC 

43 C LIN2 =NO . OF DEVICES IN MAT DD 

44 C TO MASBAL 

45 C TN, TNI = INCREMENT END & BEGINNING TIME (HRS) 

46 C CNVERR=CONVERGENCE ERROR 

47 C DD , NROW, NCOL=NAME & DIM OF MAT DD 

48 C CC, NROWl, NCOLl =NAME & DIM OF MAT CC 

49 C CDI , NROW2 , NCOL2 =NAME & DIM OF MAT CDI 

50 C LIN=NO. OF CONTAMINANTS IN MAT CDI 

51 C LIN2 =NO . DEVICES IN MAT DD 

52 C CVOL=CABIN VOLUME (CU M) 

53 C TCABIN=CAB IN TEMP (DEG K) 

54 C CINIT= INCR INIT CABIN CONT CONC (MG/CU M)=CC(I,1) 

55 C TO MCALC 

56 C CAVCLC 

57 C 

58 C ZERO MAT DD COL 17-21 
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CALL PRAFIL ( DD , NROW, NCOL , 17,21) 


59 

60 

61 C LOAD REM EFF FROM LAST TIME INCR FROM MAT CC INTO MAT DD COL 20 

62 CALL LODEFF ( I , DD , NROW, NCOL , CC , NROWl , NCOL1 , LIN2 ) 

63 

64 C FIND CAV PRED FOR THESE REMOVAL EFFICIENCIES 

65 CALL MASBAL { I, TN, TNI, DD, NROW, NCOL, CC, NROWl, NCOL 1, 

66 + CAVPRD , CDI , NROW2 , NCOL2 , CFINAL , CEQLIB , LIN , LIN2 ) 

67 C 

68 C SET CAV IN PRED DD{I,22)=CAV IN CALC DD(I,17) 

69 DO 100 J=1 , LIN2 

70 DD ( J , 2 2 ) =DD ( J , 17) 

71 100 CONTINUE 

72 C 

73 RETURN 

74 C ************ op SUBROUTINE PREDCT ************************** 

7 5 END 

NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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RM/ FORTRAN Compiler (V2.42) 

Source File: C:\FORTRAN\TCC\RCHBD. Options: /C 80 /L /BIJY 03/15/94 10:41:07 


2 ^ 

2 C * FILE: RCHBD. FOR * 

3 C * SUBROUTINE FOR REM EFF-RADIAL FLOW CHARCOAL BED * 

4 C * DOESNT ALLOW FOR DESORPTION * 

5 C * ASSUMES RELATIVELY THIN BED (OD CLOSE TO ID) * 

7 C 

8 SUBROUTINE RCHBD <TN, TNI , CIIN, TCABIN, COEXIS , BEDQ , EMAX, CARTL, 

9 +BEDOD , BEDID , DENCH , TRTTYP , DCONT , VMOL , MW, VCONC , SOL , SMR , EFF , RH ) 

10 C OUTPUT: 

11 C EFF=BED REMOVAL EFF (DEC) 

12 C INPUTS: 

13 C TN , TNI = INCREMENT INITIAL AND FINAL TIMES (HR) 

14 C Cl IN=BED INLET CONT CONC (MG/CU M) 

15 C TCABIN=CABIN TEMP (DEG K) 

16 C COEXIS=COEXISTANCE FACTOR 

17 C BEDQ=BED FLOW RATE (CU M/HR) 

18 C EMAX =MAX IMUM BED EFF (DEC) 

19 C CARTL=CARTRIDGE LENGTH (M) 

20 C BEDOD=BED OUTSIDE DIAMETER (M) 

21 C BEDID=BED INSIDE DIAMETER (M) 

22 C DENCH=DENS ITY OF CHARCOAL IN BED (KG/CU M) 

23 C TRTTYP=BED TREATMENT TYPE ( 1=CI CHAR,2=PHOS ACID, OTHER #=NONE) 

24 C DCONT=CONT LIQUID DENSITY (GM/CC) 

25 C VMOL=CONT MOLAR VOL (GM/CC) 

26 C MW=CONT MOLECULAR WGT 

27 C VCONC = CONT VAPOR CONCENTRATION AT TCABIN (MG/CU M) 

28 C SOL= HENRY 1 S LAW CONSTANT FOR WATER SOLUBILITY 

29 C (ATM/MOL FRACTION) 

30 C SMR=SUM OF CONT MASS STORED IN BED (MG) -FROM LAST INCR 

31 C 

32 REAL LPREV , LAVNl , LUTIL, LIMM, LAVAV, LADS , MW 

33 INTEGER FACID, FCI 

34 C 

35 C SET CIN=CI IN (THIS PREVENTS CIN FROM BEING PASSED BACK UP 

36 C TO OTHER SUBROUTINES IF IT IS SET TO IE-20) 

37 CIN=CI IN 

38 

39 C BED TREATMENT LOGIC 

40 C FAC I D= FLAG IF BED IS TREATED WITH PHOSPHORIC ACID (Y=l N=0) 

41 C FCI=FLAG FOR Cl CHAR IN BED (REMOVES FORMALDAHYDE) 

42 IF (NINT (TRTTYP) . EQ. 2 ) THEN 

43 FACID=1 

44 FCI=0 

45 ELSEIF (NINT (TRTTYP) . EQ. 1) THEN 

46 FACID=0 

47 FCI=1 

48 ELSE 

49 FACID=0 

50 FCI=0 

51 ENDIF 

52 C 

53 C TEST FOR NO BED FLOW ( BEDQ=<0 ) OR TN-TN1<=0 ; BEDL , BEDDIA, DENCH=0 

54 IF ( (BEDQ.LE.0) .OR. (TN-TN1 . LE . 0 ) .OR. (CARTL.LE.0) .OR. (BEDOD.LE.0) 

55 +. OR. (DENCH. LE.0) ) THEN 

56 EFF=0 

57 GOTO 199 
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58 ENDIF 

59 C TEST FOR Cl CHARCOAL AND FORMALDEHYDE {FCI=1 AND MW=30.03) 

60 IF( (MW.EQ.30.03) .AND. (FCI.EQ.l) ) THEN 

61 CALL RCICH ( EFF , EMAX , CARTL , BEDOD, BEDID, DENCH , SMR, BEDQ) 

62 GOTO 199 

63 ENDIF 

64 C 

65 C TEST FOR AMMONIA AND PHOS ACID ON CHAR(FACID=1 AND MW=17 . 0 ) 

66 IF ( (MW.EQ.17.0) .AND. (FACID.EQ.l) ) THEN 

67 CALL RACCH(EFF, EMAX, CARTL, BEDOD, BEDID, DENCH, SMR) 

68 GOTO 199 

69 ENDIF 

70 C 

71 C TEST FOR MOL VOL=0 (NO CHAR REMOVAL) 

72 IF (VMOL.EQ.O) THEN 

73 EFF=0 

74 GOTO 199 

75 ENDIF 

76 C 

77 C CHARCOAL REMOVAL EFFICIENCY CALCULATION 

78 C BED LENGTH (M) -ASSUMES THIN BED 

79 BEDL= (BEDOD-BEDID) / 2 

80 IF (BEDL.LT.O) BEDL=0 

81 C BED WGT (KG) 

82 BEDWGT=DENCH * . 785* (BEDOD**2-BEDID**2 ) * CARTL 

83 C SUPERFICIAL BED VEL ( FT/MIN) 

84 BEDVEL=BEDQ* .0348/ ( (BEDOD+BEDID) * CARTL) 

85 C TEST FOR CIN TOO SMALL IN AVAL CALC 

86 IF (CIN.LT. IE-20) CIN=lE-20 

87 AVAL= (TCABIN/VMOL) *LOG10 (VCONC/CIN) 

88 C ADS ZONE LENGTH FOR 90% REMOVAL (M) 

89 LADS=AVAL* .000275* (BEDVEL/ 1 . 3 ) * * . 8 

90 C GET QI(CC LIQ CONT/GM CHAR) 

91 CALL FQI (AVAL, QI,FACID, SOL, RH) 

92 C LENGTH OF BED PREVIOUSLY USED BY CONT AT THIS C INLET (M) 

93 LPREV=SMR*1 . 0E-6*COEXIS*BEDL/ ( DCONT*BEDWGT*QI ) 

94 C RATE OF BED USAGE (M BED/ MG CONT) 

95 LIMM=1.0E-6*COEXIS*BEDL/ (DCONT*BEDWGT*QI ) 

96 C LENGTH OF BED AVAILABLE FOR ADS ZONE AT BEGINNING OF INCR (M) 

97 LAVNl=BEDL-LPREV 

98 IF (LAVN1.LT.0) LAVN1=0 

99 C FIX HERE IF DESORPTION IS DESIRED 

100 IF (LAVNl /LADS . GT .20) THEN 

101 EFFAV=EMAX 

102 ELSE 

103 C INIT INCR EFF BASED ON C IN AND BED L AVAIL AT BEG OF INCR (DEC) 

104 EFAVNl=EMAX* ( 1-EXP ( - 2 . 302 5 851 *LAVNl /LADS ) ) 

105 C LOOP FOR EFFICIENCY 

106 EFFAV=EFAVNl 

107 DO 399 J=l, 10, 1 

108 C LENGTH OF BED UTILIZ IN INCR (M) 

109 LUT I L =C IN * BEDQ * EF FAV * (TN-TN1 ) *LIMM 

110 IF (LUTIL.GT. LAVNl) THEN 

111 GOTO 299 

112 ELSE 

113 C AVERAGE BED LENGTH AVAIL (M) 

114 LAVAV = LA VN 1 - LUT I L / 2 

115 IF ( (LAVAV/LADS) .GE. 20) THEN 

116 EFFAV=EMAX 

117 GOTO 299 
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118 ELSE 

119 C AV EFF BASED ON AV BED L AVAIL (DEC) 

120 EFFAV=EMAX* (1-EXP (-2 . 302 5 85 1 *LAVAV/LADS ) ) 

121 ENDIF 

122 ENDIF 

123 399 CONTINUE 

124 299 ENDIF 

125 C MAX EFF BASED ON C IN AND RATE OF BED USAGE (DEC) 

126 EFFMAX=LAVNl/ (CIN*BEDQ* (TN-TNl ) *LIMM) 

127 IF (EFFAV.GT. EFFMAX) EFFAV=EFFMAX 

128 IF (EFFAV.LT. 0) EFFAV=0 

129 IF (EFFAV.GT. EMAX) EFFAV=EMAX 

130 C EFF=ACTUAL EFF OUTPUT FROM SUBROUTINE 

131 EFF=EFFAV 

132 C REMOVE THIS CHECK IF DESORPTION IS ADDED 

133 199 IF (EFF.LT.O) EFF=0 

134 IF (EFF. GT. EMAX) EFF=EMAX 

135 RETURN 

136 END 

NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 

C ***************************************************** 

138 C * SUBROUTINE RACCH - CALCULATES REMOVAL EFF * 

139 C * BED WITH NH3 AND 1.22 MILLIMOLE H3P04 ON CHAR * 

140 C ***************************************************** 

141 SUBROUTINE RACCH ( EFF, EMAX , CARTL , BEDOD, BEDID, DENCH , SMR) 

142 C OUTPUTS 

14 3 C EFF =OUTPUT REMOVAL EFF (DEC) 

144 C INPUTS 

145 C EMAX=MAXIMUM BED REMOVAL EFF (DEC) 

146 C CARTL =CARTRIDGE LENGTH (M) 

147 C BEDOD=BED OUTSIDE DIAMETER (M) 

148 C BEDID=BED INSIDE DIAMETER (M) 

149 C DENCH =CHARCOAL DENSITY (KG/CU M) 

150 C SMR=SUM OF MASS OF CONT REMOVED AT BEG OF INCR (MG) 

151 C 

152 C FOR AMMONIA CAPACITY AT SMAC 

153 C CHAR USED (KG) 

154 CHRUSD=1 . 6E-4*SMR 

155 C CHAR BED WGT(KG) 

156 BEDWGT=CARTL * (BEDOD* * 2 - BEDID* *2 ) * . 785*DENCH 

157 IF (CHRUSD. LT. 0 . 8*BEDWGT) THEN 

158 EFF=EMAX 

159 ELSE 

160 EFF= EMAX* SIN ( ( BEDWGT- CHRUSD) *1.57/ (BEDWGT*0.2) ) 

161 ENDIF 

162 C PREVENTS NEGATIVE EFF FOR REACTION 

163 C IF (EFF.LT.O) EFF=0 

164 I F ( EFF . GT . EMAX ) EFF=EMAX 

165 RETURN 

166 END 

NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
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167 C 

168 C 

169 C 

170 C 

171 

172 C 

173 C 

174 C 

175 C 

176 C 

177 C 

178 C 

179 C 

180 C 

181 C 

182 C 

183 

184 C 

185 

186 

187 

188 

189 

190 

191 C 

192 

193 

194 

195 

196 C 

197 C 

198 

199 

200 

NUMBER 

NUMBER 

NUMBER 

NUMBER 


************************* 


************************** 


* * 


★ 

★ 

* 


SUBROUTINE RCICH - CALCULATES REMOVAL EFF 
FOR FORMALDELYDE AND Cl CHAR BED 

★ ★★★★★★★it********************************** 


*★★★**** 


★ 


* 


SUBROUTINE RC ICH ( EFF , EMAX , CARTL , BEDOD , BEDID , DENCH , SMR , BEDQ } 
OUTPUTS 

EFF=OUTPUT REMOVAL EFF (DEC) 

INPUTS 

EMAX=MAXIMUM BED REMOVAL EFF (DEC) 

CARTL =CARTRIDGE LENGTH (M) 

BEDOD=BED OUTSIDE DIAMETER (M) 

BEDID=BED INSIDE DIAMETER (M) 

DENCH =CHARCOAL DENSITY (KG/CU M) 

SMR = SUM OF MASS OF CONT REMOVED AT BEG OF INCR (MG) 


BEDQ=BED FLOW RATE (CU M/HR) 


BEDWGT=CARTL* (BEDOD* * 2 -BEDID* *2 ) * . 7 85 *DENCH 
PERCENT OF BED WEIGHT CONSUMED (DEC) 

PBWGT=SMR/ (BEDWGT*1E6) 

IF (PBWGT. LT. . 0012 ) THEN 
EFF=1- PBWGT* 83 . 3 
ELSE 

EFF= . 9*COS ( PBWGT* 1 . 57/ . 05) 

END IF 

BED RESIDENCE TIME (SEC) 

BREST= (BEDOD- BEDID) * CARTL* ( BEDOD+BEDID) *3600/ (BEDQ* 1 . 273 ) 
I F ( BREST . LT . 0 . 2 5 ) TH EN 
EFF=EFF* BREST/ . 2 5 
ENDIF 

PREVENTS NEGATIVE EFF FOR REACTION 
IF (EFF.LT.0) EFF=0 
IF (EFF. GT. EMAX) EFF=EMAX 
RETURN 
END 


OF WARNINGS IN PROGRAM UNIT: 0 
OF ERRORS IN PROGRAM UNIT: 0 
OF WARNINGS IN COMPILATION : 0 
OF ERRORS IN COMPILATION : 0 
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RM/ FORTRAN Compiler (V2.42) 

Source File: C:\RMFORT\TCC\REGEN.F Options: /C 80 /L /BY 05/21/92 13:00:29 

1 c **************************************************************** 

2 C * REGENERATION SUBROUTINE -REGEN * 

3 C * FOR ALL BEDS (3-15) DETERMINES IF BED IS TO BE REGENERATED * 


4 C * AT BEGINNING OF TIME INCREMENT, AND IF REGENERATION IS TO * 

5 C * DURING THE ENTIRE TIME INCREMENT- IF THE BED IS TO BE * 

6 C * REGENERATED THE MASSES STORED ARE SET TO ZERO, AND IF * 

7 C * REGENERATION IS TO OCCUR THROUGHOUT THE TIME INCREMENT THE * 

8 C * BED FLOW RATE IS SET TO ZERO; OTHERWISE IT IS SET TO THE * 

9 C * ORIGINAL VALUE * 


3_0 Q *************★*★*********★★+*★★***★****★★************★**★★***★★* 

11 SUBROUTINE REGEN (TN, TNI , DD, NROW, NCOL, CC , NROWl , NCOLl , 

12 +CDI , NROW2 , NCOL 2 , LIN, LIN2 , IMSGDN) 

13 INTEGER NROW, NCOL , NROWl , NCOLl , NROW2 , NCOL2 

14 REAL DD (NROW, NCOL) 

15 REAL CC (NROWl , NCOLl ) 

16 REAL CDI (NROW2 , NCOL2 ) 

17 

18 C NOTE: BEFORE RUNNING THIS SUBROUTINE THE ORIGINAL FLOW RATES FROM 19 C 

20 C TIME INCREMENT THEY MUST BE RESTORED TO COL 2 

21 C DIRECT INPUTS: 

22 C TN= INCREMENT END TIME (HRS); TNl=INCR BEGINNING TIME HRS 

23 C DD, NROW, NCOL=NAME AND SIZE OF MAT DD 

24 C CC, NROWl, NCOLl =NAME AND SIZE OF MAT CC 

25 C CDI , NROW2 , NCOL2=NAME AND SIZE OF MAT CDI 

26 C LIN=NO . OF CONT IN MAT CDI 

27 C LIN2 =NO . OF DEVICES IN MAT DD 

28 

29 C OTHER INPUTS FROM MAT DD 

30 C TIR=INITIAL (FIRST) REGENERATION TIME (HRS) 

31 C TRC I =REGEN /CHANGEOUT INTERVAL (HRS) 

32 C TRD= REGENERATION DURATION (HRS) 

33 C DEVICE NO. , TYPE, FLOW RATE, ETC 

34 

35 C OUTPUTS: 

36 C A) IF REGENERATION OCCURS AT THE BEGINNING OF ANY TIME INCREMENT 

37 C 1 ) FOR ANY DEVICE WHICH IS A CHARCOAL BED 

38 C FOR ALL CONT 1 TO LIN IT PUTS SUM MASS REM=0 IN MAT CC 

39 C COL 12 , 15, 18 . . . . 48 AS APPROPRIATE FOR THAT DEVICE 

40 C 2) FOR ANY DEVICE WHICH IS A LIOH BED 

41 C IT DOES 1) ABOVE, AND IN ADDITION PUTS SUM MASS REM=0 IN 

42 C MAT DD COL 16 FOR THAT DEVICE 

43 C B) IF REGENERATION IS OCCURRING THROUGHOUT THE WHOLE INTERVAL 

44 C IT SETS Q OF DEVICE=0; IF REGENERATION IS NOT OCCURRING, IT 

45 C SETS Q=THE ORIGINAL VALUE 

46 C 

47 C SUBROUTINES REQUIRED: 

48 C REGCHG 

49 C 

50 C START LOOP FOR ALL DEVICES 3 TO 15 

51 DO 100 J=3 , LIN2 

52 C IF DEVICE DOES NOT EQUAL CHARCOAL OR LIOH THEN GO TO END OF LOOP 

53 I F ( DD ( J , 3 ) . NE . 3 . AND . DD ( J , 3 ) . NE . 4 . AND . DD ( J , 3 ) . NE . 5 ) GOTO 100 

54 

55 C ASSIGN PROPER VARIABLES FOR DEVICE 

56 C DEVICE = CHARCOAL 

57 I F ( DD ( J , 3 ) . EQ . 3 . OR . DD ( J , 3 ) . EQ . 4 ) THEN 

58 C REGENERATION/ CHANGEOUT INTERVAL (HRS) 

59 TRC I =DD ( J , 15) 
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60 C 

61 

62 C 

63 

64 

65 

66 C 

67 

68 C 

69 

70 C 

71 

72 C 

73 

74 

75 

76 C 

77 C 

78 

79 C 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 C 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 C 

104 

105 

106 

107 

108 C 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 
119 


REGENERATION DURATION (HRS) 

TRD=DD ( J, 16) 

INITIAL (FIRST) REGENERATION 
TIR=DD ( J, 14) 

END IF 

DEVICE = LIOH 

IF (DD ( J, 3 ) . EQ. 5) THEN 

REGENERATION/CHANGEOUT INTERVAL (HRS) 

TRCI=DD ( J , 14 ) 

REGENERATION DURATION (HRS) 

TRD=0 

INITIAL (FIRST) REGENERATION 
TIR=DD ( J , 13) 

END IF 

CHECK AND FIX INPUT AS REQ + PRINT WARNINGS 
TIME INCREMENT (HRS) 

TINC=DD (1,11) 

INITIAL TIME NOT EQUAL TO MULTIPLE OF TIME INCREMENT 
IF ( AINT (TIR/TINC ) .NE. (TIR/TINC) ) THEN 
TIR=AINT (TIR/TINC) *TINC 
OPEN ( IMSGDN, FILE=' CON' , IOSTAT=IOVAL) 

WRITE ( IMSGDN, *) 1 INCREMENT BEGINNING TIME ' ,TN1, 

+ ' DEV NO. ' , (J) _ mTWT , 

WRITE ( IMSGDN, * ) 'INITIAL TIME NOT = MULTIPLE OF TIME 

+ INCREMENT’ 

WRITE ( IMSGDN , * ) ' TRUNCATED TO ' , TIR 
CLOSE (IMSGDN) 

END IF 

REGEN/CHGOUT INTERVAL < TIME INCR OR NOT= MULTIPLE OF TIME INCR 
IF (AINT (TRCI/TINC) .NE. (TRCI/TINC) ) THEN 
TRCI=AINT (TRCI/TINC) *TINC 
OPEN (IMSGDN, FILE='CON' , IOSTAT=IOVAL) 

WRITE ( IMSGDN, *)' INCREMENT BEGINNING TIME ',TNl, 

+ ' DEV NO. ' , (J) ^ 

WRITE (IMSGDN, *) 'REGEN/CHG TIME NOT = MULTIPLE OF TIME 

+ INCREMENT' 

WRITE { IMSGDN, * ) ' TRUNCATED TO ' , TRCI 
CLOSE (IMSGDN) 

END IF 

IF REGEN/CHGOUT INTERVAL <= O THEN GOTO END OF LOOP FOR CONT 
IF (TRCI.LE.0) THEN 
GOTO 100 

END IF 

REGEN DURATION < OR NOT = MULTIPLE OF TIME INCREMENT 
IF (AINT (TRD/TINC) .NE. (TRD/TINC ) ) THEN 
TRD=AINT (TRD/TINC ) *TINC 
OPEN ( IMSGDN , FILE= ' CON ' , IOSTAT=IOVAL) 

WR ITE( IMSGDN, *)' INCREMENT BEGINNING TIME ',TNl, 

+ ' DEV NO . ' , ( J ) 

WRITE ( IMSGDN, * ) 'REG DURATION NOT = MULTIPLE OF TIME 
+ INCREMENT' 

WRITE ( IMSGDN , * ) ' TRUNCATED TO ' , TRD 
CLOSE (IMSGDN) 

END IF 
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I C REGENERATION DURATION > REGEN/CHG INTERVAL 
IF(TRD.GT.TRCI) THEN 
TRD=TRCI 

OPEN ( IMSGDN, FILE= 'CON' , IOSTAT=IOVAL ) 

WRITE { IMSGDN, *)' INCREMENT BEGINNING TIME ' , TNI 

i + 'DEV NO. ' , (J) 

WRITE (IMSGDN, *) 'REGEN DURATION > REGEN/CHG INTERVAL' 

WRITE ( IMSGDN, * ) 'TRUNCATED TO 1 , TRD 
CLOSE (IMSGDN) 

END IF 

L CHECK TO SEE IF REGENERATION OCCURS AT BEGINNING OF TIME INCR, 
C AND IF REGEN OCCURS THROUGHOUT WHOLE TIME INCREMENT 

CALL REGCHG (TNI , TRCI , TRD, TIR, TINC , IRBFLG, IRTFLG) 

C REGENERATION OCCURS AT BEGINNING OF INCREMENT 
IF (IRBFLG. EQ. 1) THEN 

C PUT SUM MASS REM =0 IN MAT CC FOR THIS DEVICE 

C START LOOP FOR ALL CONT FOR THIS DEVICE 

K=J* 3+3 

DO 101 1=1, LIN 
CC ( I , K ) =0 

101 CONTINUE 

C IF DEVICE = LIOH BED PUT SUM MASS=0 IN DD(J,16) 

IF ( DD ( J , 3 ) . EQ . 5 ) THEN 
DD ( J, 16) =0 
END IF 
END IF 

C REGENERATION OCCURS THROUGHOUT ENTIRE INCREMENT 
IF (IRTFLG. EQ. 1) THEN 
C SET DEVICE Q=0 

DD ( J , 2 ) =0 
ELSE 

C SET DEVICE Q= ORIGINAL VALUE 

DD(J, 2) =DD(J, 7) 

END IF 

C END OF J LOOP FOR EACH DEVICE 
100 CONTINUE 

RETURN 

END 

!R OF WARNINGS IN PROGRAM UNIT: 0 
2R OF ERRORS IN PROGRAM UNIT: 0 


* AUXILIARY REGENERATION SUBROUTINE -REGCHG * 

* DETERMINES IF REGEN /CHANGEOUT IS TO OCCUR AT BEGINNING OF * 

* TIME INCREMENT- ALSO DETERMINES IF REGENERATION IS OCCURRING * 

* THROUGHOUT THE TIME INCREMENT * 

SUBROUTINE REGCHG ( TNI , TRCI , TRD , TIR , TINC , IRBFLG, IRTFLG) 

INPUTS : 

TNI = INCREMENT INITIAL TIME (HRS) 
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176 

177 

178 

179 

180 
181 
182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

218 
219 


TRC I = C HANG EOUT / R EG ENERAT I ON INTERVAL (HRS) 

TRD= REGENERATION DURATION (HRS) 

TIR=INITIAL (FIRST) REGENERATION TIME (HRS) 

TINC=TIME INCREMENT (HRS) 

OUTPUTS : 

REGENERATION OCCURS AT BEGINNING OF TIME INCREMENT (Y OR 1 
(IRBFLG=1 FOR Y & 0 FOR N) 

REGENERATION IS OCCURRING THROUGHOUT THE WHOLE INCREMENT (Y 
(IRTFLG=1 FOR Y & 0 FOR N) 

C REGENERATION OCCURS AT BEGINNING OF TIME INCREMENT 

IF (TN1.EQ.0) GOTO 10 
IF(TNl.LT.TIR) GOTO 10 
IF (TRC I . LE . 0 ) GOTO 10 

IF ( AINT ( (TNl-TIR) /TRCI) .EQ. ( (TNl-TIR) /TRCI) ) THEN 
C REGENERATION OCCURS 

IRBFLG=1 
GO TO 20 
END IF 

C NO REGENERATION OCCURS 

010 IRBFLG=0 
020 CONTINUE 

C REGENERATION OCCURRING THROUGHOUT ENTIRE TIME INCREMENT 
IF (TRCI .LE. 0) GOTO 30 

IF( (TRD.LE. 0) .OR. (TNI .LT.TIR) ) GOTO 30 
IF (TNI .GE. AINT ( (TNl-TIR) /TRCI ) *TRCI+TIR+TRD) THEN 
GOTO 30 
ELSE 

C REGENERATION OCCURS 

IRTFLG=1 
GOTO 40 
END IF 

C REGEN DOESN'T OCCUR 

030 IRTFLG=0 

040 CONTINUE 
RETURN 
END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 


) 

OR N) 
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RM/FORTRAN Compiler (V2.42) 

Source File: C:\RMFORT\TCC\RINCDD. Options: /C 80 /L /BY 05/21/92 13:00:48 
£ ★ ★*******★★★★*★**★★★★**★★★*★★★*★★★★★★**★★★★★★**★★*★★**★★★*★★*★*★★ 

2 C * SUBROUTINE RINCDD * 

3 C * SUBROUTINE TO OPERATE ON INCREMENT DEPENDENT DATA * 

4 C * READS DATA FROM MAT TT AND PUT IT IN THE PROPER PLACES IN * 

5 C * MAT CDI OR MAT DD - USED AT THE BEGINNING OF EACH TIME INCR * 

c. : p ***************************************************************** 


7 

8 SUBROUTINE RINCDD ( I , TN, TNI , DD , NROW, NCOL , LIN2 , 

9 +CC , NROW1 , NCOL1 , CDI , NROW2 , NCOL2 , LIN , TT , NTTROW , NTTCOL , LINl ) 

1 0 INTEGER NROW , NCOL , NROWl , NCOLl , NROW2 , NCOL2 , NTTROW , NTTCOL 

11 REAL DD (NROW, NCOL) 

1 2 REAL CC ( NROWl , NCOL 1 ) 

13 REAL CDI (NROW2 , NCOL2 ) 

14 REAL TT (NTTROW, NTTCOL) 

15 

16 C SUBROUTINES REQUIRED: 

1 7 C NONE 

18 C DD , NROW , NCOL , LIN2 =NAME , DIM & NO DEV IN MAT DD 

19 C CC, NROWl, NCOLl=NAME & DIM OF MAT CC 

20 C CDI, NROW2, NCOL, LIN=NAME, DIM & NO CONT IN MAT CDI 

21 C TT, NTTROW, NTTCOL, LIN1=NAME, DIM & NO ITEMS IN MAT TT 

22 

23 IF (LIN1.EQ.0) GOTO 999 

24 C BEGIN LOOP FOR ALL LINES IN MAT TT 

25 DO 100 K=1 , LINl 

26 

27 C CHECK FOR TIME >= TNI AND < TN 

28 I F ( ( TT ( K , 1 ) . LT . TN 1 ) . OR . ( TT ( K , 1 ) . GE . TN ) ) THEN 

29 GO TO 100 

30 END IF 

31 

32 C IDENTIFY VARIABLES 

33 ICONTN=NINT(TT(K, 2) ) 

34 GENRT=TT(K, 3) 

35 IDEVNO=NINT(TT(K, 4) ) 

36 DEVQ=TT (K , 5 ) 

37 ICOLNO=TT(K, 6) 

38 VAL=TT ( K , 7 ) 

39 

40 C ICONTN= INTEGER CONTAMINANT NO. - TT(K,2) 

41 C GENRT=CONT GENERATION RATE (MG/HR) - TT(K,3) 

42 C I DEVNO= INTEGER DEVICE NUMBER - TT(K,4) 

43 C DEVQ=DEVICE FLOW RATE (CU M/HR) - TT(K,5) 

44 C I COLNO= INTEGER COLUMN NUMBER IN MAT DD - TT(K,6) 

45 C VAL=NEW VALUE IN MAT DD - TT(K,7) 

46 

47 C CASE NO. 1 - CHANGE CONTAMINANT GENERATION RATE 

48 

49 IF ( ( ICONTN. GT . 0 ) . AND. ( ICONTN . LE . LIN) ) THEN 

50 IF ( IDEVNO. EQ. 1 ) THEN 

51 CDI (ICONTN, 1) =GENRT 

52 ELSEIF ( (IDEVNO. GE. 3) .AND. (IDEVNO. LE.LIN2) ) THEN 

53 CDI (ICONTN, (7+IDEVNO) ) =GENRT 

54 ENDIF 

55 GOTO 100 

56 ENDIF 

57 

58 C CASE 2 - CHANGE DEVICE FLOW OR OTHER DD DATA 
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59 C THIS CASE WORKS ONLY IF ANY CONT NO. <=0 

60 C MUST USE -1 FOR ANY Q OR NEW VALUE NOT TO BE CHANGED 

61 C MAT DD COL NO. <0 ALSO STOPS NEW VALUE FROM BEING CHANGED 

62 

63 IF ( ICONTN.LE. 0} THEN 

64 IF ( ( IDEVNO . GE . 1 ) . AND . ( IDEVNO . LE . LIN2 ) ) THEN 

65 IF(DEVQ.GE.O) THEN 

66 C CHANGE DEVICE FLOW IN MAT DD 

67 DD (IDEVNO, 2) =DEVQ 

68 END IF 

69 IF ( ( ICOLNO . GE . 1 ) . AND. ( ICOLNO. LE . 16 ) ) THEN 

70 IF (VAL.GE.O) THEN 

71 C CHANGE VALUE IN MAT DD 

72 DD ( IDEVNO, ICOLNO) =VAL 

73 END IF 

74 END IF 

75 END IF 

76 END IF 

77 100 CONTINUE 

78 999 RETURN 

79 END 

NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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RM/FORTRAN Compiler (V2.42) 

Source File: C:\RMFORT\TCC\RRIN.FO Options: /C 80 /L /BY 


05/21/92 13:00:57 


\ IXi'J L Ui\ X \ \iuvj. iN.i u ^ ***-'• / ' — ' 

^ £ ****************************************************************** 

2 C * SUBROUTINE RRIN * 

3 C * SUBROUTINE TO READ REAL DATA INTO MAT XX (ROW, COL) 

4 C * RETURNS NUMBER OF LINES OF DATA READ FROM FILE * 

5 C * READS FROM COL 1 TO COL LSTCOL 

6 c ****************************************************************** 

7 C NOTE: INPUT NUMBERS MUST BE SEPARATED BY BLANKS 

8 SUBROUTINE RRIN (XX, NROW,NCOL, LSTCOL, LIN) 

9 INTEGER NROW, NCOL , IOVAL , LSTCOL , LIN 

10 CHARACTER FNAME*24 

11 REAL XX (NROW, NCOL) 

12 IF (LSTCOL. GT. NCOL) LSTCOL=NCOL 

13 010 READ ( * , 1 (A) ’ ) FNAME 

14 OPEN ( 1 , FILE=FNAME , STATUS= 1 OLD 1 , IOSTAT=IOVAL ) 

15 IF ( IOVAL. NE.0) GOTO 900 

16 LIN=0 

17 DO 100 1=1, NROW 

18 READ ( 1 , * , IOSTAT= IOVAL , END=500 , ERR=900 ) (XX (I , J) , J=1 , LSTCOL) 

19 LIN=LIN+1 

20 100 CONTINUE 

21 500 WRITE (*, ’ (A) ’ ) ‘ DONE WITH FILE INPUT’ 

22 WRITE (*,*) ’ ' 

23 CLOSE (1) 

24 GOTO 990 

25 900 WRITE (*,*)’ IOERROR= ’, IOVAL 

26 CLOSE (1) 

27 WRITE (*,*) ’WHAT IS THE INPUT FILE NAME? ’ 

28 GOTO 10 

29 990 RETURN 

30 c ************ END OF SUBROUTINE RRIN ****************************** 

31 END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 

NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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RM/FORTRAN Compiler (V2.42) 

Source File: C:\RMFORT\TCC\RROUT.F Options: /C 80 /L /BY 05/21/92 13:01:03 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 


★ ★★*★★★★★★*★★****★★*****★★★★★**★**★*****★★**★★★★★★****★*★******★** 


c 

c 

c 

c 

c 


SUBROUTINE RROUT 

SUBROUTINE TO WRITE DATA TO CONSOLE, OR PRINTER 
WRITES REAL DATA FROM MAT XX (ROW, COL) 

WRITES FROM FSTCOL TO LSTCOL 


****************************************************************** 

SUBROUTINE RROUT (XX , NROW, NCOL , FSTCOL, LSTCOL , LIN , IMSGDN, FNAME , 

+ IDEVNO, IOVAL) 

INTEGER NROW, NCOL, IOVAL, FSTCOL, LSTCOL, LIN, IDEVNO 
CHARACTER FNAME *2 4 , DES* 1 
REAL XX (NROW, NCOL) 

IF ( FSTCOL. GT. NCOL) FSTCOL=NCOL 
IF (LSTCOL. GT. NCOL) LSTCOL=NCOL 
IF (FSTCOL. GT. LSTCOL) FSTCOL=LSTCOL 
010 OPEN ( IMSGDN, FILE=* CON' , IOSTAT=IOVAL) 

WRITE (IMSGDN, ' (A) * ) ' WRITE TO LPTl OR CON OR END 1 

CLOSE (IMSGDN) 

READ ( * , 1 (A) 1 ) FNAME 

QUIT IF FNAME = END 

IF (FNAME. EQ. ' END' ) GOTO 990 

IF ( (FNAME. NE. 'LPTl ') .AND. (FNAME. NE. 'CON' ) ) GOTO 10 
OPEN ( 1 , FILE=FNAME , IOSTAT=IOVAL) 

IF ( IOVAL. NE.0) GOTO 900 
DO 110 1=1, LIN 

WRITE (1,70, IOSTAT= IOVAL , ERR= 900) { XX ( I , J ) , J=FSTCOL , LSTCOL ) 

WRITE (IDEVNO, 70, IOSTAT= IOVAL, ERR=9 00 ) (XX (I, J) , J=FSTCOL, LSTCOL) 
FORMAT ( IX, 7G1 1.4) 

WRITE (1, *, IOSTAT=IOVAL , ERR=9 00 ) 

WRITE ( IDEVNO, * , IOSTAT=IOVAL , ERR=900 ) 

110 CONTINUE 

CLOSE (1) 

GOTO 990 

900 OPEN( IMSGDN, FILE= 'CON' , IOSTAT=IOVAL) 

WRITE ( IMSGDN, * ) ' IOERROR= 1 , IOVAL 
CLOSE (IMSGDN) 

CLOSE (1) 

CLOSE (IDEVNO) 

GOTO 10 


070 


990 RETURN 


*************** 


END OF SUBROUTINE RROUT 


************************ 


END 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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RM/ FORTRAN Compiler (V2.42) 

Source File: C:\RMFORT\TCC\RROUT2. Options: /C 80 /L /BY 05/21/92 13:01:08 

1 c *********************★*****★★★★★★*****★*★*★★****★★**★*★*★*★******* 

2 C * SUBROUTINE RROUT2 * 

3 C * SUBROUTINE TO WRITE DATA TO CONSOLE, OR PRINTER * 

4 C * WRITES REAL DATA FROM MAT XX (ROW, COL) * 

5 C * WRITES FROM FSTCOL TO LSTCOL * 

6 C ****************************************************************** 

7 SUBROUTINE RROUT2 (XX, NROW, NCOL, FSTCOL, LSTCOL, LIN, IMSGDN) 

8 INTEGER NROW, NCOL , IOVAL , FSTCOL, LSTCOL , LIN , IDEVNO 

9 CHARACTER FNAME *24, DES * 1 

10 REAL XX (NROW, NCOL) 

11 IF ( FSTCOL. GT. NCOL) FSTCOL=NCOL 

12 IF (LSTCOL. GT. NCOL) LSTCOL=NCOL 

13 IF ( FSTCOL. GT. LSTCOL) FSTCOL=LSTCOL 

14 010 OPEN ( IMSGDN, FILE- ' CON ' ,.IOSTAT=IOVAL ) 

15 WRITE (IMSGDN, ' (A) ' ) ' WRITE TO LPT1 OR CON OR END ' 

16 CLOSE (IMSGDN) 

17 READ ( * , ' (A) ' ) FNAME 

18 QUIT IF FNAME = END 

19 IF (FNAME. EQ. ' END' ) GOTO 990 

20 IF ( (FNAME .NE. 1 LPTl ') .AND. (FNAME. NE. 'CON 1 ) ) GOTO 10 

21 OPEN ( 1 , FILE=FNAME, IOSTAT=IOVAL ) 

22 IF ( IOVAL. NE. 0 ) GOTO 900 

23 DO 110 1=1, LIN 

24 WRITE (1, 70, IOSTAT=IOVAL,ERR=900) (XX ( I , J) , J=FSTCOL , LSTCOL) 

25 070 FORMAT (IX, 7G11 . 4 ) 

26 WRITE (1, *, IOSTAT= IOVAL , ERR=900 ) 

27 110 CONTINUE 

28 CLOSE (1) 

29 GOTO 990 

30 900 OPEN ( IMSGDN, FILE= ' CON ' , IOSTAT=IOVAL ) 

31 WRITE (IMSGDN, * ) ' IOERROR= ' , IOVAL 

32 CLOSE (IMSGDN) 

33 CLOSE (1) 

3 4 GOTO 1 0 

35 990 RETURN 

36 C *************** END Q F SUBROUTINE RROUT ************************ 

3 7 END 

NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
NUMBER OF ERRORS IN COMPILATION : 0 
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RM/ FORTRAN Compiler (V2.42) 

Source File: C:\RMFORT\TCC\SLIOH.F Options: /C 80 /L /BY 05/21/92 13:01:15 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 


Q *************************************************************** 

C * SUBROUTINE -SLIOH 

C * SUM LIOH USED IN TIME INCREMENT FOR EACH BED ONE AT A TIME 

C * AND FOR ALL CONTAMINANTS FOR EACH BED 

Q *************************************************************** 

SUBROUTINE SL IOH ( TN , TNI , DD , NROW , NCOL , CC , NROWl , NCOLl , 

+CDI , NROW2 , NCOL 2 , LIN, LIN2 ) 

INTEGER NROW, NCOL , NROWl , NCOLl , NROW2 , NCOL2 
REAL DD (NROW, NCOL) 

REAL CC (NROWl, NCOLl) 

REAL CD I (NROW 2 , NCOL2 ) 

C 

C SUBROUTINES REQUIRED: NONE 
C DIRECT INPUTS: 

C TN= INCREMENT END TIME ( HRS ) ; TNl=INCR BEGINNING TIME HRS 
C DD , NROW , NCOL=NAME AND SIZE OF MAT DD 

C CC, NROWl , NCOLl =NAME AND SIZE OF MAT CC 
C CD I , NROW2 , NCOL2=NAME AND SIZE OF MAT CDI 
C LIN=NO. OF CONT IN MAT CDI 
C LIN2 =NO . OF DEVICES IN MAT DD 

C OTHER INPUTS FROM MAT DD 
C DD (J, 3 ) =DEVICE NUMBER 

C DD ( J , 16) = AMT OF LIOH PREVIOUSLY USED BY DEVICE 

C CDI (1,17) =LB LIOH UTIL/LB CONT ADSORBED IN BED (FOR ONE CONT) 

C OUTPUTS (STORED IN MAT DD) : 

C DD ( J , 16) = AMOUNT OF LIOH UTILIZED BY DEVICE THROUGH THE END OF 

C THIS TIME INCREMENT 

C DD ( J, 15) =RATE OF LIOH USAGE FOR DEVICE 
C 

K=1 1 

C START LOOP FOR ALL DEVICES 3 TO 15 

DO 100 J=3 , LIN2 

C CHECK FOR DEVICE = LIOH BED 

IF ( DD ( J , 3 ) . EQ . 5 ) THEN 
C RATE OF LIOH UTILIZATION (KG/HR) 

RWUTLI=0 

C BEGIN LOOP FOR ALL CONTAMINANTS 

DO 110 1=1, LIN 

RWUTLI=RWUTLI+CC (I , K) *CDI (1,7) *lE-6 
110 CONTINUE 

C STORE RATE OF LIOH UTILIZATION IN MAT DD FOR THIS DEVICE 

DD ( J , 15 ) =RWUTLI 

C UPDATE AMOUNT OF LIOH UTIL THROUGH THE END OF TIME INCR(KG) 

DD ( J, 16) =DD(J, 16) +RWUTLI* (TN-TN1 ) 

END IF 
K=K + 3 

C END J LOOP 
100 CONTINUE 
RETURN 
END 


★ 


★ 


NUMBER OF WARNINGS IN PROGRAM UNIT: 0 
NUMBER OF ERRORS IN PROGRAM UNIT: 0 
NUMBER OF WARNINGS IN COMPILATION : 0 
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NUMBER OF ERRORS 


IN COMPILATION : 0 
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APPENDIX B 

TOXIC HAZARD INDEX DESCRIPTION 




V 


.MS 
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The toxic hazard index, or T-value, is the method used by toxicologists to assess the acceptability 
of an atmosphere containing a mixture of contaminants. This approach is derived from the American 
Conference of Governmental Industrial Hygienists guidelines for setting threshold limit values for con- 
taminant mixtures. Since the effects on humans of many atmospheric contaminants are considered to be 
additive, this mixture approach is applied to 16 contaminant groups. The groups considered in the T- 
value calculation used in the TCCS computer program are the following: 

1. Alcohols 

2. Aldehydes 

3. Aromatic hydrocarbons 

4. Esters 

5. Ethers 

6. Chlorocarbons 

7. Chlorofluorocarbons 

8. Fluorocarbons 

9. Hydrocarbons 

10. Inorganic acids 

1 1 . Ketones 

12. Mercaptans and sulfides 

13. Nitrogen oxides 

14. Organic acids 

15. Organic nitrogens 

16. Miscellaneous 

The group numbers used in the computer program output correspond to the above group listing. 

The T-value is calculated for each group by calculating the sum of the ratios of the contami- 
nants’ concentrations to their maximum allowable concentration, while the overall T-value is the sum of 
the group T-values for the alcohols, aldehydes, aromatic hydrocarbons, esters, ethers, hydrocarbons, 
inorganic acids, ketones, nitrogen oxides, organic acids, and miscellaneous groups. These calculations 
are conducted according to the following equations: 


^group — E C c /C m , 

(Bl) 

^overall = £ ^group > 

(B2) 
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where C c is the contaminant concentration in the atmosphere in mg/m 3 and C,„ is the maximum allow- 
able concentration in the atmosphere in mg/m 3 . 

The criteria for acceptability are the following: 

1. The T-value for each group must be less than one 

2. The overall T-value must be less than one. 

If either of these criteria are exceeded, the atmosphere is considered unacceptable. 
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